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Summary. — It is shown that the form of momentum spectrum of jet 
particles is essentially dependent on the degree of anisotropy of the angular 
distribution, and that e.g. the form dp/p? corresponds to a definite 
«medium » anisotropy. Integral spectra of directly measured values of p 
indicate a systematic decrease of transversal momenta in the vicinity of 
the axis and of the z/2-plane in C.M. system. 


The momentum spectrum of particles emitted in nuclear interactions with 
primary energies above 5-10! eV has been investigated in nuclear emulsions 
by direct measurements of the momenta p in individual events (e.g. (1*)) 
as well as indirectly using the angular distribution of large collections of jets (*). 
As a generally accepted result the differential momentum spectrum in C.M. 
system is of the Heisenberg form dp/p?. 
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We should like to point out that the form of the momentum spectrum de- 
pends strongly on the degree of anisotropy of the angular distribution and we 
shall try to show that the above result does not hold generally. 

At first, let us suppose that the transversal component of the momentum p, 
is constant and independent on the primary energy as well as on the angle of 
emission 09. Then it follows from the definition p = p, cosec 0 that any form 
of the momentum spectrum is connected with one definite type of the angular 
distribution; e.g. for the Heisenberg type spectrum we have in the C.M. system 
the angular distribution cotg 0 d(cos 0) representing a «medium » anisotropy. 


dp d(cosec @) 
p cosec? 0) 


(1) dn ~ = cotg 0 d(cos 6) . 


Expressing the angular anisotropy in the more usual way as cos” 6 d(cos@) 
we have a spectrum 


DI (m—1)/2 IL 
(2) dn ~ cos” 0 d(cos 0) ~ {1 — — = OD 


DE D 
going over to da — dp/p* for large values of p. For a «medium» anisotropy 
about m = 4 the distribution eq. (2) is similar to that of the formula (1). * 
To investigate the influence of the angular distribution on the form of the 
momentum spectrum we divided 21 jets with primary energy above 5-101! eV 
from our laboratory and several others from the literature (4°) into three 
groups according to their anisotropy expressed by the value of y= <cosec 0),,= 
= 1.254/m (for large values of m) ace. to (4) and (9): 


(I) y<1.4 isotropy > = 13 1227 tracks 
(II) IO medium anisotropy <y>,, = 2.3 116 tracks 
(III) 4< y= 16 strong anisotropy <y>,,=5.5 40 tracks. 


In Fig. 1 we plotted the integral spectra of the values (in the C.M. system) 
of cosec 0 = const. p for individual groups. The group I of isotropic jets 
follows fairly well the curve corresponding to m = 0 in eq. (2), the spectrum 
of the medium anisotropic events (group II) is very near to the form dp/p*, 
while the sequence of points related to the strongly anisotropic group IIT 
decreases much more slowly. The hard tails of all the spectra seem to have 
a parallel slope approximately as dp/p°°. 
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In Fig. 2 we collected together the data of all events from group I, 1I and IIT. 
The average value of y of all these jets is “Pay = 2.0 corresponding to a 
«medium » anisotropy (roughly m = 4 in eq. (2)); we were therefore not sur- 
prised that we obtained a rather good agreement with the line dp/p? (4). 


1%al : 2 \ \ Tey 
1 10 cobec © 10 10 cosec@ 100 
Fig. 1. — Integral momentum spectra of three Fig. 2. — Integral momentum spectrum 


groups of jets with different anisotropy. of all jets from Fig. 1 together. 


It is shown on Figs. 1 and 2 that the dashed line drawn for the fraction 
F—0.5 crosses the spectrum of any group of jets in the point cosec d= 
= (cosec 9>, == y. Consequently the spectra of isotropic events lie to the left 
of the line dp/p? while the strong anisotropic jets have the momentum spectrum 


shifted towards higher values of p. 


Similar argumentation holds also 100% _ 0.5cosec © 
for the directly measured spectra. E.g. 
the Torino jet (1) with the experimen- 50% 
tal momentum spectrum dp/p? has a À | 
medium anisotropy y = 1.8 which is | N 
very near to our group II or to the | 
collection of our Fig. 2. An eventual ol ‘sil 
difference between the integral spectra ie 
DE dp \ 
| BA x 
Fig. 3. — Integral spectrum of directly : 
measured values of p (points) and of mo- \ 
menta defined as p=0.5 cosec 0 (histo- 1% VA _— 
gram); average anisotropy y=2. 01 ay ie p(Gev/c) 10 
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of cosec 0 and of the measured p values might come from some kind of deviation 
from our assumption p, = const (7°). In Fig. 3 we compared the sequence of 
points related to the experimental 
values of p of four jets measured 
in this laboratory (the mean 
y= 2) with the histogram of their 
0.5 cosec 0 GeV/c values. While 
the agreement is satisfactory be- 
tween 7 = 0.7 and 0.3 there seems 
to be a shift of the experimental 
points to the left in the regions 
corresponding to the vicinity of 
the x/2-plane and of the axis, in- 
dicating there an experimental 
decrease of the p, values. 

In Fig. 4 we therefore plotted 
the histogram of the mean values 
(Lee ai de O) of log p, of the same jets (full line) 

0° 20° 40° 60° 80° 100° 120° 140° 160° 180° and also together with the S-star 

Fig. 4. — Angular dependence of transversal and the Torino jet (dashed line). 
momentum of a group of jets from Fig. 3; There appears to be a decrease of 
points: individual values; full lines: mean p, values in the neighbourhood of 
values; dashed lines: mean values of lg Pr to- the axis and of the z/2-plane as it 
gether with jets (1) and (5). was already indicated on the inte- 

gral momentum spectrum (Fig. 3). 

The full curves in Figs. 3 and 4 are drawn as predictions (7) of the two- 
centres model (*) of particle emission. Adhering tentatively to this model we 
tried to guess the form of the energy spectrum in the rest-system of the emitting 
centres O*. The hard tails of the momentum spectra are nearly insensitive to 
transformation into another system; we may deduce from their parallel slopes 
in the C.M. system (Fig. 1) thatin the C*-system the momentum spectrum near 
to the axis is also of the form 


(3) <P with a>2 


(see Figs. 1 and 3). Because of the isotropy of the angular distribution in the 
C* system (°) and having in mind the approximate dp/p* law for the isotropic 
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(9) J. G. MCEVEN: preprint (we are greatful to prof. GIrERULA for kindly sending 
us the preprint). 


2340 


ON THE MOMENTUM SPECTRUM OF PARTICLES IN COSMIC RAY JETS 133 


group I (Fig. 1) we may deduce that a similar spectrum as (3), is valid in the 
whole space of the C*-system and that the interaction (proportional to x) 
among created particles is very strong in this emitting centre. The resem- 
blance of eq. (3) with the laboratory spectrum of photon-electron cascades 
with primary energy above 10! eV (19-11) seems to be interesting. 


(10) J. G. DUTHIE, C. M. FISHER, P. H. FowLER, A. KADDOURA, D. H. PERKINS 
and K. PINKAU: Moscow Conference on Oosmic Rays (1959). 

(11) Y. Fupyimoro, S. HAsEGAWA, M. Kazuno, J. NisHimura and K. Niu: Moscow 
Conference on Cosmic Rays (1959). 


RIASSUNTO (*) 


Si mostra che la forma dello spettro della quantità di moto dei getti è essenzial- 
mente dipendente dal grado di anisotropia delle distribuzioni angolari, e che, p. es., 
la forma dp/p? corrisponde ad una definita anisotropia del mezzo. Spettri integrali 
di valori di p misurati direttamente indicano una diminuizione sistematica dei momenti 
trasversali nella vicinanza dell’asse e del piano 7/2 nel sistema del c.m. 


(*) Traduzione a cura della Redazione. 
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Summary. — 117 interactions of cosmic ray particles with aluminium 
nuclei have been studied. 101 showers containing no apparent unstable 
particles yielded the following multiplicities: 


n,> = 7.33 ny y= 3.0% (RICO CR 10:08 


The momentum and angular distributions of the secondary particles were 
examined, and an attempt was made to study the momentum-angle 
correlation. The mean transverse momentum was 0.31 GeV/c, the RMS 
deviation about this mean being 0.23 GeV/c. A quantity related to the 
(U — p,) of the recoil nucleon was also investigated. 12 showers containing 
one V° each and 1 shower containing two V° particles were studied; this 
showed that the 117 showers contained: 29.541 A° and 18.0*150 0° par- 
ticles. 3 showers containing V= particles were also studied; one of them 
was probably a E particle. 


1. — Introduction. 


Two groups of physicists from the Imperial College (London) and the Ecole 
Polytechnique (Paris) took part in this experiment. It was carried out with 
the equipment of the Ecole Polytechnique at the Pic du Midi de Bigorre 
(alt. 2860 m). This equipment, originally used for work on V° particles, has 
a maximum detectable momentum of about 30 GeV/c, and is therefore not 


unsuited to the study of nuclear interactions with primary momenta around 
40 GeV/c. 
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An aluminium producer was chosen, as it satisfied the following conditions: 


a) its low atomic weight favoured the occurrence of a high proportion 
of nucleon-nucleon interactions; 


b) at the altitude of operation the producer could be made thin enough 
(4 em) for the loss of wide-angle particles to be small, the detection efficiency 
of unstable particles to be relatively high, and yet the frequency of occurrence 
of well measurable interactions to be sufficient for a statistical study. 


2. — Experimental procedure. 


21. Apparatus. — The installation at the Pie du Midi was described in detail 
in a previous article, ref. (1), We mention here its essential characteristics 
and any modifications made during the present experiment. 

Two Wilson cloud 


chambers were situated os 

one above the other as il- RE = 6 
lustrated in Fig. 1. The 340m 

upper chamber, of useful | ZZ 0000 FF, C 
volume (68X68 X30) em*, M ia 
was situated in a mag- 

netic field of about | 

5000 gauss whilst the lower Ì 3 | | 
chamber, of similar useful RS 

volume, contained 13 cop- | | 

per plates each 1 em thick. i | | 

A lead slab 22 cm in thick- === A Ei ig D 
ness separated the two 2 S| ZZ. 
chambers. The aluminium | = È 
producer, 4 em thick, GETZ. ee esi 
formed the roof of the | e 
magnet chamber. The | ease See 
expansion of the cham- e 

bers was initiated by a 

coincidence pulse from the i 

counter trays B, €, D : PE 

and Æ. Two combina- L ii 
tions BCD,E,—B, and | 

BCD,E,— B, where used. Fig. 1. 


È Y 
(1) B. P. GRÉGORY, A. LAGARRIGUE, L. LEPRINCOE-RINGUET, F. MurLLeR and C. 
Peyrou: Nuovo Cimento, 11, 292 (1954). 
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Each chamber was photographed by 3 cameras, and a supplementary 
camera placed on a level with the underside of the producer improved visi- 
bility in the important region near the origin of the interaction. 


2°2. Measurement. — The three magnet chamber photographs of each event 
were projected in turn at unit magnification and the curvature of each track 
measured by comparison with standard curves. 

The maximum detectable momentum (M.D.M.) of the chamber was deter- 
mined by measuring 145 no field u-meson tracks which passed through the 
lead slab and the multiplate chamber. From the distribution of the measured 
sagittae of these u-mesons tracks we have obtained a value of M.D.M. of 
30 GeV/c for a track length of 40 em. 

The reconstruction of the tracks in space was obtained from full scale 
drawings of the three projections. The relevant calculations using these measu- 
rements were then undertaken with the aid of an electronic computer. 


3. — Selection of events. 


The events analysed were selected from the non-electronic showers using 
the following criteria: 


i) Multiplicity: the number of minimum ionizing tracks with mo- 
mentum > 110 MeV/c, n,, must be > 3. 


li) Energy: the total resultant momentum of the charged secondaries, 
Yp,, must be >10 GeV/c. 


ili) Geometry: the events must have good geometry, that is: 
— be visible on all three photographs, and 
— have > 60% of all tracks measurable over 25 em of track length. 


iv) Origin: we have rejected all events which had two or more particles 
incompatible with the origin defined from the drawings of the three photo- 
graphs. 

In 9600 photographs, we found 142 showers satisfying the first three cri- 
teria, as follows: 


84 without V° 


Tan fr 8 with V° 
( nique origin ii 
142 dpr 10 GeV 3 with V~ 
CARI 17 without V° 
1 without 
ood geometry 21 3 ti À 
g g À one incompatible HN 


om f 16 shower origins detectable 
| 9 multiple origin 
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We have therefore chosen for study 101 showers apparently unaccompanied 
by unstable particles, and 16 showers with unstable particles. 

In addition we have three showers of energies certainly greater than 
100 GeV, which we have been unable to study because of superposition of the 
tracks. 


4. — Showers with no apparent unstable particle. 


41. Introduction. — Among the 117 showers selected as described in Sec- 
tion 3 there were 101 without apparent unstable particles. 
In this Section, we shall study the properties of the 101 showers under the 
following headings: 
— the primary particles and the nature of the interaction; 
— the secondary particles: nature, multiplicity, distributions of mo- 
menta, angles and transverse momenta; 
— the quantity 1.5 } (U—p,). 


42. Primary particles. 


42.1. Mean momentum. — 
The primary momentum was defined DI 
from the vector sum of the mo- 
menta of the charged secondary 
particles (> p,). Fig. 2 shows the 
distribution of > p, for the 101 
showers. The mean and median 
values are about 28 and 27 GeV/c, 
respectively. 

Assuming that there is charge in- Fig. 2. — Total visible momenta of the 
dependence in the production of se- showers. 
condary particles, and that in each 
shower the primary particle becomes the most energetic of the visible parti- 


number of showers 
5 


Ne i] mal 
1020 30 40 50 60 70 80 90 Ey, Ge 


cles, we can write: 


pi D Pz + 3 x (pr —D5) ; 


where p is the mean primary momentum, 
Sp, is the mean value of the vector sum of all the charged 
secondaries. 
Sp, —p,) is the mean difference between Sy and the momentum 
of the most energetic visible particle. 
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Among our showers, 
> 2; == 23 Gen ic, 


> pc 
D 30, CON Ce 


This is compatible with the mean energy of primary protons between 10 
and 100 GeV, obtained from ref. (?). 


42.2. Nature of the interaction. — For our experiment we adopted 
a triggering system designed to eliminate neutral primary particles (cf. Sec- 
tion 2:1). At the altitude of the Pie du Midi, the interacting charged com- 
ponent of cosmic rays contains protons, x+-mesons, x-particles. The propor- 
tion of «particles at our altitude is negligible (ref. (?2)). 

From ref. (4), the ratio 


number of charged 7-mesons 
number of protons 


e 0.4 ZE 0.1 . 


From the ratio given above, we have calculated the numbers of our primary 
particles which may have been protons and z+-mesons. Assuming that these 
two types of particles have the same interaction cross-section, we have found 
the following proportions: 


70% of interactions were caused by protons, 


and 


30% by charged x-mesons. 


We have attempted to use the distribution of the mean excesses of posi- 
tive particles among our showers (Fig. 3a) in order to obtain some information 
about the nature of the interaction of the primaries with the nucleus of alu- 
minium. For each shower, the positive excess is the number of positive par- 


ticles minus the number of negative particles. Electrons identified as in 
Section 43 (below) are ignored. 


(7) 5. F. SINGER: Progr. Cosmic Ray Phys., 4, 282 (1958). 
©) H. V. Newer: Ann. Rev. Nucl. Sci., 8, 225 (1958). 
(4) K. GREISEN: Progr. Cosmic Ray Phys., 3, 108 (1956). 
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The mean positive excess is 1.7, after account has been taken of the par- 
ticles oî unknown sign. The mean positive excess depends not only on the 
mean positive excess in the primaries but also on the mean number of target 


* one particle of indeterminate sign 
xx (wo particles ” » 


® number of showers 


1 


x XxX X x X x xX 


15] 


By ene Pod Pas oS CE EN ERIN DE 


107 


ISAIA EY ET BEES Gente 


a 6 2 


Positive excess per shower 


Fig. 3a. — Excess of positive particles. 


nucleons involved in the collision. If there are equal numbers of positive and 
negative 7-meson primaries then the mean positive excess in the primaries 
for the quoted ratio of z+-mesons to protons (~ 0.4), is 0.7. Since there are 
14 neutrons and 13 protons in the aluminium nucleus we should therefore have 
an average of two nucleons involved in the collision in order that the struck 
protons may yield the residual positive excess of unity. If the nuclear radius 
is given by 
where 

iit, = od Mr” Gan 
and 


A — atomic weight, 


a calculation based on a simple model of the nucleus (ref. (5)) shows that a 
single particle travelling through the nucleus of aluminium .will make an aver- 


(5) B. Rossi: High Energy Particles (1952), p. 360. 
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— observed distribution ss ie allowance 
--.theorelical distribution if all primaries positive | made for 


; ; 85% positive and 15% |+ secondaries 
negative primaries 


® number of shower 
I] 
I 


3 


—.— È 2 = 
22 (0) 2 4 6 8 positive excess per shower 
Fig. 3b. — Excess of positive particles: primary-nucleon collisions. 

age of 1.8 collisions, provided the nucleon-nucleon cross-section is taken as 


20 mb. Since this is close to the experimental values of the nucleon-nucleon 
cross-section at 10 GeV, ref. (f), we have calculated the expected distribution 


si bam 

è| — observed distribution 

<= 

<4 Lasa theoretical distribution if85% positive and 15%} Allowance made 
x negative primaries for secondaries 
Q 

[= 

2 

30 


205 


10 


positive excess per shower 


Fig. 3c. — Excess of positive particles: primary-Al nucleus collisions. 


(0) V. S. BARASHENKOV, V. A. BELYAKOv, E. G. BUBELEV, WANG SHOU FENG, 
W. M. MaLrsEv, TEN Gyn and K. D. Torsrov: Nucl. Phys., 9, 74 (1958). 
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of the positive excess for the simple case of a single particle traversing an 
aluminium nucleus. The broadening effect due to the particles of unknown 
sign was taken into account and the nucleus was assumed of uniform density. 
Fig. 3b and 3c show the results. | 

Fig. 3b demonstrates the expected distribution if the primaries consisted 
of 30% x*-mesons and 70% protons and collided with a single nucleon in 
the aluminium nucleus (dashed curve). For comparison we have shown the 
distribution we would have expected if the primaries had consisted solely of 
protons (dot-dash curve). 

Fig. 3¢ then gives the expected distribution when the probabilities of col- 
liding with several nucleons in the aluminium nucleus is taken into account. 
The disagreement between the expected distribution and the observed distri- 
bution is not large. This may indicate that there is a mechanism which pre- 
vents the m-meson secondaries produced in the collision from interacting again 
in their parent nucleus, for otherwise we might expect a considerably larger 
positive excess than observed. 


43. Nature of the secondary particles. — These may consist of protons, 
charged hyperons and K-mesons, charged 7-mesons, and electrons. 

At momenta below a few GeV/c, the ionizing effect of the particles increases 
as the momentum falls. The ratio of track ionization to minimum ionization 
is >2 for protons and hyperons of momenta < 750 MeV/c, K-mesons 
< 350 MeV/c, and r+-mesons < 110 MeV/c. We therefore estimated that we 
could identify baryons, K-mesons, and 7*-mesons, below the corresponding 
momenta given above. All particles below 110 MeV/c with minimum ioni- 
zation have been considered electrons. 

In addition we have thought it possible to identify among the minimum 
ionizing particles with momenta over 110 MeV/c a number of electron pairs 
from considerations of geometry and scattering. 

This method yielded 45 pairs of electrons. Assuming that the momentum 
spectra of the negative shower particles and the 7°mesons were similar, we esti- 
mated that in addition to these 45 pairs with both electrons over 110 MeV/c, there 
were about 20 electron pairs with one electron below and one electron above 
110 MeV/c. This means that about 20 single electrons are probably included 
in our shower particles, since we have found it impossible to associate them 
with their low-energy partners because of the large angle through which the 
latter were scattered. 

In our study of the secondary particles in the paragraphs which follow, 
we have used all particles which appear to come from the shower origins, less 
the identified electrons and heavily ionizing particles with momenta >110 MeV/c. 
Separate histograms will be drawn for these latter particles. 
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Hence the particles we have studied and which are shown in our various 
histograms are 7+-mesons, possibly contaminated with: 


— charged nucleons and hyperons, of momenta > 750 MeV/e; 
— charged K-mesons, of momenta > 350 MeV/c; 
— electrons which may have esca ped detection, of momenta > 110 MeV/e 


According to various observers, the proportion of heavy particles (nucleons, 
K-mesons, hyperons) among charged shower particles may be of the order 
Ofe20 97,4 (Reta) 


= 


4°4. Multiplicity. — In Fig. 4 we 
have plotted the distribution of the 
multiplicity of our showers. This 
multiplicity, n,, was defined as the 
number of particles described at the 
end of Section 43. The mean and 


fe c | o median multiplicities are 7.3 and 7.0 
3 6 9 2 5 #8 A 2% 2 ns per shower. 


on 
=A 


number of showers 
© 


Fig. 4. — Particle multiplicities of the The mean multiplicities are given 
showers. in the following table: 
mean 7, = KM) = 7.28 particles per shower 
mean multiplicity of positives in = 3.93 
mean multiplicity of negatives =) = 2070 
mean multiplicity of particles of indeterminate 
sign = (ny = 0.64 
mean multiplicity (including heavily ionizing 
particles) = ‘Mg + n> = 7.63 
mean multiplicity of positives (including heavily 
ionizing particles) — (ni n> = 4,28 


If we distribute the n> between wi and n; in the proportion ni :n7, i.e. 
1.44:1, we shall have: 


<n, =, — 4.66. ne 2.907 


In order to estimate the number of particles produced in a collision of a 


primary with an aluminium nucleus in our energy region we have proceeded 
as follows. 


() B. Epwarps, J.;Losty, D. H. PERKINs, K. PINKAU and J. ReyNoLDS: Phil. 
Mag., 3, 237 (1958). 

(8) L. F. Hansen and W. B. FRETTER: preprint. 

(?) A. I. ALIKHANOV, G. P. ELiseev, V. $. KAMALJAN, B. A. Liusimov, B. I. 
Morseev i A. B. KHRIMJAN: Zurn. Eksp. Teor. Fiz., 36, 404 (1959) 
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Since z+-mesons and protons formed the major fraction of incident pri- 
maries and on the average a neutron and a proton in the aluminium nucleus 
were probably involved in the collision (see Section 42) we have assumed 
that the following reactions took place: 


n= + (p +n) > 0.33 r= + 0.33 27 


PIEDE SD) 0.07 p= 0-5 . 1 


where the first two occurred 
30% of the time, the third 70% 
of the time, and it is assumed 
that charge independence holds. 
Hence the charge remaining 
on the degraded primary and 
the target particles after colli- 
sion was 1.45 positive and 0.10 
negative per interaction. Since 
4.46 positive and 2.96 negative 
secondaries were observed per 
interaction then 3.21 positive 
and 2.86 negative were really 
produced. 

The distribution of », as 
a function of > p, is shown 
in Fig. 5. There is a tendency 
for high energy events to have 
a larger number of secondary 
particles than low energy events. 


Fig. 5. — Momentum-multiplicity 
correlation of the showers. 


+ 0.33 7°4 (p+ n) + produced particles 


+ (p + n) + produced particles , 


30: 


Xp, Gev/c 


20 30 40 50 60 7080 100 


in 


45. The momenta of the secondary particles. — Fig. 6 shows the distribution 
of the momenta (p) of all shower particles (as defined at the end of Section 4°3) 
and Fig. 7a those of the positive and negative particles. There are also 
32 M.D.M. particles with p> 20 GeV/c which are not shown in these histo- 
grams. For the calculation of primary momenta (Section 42.1) and axis 
directions (Section 4°6) we have assigned a value of 20 GeV/c to these particles. 
(Ideally, all particles with measured momenta approaching the M.D.M. should 
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be corrected for the effect of simulation which, because of the rapidly falling 
secondary momentum spectrum, makes it possible for a given particle to si- 
mulate one of higher or lower momentum depending on whether the diffe- 
rential secondary momentum spectrum has a slope greater ot less than — 2. 
But in our case the slope is too badly defined in the required momentum 


region.) 


number of particles 


1 DIRE) 5 10 20 p Cave 


Fig. 6. — Particle momenta. 
For 1< p< 15 GeV/c, the distribution in Fig. 6 is given by 
= = pee with 1<q<1.6, 


where V,= > n,. The positive and negative particles, taken separately, also 
seem to be distributed according to the same law. 

There is some difference between our results and those of ref. (9): this may 
be due to the different range of primary energies in our case. 
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The ratio of the numbers of positive and negative particles (n° /n,)as a 

. . . . . - . À . ù a 
function of momentum is given in Fig. 7b. This diagram has a maximum at 
0.5<p<1.0 GeV/c, which is probably due to recoil protons with 0.75 <p < 
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Figs. 7a, b. — Distribution of particle momenta. 


<1.0 GeV/c. (Recoil protons with p< 0.75 GeV/c, identified by the ioni- 
zation of their tracks, are not shown in the figure.) There appears to be 
another increase in the positive excess at about 3 GeV/c. 


46. Angular distribution of secondary particles. — In the absence of more 
precise information, we have taken the direction, of the resultant momentum 
of all the secondary particles as an estimate of the primary direction. The 
error due to fluctuation in the number of z°-mesons is negligible, while that 
due to fluctuation in their angular distribution is (cf. Appendix I): 


0.44/2 
Ag = rado 
Beye Ms 


For example, the axis of a shower of 7 charged particles with > p, of 
28 GeV/c, is defined with an error, Ag, of about 2° (assuming a nucleon-nucleon, 


interaction.) 
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We have studied the detection efficiency in the illuminated region of the 
chamber for all the showers. We found that: 


all particles with 0 < 40° 
two-thirds of particles with 40°< 6 < 60° 
one-third of particles with 60°< 9 < 90° 


would be detected. 0 is the angle at the origin between the shower axis and 
the particle’s track. 

Fig. 8a shows the angular distribution of positive and negative particles. 
The median angle is 8°. There is a noticeable positive excess below 5°. Angular 
distributions of particles of indeterminate sign (including M.D.M. particles) 
and of heavily ionising particles are given in Fig. 8b and &c. 
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à 
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1 
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Fig. 8a, b,c. — Angular distribution of particles. 


The angular distribution of all our minimum ionizing particles is similar 
to that in ref. (°). 
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47. Momentum-angle correlation. 


HTC. 147 


— We examined our showers in an at- 


tempt to discover whether there was any correlation between the angles of 
emission and the momenta of the particles. Fig. 9 shows the results. 


ber of particles 
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on 
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Fig. 9. - Momentum-angle correlation. 


As one would expect, the mean momentum of the particles decreases as 
the angle of emission inereases. There is an excess of positive particles at most 
angles as we saw in Section 46, but this excess is unusually great between 0° 
and 5°. Tt is difficult to explain this observation in terms of either 
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a) recoil protons which are unlikely to be found preferentially at 0<5° or 


b) a simple mechanism in which the primary particles (mainly protons) 
lose only a small fraction of their energy, since the ratio n, /n, is greatest at 
low particle momenta. 


Particle momenta = ( 0= 5) GeV/c ni ln, = 2.2 
25210) = ME 
20515) St hg 


48. Transverse momenta of secondary particles. — The transverse momentum 
(p,) of a secondary particle produced in a high energy interaction is an im- 
portant experimental quantity since it is invariant under a transformation 
from the laboratory to the centre of mass system of the interaction. The dis- 
tribution of p, observed in the laboratory system should therefore represent 
the true distribution in the C.M.S. and as such may be readily compared 
with theory. 

Since in the laboratory system the particles emitted in the C.M.S. backward 
hemisphere make relatively larger angles with the primary direction than do 
the particles emitted in the ©.M.S. forward hemisphere, one might expect 
the p, values of the former to be biased to slightly higher values than those 
of the latter. In order to test this hypothesis, each interaction, was transformed 
to an effective C.M.S. in which the resultant momentum of the charged par- 
ticles was zero. The p, distributions of forward and backward particles were 
found to be similar. We surmise that the error in the primary direction 
(cf. Section 46) has little effect on the distribution of p,. 

The p, distributions of positive and negative secondaries (defined at the 
end of Section 4°3) have been normalized and superposed in Fig. 10, which 


—— positive particles 


--—- negative particles normalised to 
the same area 


number of particles 


is 
14 16 pr Gev/c 


Fig. 10. — Transverse momenta. 


indicates that there is no significant difference between them. The mean value 
of p, is 0.31 GeV/c, and the R.M.S. deviation is 0.23 GeV/c. From the simi- 
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larity of the two distributions we infer that since particles of nucleonic mass 
must be present in reasonable number they must either 


i) have a p, distribution not very different from that of the x*-meson 
secondaries (if they are present predominantly as particles of the same sign, 
e.g. protons); or 


ii) exist in approximately equal numbers with positive and negative 
charges (if their p, distribution is appreciably different from that of the 
T*-Mesons). 


As shown in Section 4°4, it seems that there were very few charged hype- 
rons among our showers. This indicates that the particles of nueleonic mass 
were in a large majority protons. The similarity of the positive and negative p, 
distributions is therefore good evidence that the p, distribution of the protons 
is not markedly different from that of the 7*-mesons. We note in passing that 
the p, of the heavily ionizing particles has a mean value of 0.21 GeV/c, with 
a R.M.S. deviation of 0.12 GeV/c. 

Fig. 11 is a histogram of the p, values of shower particles of both signs. 
Two curves have been superposed on it for purposes of comparison. The dotted 
curve is a theoretical p, distribution taken from ref. (!°), for a final thermal 
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Fig. 11. — Transverse momenta. 


equilibrium temperature of 7,= m,c?/2k. The other curve is derived from the 
experimental momentum distribution of the 7*-mesons created in antiproton 


annihilation events, ref. (11). 


(19) G. A. Miceëin and I. L. RozenTAL: Suppl. Nuovo Cimento, 8, 770 (1958). 
(11) E. SEGRÈ: Report of Kiev Conference, (1959). 
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Fig. 12a. — The correlation of the median transverse momentum of the showers and 
their multiplicity. 
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The fact that we were able to measure the momenta of most particles of 
each shower allows us to consider each interaction separately. Fig. 12a shows 
the median p, value of each shower as a function of its multiplicity. These 
quantities are independent. For each value of », the values of median Pi are 
widely dispersed. This diagram may be of theoretical interest if p, really 
reflects the C.M.S. momentum of the particle. This would be true if the ++- 
mesons were emitted isotropically from their centre of production. We also 
found that the median p, is independent of the total visible momentum, dp, 
(Fig. 120). 


49. The quantity 1.5 Y.(U—p,). — We have thought it if interest to con- 
struct the distribution of the quantity 1.5 Y (U—p,) for the minimum ioniz- 
ing particles of each shower (Fig. 13), where U is the laboratory energy of 


lone ina two showers 
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Fig. 13. — 1.5(U — p;) for all minimum ionizing particles; all particles taken as 7-mesons. 


each shower particle, assumed of pion rest mass, and p, the longitudinal com- 
ponent of particle momentum along the primary direction. The factor 1.5 
is introduced to take account of the unseen r°-mesons. In Fig. 14 the same 
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quantity has been plotted except for those events where the most energetic 
particle had a momentum greater than twice that of any other particle in 
the shower. In these cases the fastest particle was assigned the mass of the 
nucleon and its (U—p,) was not multiplied by the factor 1.5. This pro- 
cedure was followed since it was thought not unlikely that this particle was 
the degraded primary nucleon. 


Ve en a al 
CRC Eee A EE 15(L-9)6eV 


[area for each shower 


ee ORE: 


RATTO eds ST OL GeV 


Fig. 14. —1.5(U — p,) for all minimum ionizing particles; all particles taken as -mesons 
(except that any particle with momentum at least twice that of any other in its shower, 
is taken as nucleon). 


In these figures the upper diagram represents the distribution of 1.5 
> (U—p,) where each event is indicated by a vertical line. In order to bring 
out the concentration of these events along the (U—p,) axis each event has 
been assigned the same area and the lower diagram has been constructed. 
Both figures show evidence of structure rather than a continuous distribution. 
We note that the first three maxima of Fig. 14 are situated at approximately m_, 
2m, and 3m,. The distribution of 1.5 Y (U—p,) below 0.42 GeV (~ 3m.) 
is of special interest since any recoil proton associated with these events 1. 
likely to be identified by its ionization and therefore not included in the sum- 
mation 1.5. Y (U—p,). The reason for this is that there exists, for nueleon- 
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nucleon or z-nucleon collisions, a direct correspondence between the quantity 


1.5 Y(U—p,) and the (U—p,) of the recoil nucleon given by the equation 
(cf. Appendix IT): 


1.5 > (O— p,) + (U— Py) seco = M, = Nucleon mass . 


Thus for events with 1.5 Y (U—p,)< 0.42 the (U—p,). is greater than 
0.52 GeV. Reference to Fig. 15, which shows the relation between the mo- 
mentum and angle of the recoil proton for various values of its (U—p,), then 
indicates that recoil protons with (U—p,)>0.52 GeV and transverse mo- 
menta in the region below 0.3 GeV/c will have momenta < 750 MeV/e and 
therefore be detected by ionization. 
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The apparent structure in the distribution of 1.5 > (U— p,) may thus re- 
flect the fact that the (U—p,) of the recoil nucleon is restricted to certain 
values. . | | 

In Fig. 16 a plot of the multiplicity of each event against its 1.5 > (U—p,) 
shows that a correlation exists between these two quantities. 


2361 


207 a baie 
®*N=28 
è 151 : 
(©) 
î (AE 
= fee 
q 10] 5 . * 
2 AE IR 
ely a ; 
DO“ ; ; : 
5 4 7 
F } | 2 x Ns 
5 7 9 11 13 15 17 19 21 


Fig. 16. — The correlation of 1.5 } (U—p,) and multiplicity. 


5. — Showers with unstable particles. 


51. Introduction. — Among the 117 showers selected as described in Sec- 
tion 3, we observed: 


12 showers each associated with the production of one neutral unstable 
particle ; 
1 shower associated with two neutral unstable particles; 
3 showers each associated with a charged unstable particle. 


Of our neutral unstable particles 12 had both decay products with measu- 
rable momenta. 

In the remainder of this Section we shall refer to neutral unstable par- 
ticles as V° particles, and to charged unstable as V+ particles. 

In this Section, we shall study: 


— the V° particle alone; 

— the showers with V° particles, after the V° particles have been incor- 
portated into their respective showers; 

— three showers and their associated V+ particles. 


52. The V° particles. 


52.1. Compatibility of the V° particle with the shower ori- 
gin. — We tested the vector sum of the two decay products of each V° to see 
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whether the vector direction passed through the shower origin. By taking 
into account the errors in 


— the position of the shower origin, 

— the position of the decay point of the V°, 

— the direction of the decay products, we found that the measured de- 
viations were always smaller than the calculated values. This means that all 
our V° particles can be considered as: 


i) originating at the shower origin, 


li) decaying into two particles only. 


52.2. Efficiency of detection of V° particles. — An attempt was 
made to estimate the number of V° particles which escaped detection. For 
this purpose we defined a « visible region » in the chamber, believing that any 
V particle decaying in this region was bound to be observed. 

For each of our observed V° particles we measured the distance it would 
have travelled inside the visible region if it had not decayed on the way. We 
called this distance Z. We then calculated the probability that it should have 
decayed within the region, 4.e. 


P_exp| A exp| il: 


where: À is the mean distance travelled before decay by a V° of the same 
energy, 


l’ is the distance between the shower origin and the V° particles 
entry into the visible region, and 


RTE), 


For each V° particle we have calculated P on both 0) and A° hypotheses, 
assuming that the mean lifetimes of 0° and A° particles are 10-1° s and 3-10- s, 
respectively; ref. (2). This yielded: 


A° hypothesis CLP yee) 
0° hypothesis CIRO 


The R.M.S. deviations of the distributions of 1/P are 0.5 and 0.7, respectively. 

Our detection efficiencies were therefore similar for A° and for 0° particles, 
if we suppose that both particles are created with the same momentum and 
angular distributions. 


(2) C. Franzinertr and G. Morpurco: Suppl. Nuovo Cimento, 6, 469 (1957). 
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52.3. Momenta of the V° particles. — The distribution of the momenta of 
the V° particles is given in Fig. 17a, which shows that 8 out of the 12 V° particles 
have momenta greater than 3 GeV/c. In addition to the 12 V° particles in- 
cluded in the figure, there were two whose decay products were visible along 
10 cm only and whose momenta could not be measured. 
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Fig. 17. — a) Momenta of particles from showers with V° particles. 


52.4. Armenteros-Podolanski diagram. — Tf a V° particle decays to give 
two charged secondaries, it may be represented by a point on. this diagram, 
ref. (1°). The position of the point is defined by the transverse momentum (p,) 
of the secondaries, and by the «-value of the V°, where 


2 2 
Pi — PÎ 
C= E 


ne 
Po; P+, P_ are the momenta of the V° and its positive and negative secondary 
particles, respectively. The points representing different V° particles lie on 
families of ellipses. There is one family for each decay scheme of the V°, and 
one member ellipse of each family for each value of po. If p, is greater than 
about 2 GeV/c, then for each decay scheme the ellipses become indistinguishable. 

Hence, for p, > 2 GeV/c, we can draw a single ellipse (shown on the dia- 


(5) J. PODOLANSKI and R. ARMENTEROS: Phil. Mag., 45, 13 (1954). 


2364 


AN EXPERIMENT ON NUCLEAR INTERACTIONS ETC. 157 


gram, Fig. 18) for V° particles (A° particles) with the decay scheme: 


A°>p+x-+37.9 MeV, Ref. (14). 


-09 -06 -03 0 03 06 09 a 


Fig. 18. — Armenteros-Podolanski diagram. 


All A° particles have a positive value of «. For p,> 2 GeV/c, we can also 
draw a single ellipse (Fig. 18) for V° particles (0, particles) with the decay 
scheme 

0° > + +-+ 217 MeV, Ref. (14). 


The values of « of 0} particles have equal a priori probabilities of being 
either positive or negative. 

The values of « and p, are plotted in Fig. 18. The errors given are calcu- 
lated from the extreme values of the angle between the secondaries, and of 
their momenta. This figure shows that all 12 V° particles have positive values 


(14) C. p’ANDLAU, R. ARMENTEROS, A. AsTIER, H. C. DESTAEBLER, B. P. GREGORY, 
CG. LeprINcE-RincueT, F. Mutter, C. Peyrou and J. H. TINLOT: Nuovo Cimento, 
62135 (1957)- 
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of x. Since the «-values of 0° particles are equally likely to be positive or 
negative, this result inclines us to suppose that a high proportion of our ob- 
served V° particles are A°. 

This conclusion is borne out by the distribution of 6* for A° and 6) hypo- 
theses, shown in Fig. 19. 0* is the angle in the rest mass system or the sve 
particle between the positive decay product’s direction and the line of flight 
of the V°. 
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Fig. 19. — Distribution of 0*. 


We may note that of the 5 V° particles with showers of Sp, below 


10 GeV/e which occur among our 9600 photographs, three have « positive 
and two have « negative. 


5.2.5. Number of A° and 0° actually produced. — In order to calculate these 
numbers we had to know two things: the numbers of A° and of 6° particles 
among our 14 V° particles, and the proportions of the total numbers produced 
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(allowing for the undetected V° particles) which these numbers represented. 

Of the 12 V° particles on the Armenteros-Podolanski diagram, 3 are com- 
patible with 0% and lie in a region where the probability of finding a A° is nil. 
The other 9 are in a region of p,= 50 to 100 MeV/c and «> 0. If we sup- 
pose that decay products of a V° particle are on the average emitted isotro- 
pically in the rest mass system of the V°, a given A° has a probability of 86%, 
of lying here, and a given 0! has a probability of 8%. Hence in a given batch 
of V° particles we would have expected about 10.5 times as many A° as 0? par- 
ticles in this region. This would indicate that there were about 8 A° and 
4 0% among these 12 V° particles. Each of the 2 unmeasurable V° particles was 
thus twice as likely to be a A° than a 0° particle. 

We therefore supposed that 10 + 1 of our V° particles were A° and 47 1 
were 0°. 

Using the detection efficiency determined in Section 52.2 and the fol- 
lowing branching ratios, ref. (1): 


te SETT (32 +5)% 
A 0; (50%) artt (854 6)% 

0° 0° 
“Of (50%) Nate (15 4 6)% 


we estimated the total numbers of A° and 6° particles produced in our 
117 showers. We found: 


29.5+140 A° particles, 
9.0+75 6° particles, or 18*%, 6° particles. 
The branching ratios quoted in ref. (1°) are slightly different and would 


have given us mean numbers of particles produced of 31.5 for A° and 19.0 


lora0°: 

It is interesting to note that if the 4 V° particles interpreted as 0! were 
indeed 0° particles, then the probability of observing and « — p, distribution 
similar to theirs was only 6%. 


513. Showers with V° particles incorporated. 


53.1. Momenta of the primary particles. — The total visible mo- 
menta (> p,) of our 11 showers with measurable V° particles are shown in 


(5) Report of Annual International Conference on High Energy Physics (CERN, 1958)- 
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Fig. 20a and 20b. The median value is 24 GeV/c. The presence of the neutral 
particle has a notable effect on the distribution of > p,. 
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Fig. 20. — Total visible momenta of showers with V° particles. a) Including V° particles. 


b) Excluding V° particles. 


53.2. Multiplicity, momenta, angular distribution and trans- 
verse momenta of the particles. — 


a) The mean multiplicity of charged particles (other than heavily ioniz- 
ing particles and identified electrons) is 6.9 particles per shower. ‘This is si- 
milar to the value obtained for the showers without apparent unstable par- 
ticles (cf. Section 4°4). 


b) The momenta of all the measurable particles are shown in Fig. 17). 
It will be seen that below 10 GeV/c the charged shower particles are distri- 
buted roughly according to the formula given in Section 45. The momenta of 
the V° particles show no such distribution. 

The relatively small number of secondary charged particles with sufficient 
momentum to produce V®s of the observed momenta indicates that, unless 
their cross-section for producing V° particles on their way out of the alumi- 
nium nucleus was high, then the V° particles must have been produced directly 
in the collision of the primary particle with a target nucleon. 


c) We have defined the primary directions by the vector sum of the 
momenta of all particles, including the V° particles. The angular distribution 
of V° particles is shown in Fig. 21; all are emitted with an angle < 20° with 
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the direction defined above. The angular distribution of charged particles is 
compatible with that of the particles from showers without V° (cf. Section 4°8). 


number of particles 


St Tec raso Se ESE 
30° 40° 50° 60° 0 


Fig. 21. — Angular distribution of V° particles. 


d) The mean transverse momentum of the V° particles (Fig. 22a) is 
0.59 GeV/e with an R.M.S. deviation of 0.15 GeV/c, and of the other particles 
0.31 GeV/c with an R.M.S. deviation of 0.16 GeV/c. There is no apparent 


a) 

44 

N 
24 K 
M F i 

0 CI MLA Mal onal. 
3 | 
3 
rs 
© 
2) 
& 
D 
& 


__- normalised transverse momenta 
of particles from showers without Ves 


OAIO 208104 05106 0/7 FOLNOS TO Mins 2 pi Gev/c 


Fig. 22. — a) Transverse momenta of V° particles. b) Transverse momenta of charged 
particles from showers with V° particles. 


difference between, the transverse momenta of positive and negative particles. 
In Fig. 22b, the transverse momentum distribution of the charged particles 
from the V° showers is superposed on the normalized distribution for the non- 
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V° showers. The 7? test indicated that the two distributions are significantly 
different. The main contribution to y? comes from the excess of particles with 
p, values in the region of 0.5 GeV/c in the case of the V° showers. 


53.3. Positive excess. — The mean positive excess of the 11 showers 
with V° particles of measurable momenta is 2.5 particles per shower. If we 
add the two V° showers with unmeasurable V®s the mean excess does not 
change. This is significantly greater than the mean positive excess of the 
101 showers without apparent V° particles (cf. Section 5°2.2), and may indi- 
cate that on the average more nucleons of the aluminium nucleus are involved 
in interactions in which V®s are produced than in others. 


54. Oharged unstable particles (V*). 


541. Experimental data. — We observed three showers with charged 
unstable particles among our 117 showers. Their chief characteristics are given, 
in the Table: 


| Momentum of. Momentum of 
; 7 T+ 
Nature | Momentum of V= charged sec- eni 


| \ondary (GeV/e)) secondary 


Angle of 
decay of V= 


Wee p > MDM along 25 cm too short 2 ~1° 
Ve too short 0.61 0.02 ? 6.1° 
Wi | too short 0.2 +0.01 2182203 TIRI” 


As the mean lifetimes of K* and charged hyperons are of the order of 10-8 
and 10-1° s, respectively, the Ve particles may be considered hyperons under 
our experimental conditions. 

The third particle (V~) is accompanied by a V° particle which decays in 
the chamber about 1 em from the point of decay of the V7. The vector sum 
of the momenta of the charged secondaries of the V° coincides with the line 
passing through the point of decay of the VT and the shower origin. On the 
Armenteros-Podolanski diagram, this V° seems compatible with both A° and 6?. 
It lies, however, in a part of the diagram where a A° has a higher probability 
of occurring than a 0°: The best explanation of this event seems to be: 


Pao 


54.2. Number of charged hyperons. — It is difficult to calculate 
the detection efficiency for charged hyperons as it.is difficult to detect small 
angles of decay, especially 
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i) if the particle decayed near the top or bottom of the visible region 
of the chamber; 


ii) if the angle of decay were in the plane perpendicular to the front 
glass of the chamber.. 


(However, one of our three charged hyperons was detected with a decay 
angle of about 1°.) Consideration of the data from V®s suggests that 
the mean decay angle of Y* particles is of the order of 10°. Taken in con- 
junction with what is known of the lifetimes of charged hyperons, ref. (15) 
and (!*), this makes it reasonable to suppose that our detection, efficiency for 
charged hyperons was not very different from that for As of similar energy. 
It therefore seems that in our interactions charged hyperons were produced 
in significantly smaller numbers than A° particles. 


* * * 


We would like to thank Dr. J. ROScH, director of the Observatoire du Pie 
du Midi, and his staff, who made this experiment possible; Dr. R. DE PossEL, 
director of the Institut Blaise Pascal, who kindly placed à Bull electronie 
computer at our disposal; and Dr. B. W. POWELL. 


APPENDIX I 


The error in the estimated primary direction due to non-observation of 7°-meson 
secondaries. 


The variance in p, for isotropic emission is given by 


pi 5 |e" sin? 0* d (cos 0*) = 


+ 


€ 
dl 


p* = 0.M.S. momentum, 0* = angle in C.M.S. with primary direction. 


The projected value of p, is given by 0.82V/p*?/2. If then n, charged 
mesons are produced the deviation of the resultant transverse momentum 


from the true value due to the unseen 7°mesons is 0.41V: n.p*?, The corre- 


(15) P. H. FoWLER et al.: quoted by A. ALVAREZ in Report of Kiev Conference (1959). 


+ 
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sponding angular deviation of the estimated primary direction from the true 
primary direction, in one plane, is then given by 


0.41V np 0.41 np” 


Dp:cos By > E* 


Ns 


tg Ag = 


Putting B,=1 and p*? = E* we have the total error in the primary direction: 


0.412 


RD Gray 


AToDia Wim xe ITAL 


By energy conservation we have: 
Up+ Ur = a U, 


where U>, is the energy of the primary particle, 
U, is the energy of the target particle, 
YU is the total energy of the secondary particles in the laboratory 
system after the collision. 


By momentum conservation we have: 


Prtpr= > p cos0 = > px, 


where pp is the primary momentum, 
Pr is the Fermi momentum of the target particle, 


p, is the total longitudinal component of momentum of the sec- 
ondary particles along the primary direction. 


=" 


(Op— Pr) +> Un tpr= 
=1.5 > (U—p,)+(U—p,) degraded primary + (U— p,) recoil . 


If the recoil nucleon is either neutral and unseen or heavily ionizing we have: 
1.5 D (U —p,) + (U— px) degraded primary = (U — p,) recoil + My. 
Since for a 30 GeV primary nucleon 


(Ur — pr) =0.01GeV, pr=0, U, = mass of target = M, . 
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RIASSUNTO (*) 


Si sono studiate 117 interazioni di particelle da raggi cosmici con nuclei di alluminio. 
101 sciami non contenenti in apparenza particelle non stabili d'edero le seguenti 
molteplicità: 


‘n= 7.3; nt>=39; n7)= 2.7; nt)= 0.6. 


Si sono esaminati gli impulsi e le distribuzioni angolari delle particelle secondarie, e 
si è fatto un tentativo di studiare la correlazione tra impulso ed angolo. L’impulso 
trasversale medio è 0.31 GeV/c; la deviazione RMS su questa media è 0.23 GeV/c. 
Si è anche studiata una quantità riferentesi all’(U —p,) del nucleone di rinculo. Si 
sono studiati 12 sciami contenenti un V° ciascuno ed uno sciame contenenti due parti- 
celle V°; ciò ha dimostrato che i 117 sciami contenevano: 29.519 particelle A° e 
18.0*1%% particelle 0°. Si sono studiati anche sciami contenenti particelle V~; una di 
esse era probabilmente una particella E7 


(*) Traduzione a cura della Redazione. 
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Summary. —- A cosmic ray experiment to study 100 GeV nuclear inter- 
actions with a magnet cloud chamber is reported. The analysis of the 
interactions is described and the hypotheses necessary for interpreting 
the measurements are discussed. Results from a first sample of 41 inter- 
actions are given. They show that 1) the transverse momenta of the 
secondary particles are generally low (< 0.5 GeV/c) and 2) in the centre 
of momentum frame of reference the emitted particles mostly have low 
momenta and their angular distribution is anisotropic. There is also 
some evidence that in a nucleon-nuclear interaction at these energies 
the number of nucleons involved is small. The proportion of strange 
particles among the secondaries is similar to that found in interactions 
of much lower energy. 


Introduction. 


This paper reports the first results of an experimental study of the nuclear 
interactions of cosmic rays with energies of about 100 GeV. The interactions 
occurred above the sensitive volume of a magnet cloud chamber situated at 
an altitude 3500m. No direct observations of the primary particles were 
made. 

The experiment to be described was essentially exploratory. The appa- 
ratus was not specifically designed for the experiment and was, in some ways, 
unsuitable. One important aim of the work was to find, from experience, the 
difficulties that arise in analysing 100 GeV interactions and to see what appa- 
ratus would be required if a further experiment were envisaged. 
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In the following sections we describe the way in which the cloud chamber 
photographs were selected for analysis, the methods of analysis, and the re- 
sults obtained from a sample of 41 high energy interactions. 


1. — Apparatus, selection of photographs, yield. 


The cloud chamber and magnet installed at the Jungfraujoch have been 
described previously (1). The useful volume of the cloud chamber is 
(5>x 55x 16) cm* and is 
placed in a magnetic field 
of 5000 gauss. The max- 
imum detectable momentum 
for tracks longer than 30 cm 
was about 10 GeV/c. 

A vertical section through 
the apparatus is shown in 
Fig. 1. The wedge-shaped 
volume at the top of the 
cloud chamber that is not 
«seen» by the cameras was 
filed by a block of paraffin 


BEE 

wax of mean thickness 6 cm. 
This paraffin target was 
chosen to give a high pro- 
bability of observing a sin- 

; : x 
gle nuclear interaction free 
from complication due to O A Ae AA © O 
N è © È SS À j F GI] MY b FT A 
secondary processes (either O SO: 
nuclear interactions or ma- AN LEE 
terialisation of high-energy > DRE (0) 
y-rays from 7°-decay). The 0 25 cm 


« geometrical » interaction 


: x Fig. 1. — Vertical section through the a aratus. 
length in paraffin is 63 cm s a Le 


1) Magnet coils. 2) Iron pole pieces. 3) Paraffin 
and the cascade length is target. 4) Stereoscopic camera. 5) Lead. A, B) Counter 
52 cm. Moreover, one third trays containing 8 counters 2 cmx40 em. €, D) 
of the geometrical cross- Counter trays containing 8 counters 3 em x 60 em 


section is due to simple hy- separated by 0.5 cm lead. 


drogen nuclei. 
The Geiger counters used to control the expansion of the cloud chamber 
are also shown in Fig. 1. A number of different counter combinations were 


(1) J. A. NEWTH: Suppl. Nuovo Cimento, 11, 297 (1954). 
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a 


168 L. MONTANET, J. A. NEWTH, G. PETRUCCI, R. A. SALMERON and A. ZICHICHI 


initially tested for rather brief periods and, finally, the simultaneous discharges 
necessary to expand the chamber were selected to be one or more from tray A 
and three or more from tray B and two or more from each of trays C and D 
(4,,+B,,+C.,+D-;,). With this combination, the counting rate of the 
selection system was 7.6 per hour and, with a dead time of 4 minutes after 
each photograph, the rate of photography was 4.9 per hour. 

The photographs were scanned by projecting to about full size and those 
showing high-energy nnelear interactions apparently occurring in the paraffin 
target were selected for careful measurement. The wide-angle photographic 
system makes it difficult to locate the site of an interaction by simple scanning; 
after measuring the stereoscopic pair of photographs it was frequently found 
that the interaction, in fact, occurred in the brass roof of the cloud chamber 
or even above the roof. 

Initially, the criteria used to select interactions for analysis were rather loose 
but, after some experience the following strict criteria were adopted: 


1) The photograph should show at least three copunctual tracks of fast 
particles (i.e. particles with momentum greater than 1 GeV/c) of which one 
had a momentum greater than 10 GeV/c or three had momenta greater than 
5 GeV/c or five had momenta greater than 1 GeV/c. 


2) There should be no evidence for a second nuclear interaction in or 
near the target. 


3) The interaction should not be explicable as pure electron-photon 
cascade; this meant that the angles between the trajectories should not all be 
compatible with zero. 


4) The majority of the tracks should lie in the well-illuminated volume 
of the cloud chamber. 


5) The photograph should not show evidence for anomalously large dis- 
tortion of the tracks. 


Interactions that satisfied these criteria were termed «class A » interactions 
and a summary of the yield of such interactions over a period of three months 


TABLE I. — Yield of interesting interactions. 
: ; : Number of 
Duration of Number of Sr: Interactions Interactions 
experiment photographs RE with n7> 15 with m7< 5 
85 days 9153 188 19 40 


The table shows the resulis of careful scanning, with uniform selection criteria, of the main 
body of photographs. Class « A » interactions are those with good geomeiry in the illuminated 
region of the cloud chamber and where there is no confusion with particles not directly asso- 
ciated with the interaction. 
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is shown in Table I. The distribution of total multiplicity (n,) for these inter- 
actions is shown in Fig. 2. 


D 
oO 


WwW 
(©) 
le 


188 events 


Number of particles 
N 
(=) 
À 


i +8 events 
1 with ny >20 
ol ni Sl 
OT 10 15 20 
Multiplicity (nr) 
Fig. 2. — The multiplicity distribution for 188 class « A» interactions selected from 


9.153 photographs. There is a cut-off imposed by the scanning criteria that eliminates 
multiplicities below 3. The small numbers of interactions with n7—4 or 3 is probably 
due to the discrimination of counter tray B against low multiplicities. 


When the total multiplicity is greater than about fifteen the measurements 
on the tracks become inaccurate because of overlapping and confusion. For 
very low multiplicities (n,< 5) the analysis becomes inaccurate owing to the 
difficulty of determining the axis of the interaction. In our experiment we 
concluded that only the interactions with », between five and fifteen could 
profitably be analysed. From Table I it is clear that about 1.5 such inter- 
actions are recorded per day. Of these, about one half occur in the paraffin 
target; the remainder are nuclear interactions in the roof of the cloud chamber 
or other neighbouring material. This yield of interactions is obtained with 
an «acceptance » of about 40 cm?-sr for the primary particles. 


2. — Measurements and analysis in the laboratory system. 


About 50 interactions (*) were selected for measurement. The photographs 
were projected through a photographic enlarger on to a horizontal table and 
the curvatures of the tracks were measured by comparison with standard 


engraved circular arcs. Drawings of the two stereoscopic views were then 


(*) Most of these were selected before the criteria given in Section 1 were estab- 
lished. Some of the 50 were later rejected for various reasons. An analysis of the 
remaining interactions is being made (?). 

(2) L. MonTANET: to be published. 
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made and the relevant tracks were extended to find the co-ordinates of the 
nuclear interaction. The direction cosines of the tangent to each track at 
the interaction were then calculated. 

In order to define an axis for each interaction, the vector sum of the mo- 
menta of all the charged particles having momenta above 1 GeV/c was found. 
(In the case of a particle whose momentum was too high to measure the ap- 
propriate maximum detectable momentum was used in the calculation.) 

The whole interaction was then referred to a system of spherical co-ordinates 
with the axis as pole. Every track was specified by its co-latitude (0) and 
its longitude (y). A target diagram of each interaction was then drawn on 
polar graph paper. 

The above operations were all routine in character. Before proceeding with 
the analysis of the interactions an attempt was made to exclude the tracks 
due to electron pairs, to knock-on electrons and to evaporation protons. Low- 
energy electrons were rejected by excluding any particle with minimum ioniz- 
ation and momentum less than 100 MeV/c. Protons that could be identified 
from ionization and momentum (< 500 MeV/c) were also excluded. Electron 
pairs were more difficult to recognize and no really satisfactory identification 
was made. However, on a target diagram an electron pair that is not close to 
the axis of an interaction stands out as a pair of tracks with opposite curva- 
tures and with a negligible included angle. Ten such pairs were observed and 
rejected while a rough calculation based on the number of 7°-mesons expected 
in the interactions predicted that fourteen electron pairs should be produced. 
It seems clear that the majority of the pairs were identified but, equally, that 
a small number probably remained among the particles used in the later ana- 
lysis. 

To summarize, at the conclusion of this stage of the analysis each inter- 
action was specified by its origin (P in the paraffin target, M in the doubtful 
region near the roof of the chamber) and by the multiplicity (,) of fast charged 
secondary particles. Each fast secondary particle had a momentum (p;) and 
angular co-ordinates (0,, g,) referred to the axis of the interaction. 


3. — Determination of the C-system. 


The next step in the analysis was the determination of a Lorentz trans- 
formation along the axis of each interaction such that, after transformation, 
the interaction possessed some specified forward-backward Symmetry. We call 
the frame of reference in which this symmetry exists the O-system and its 

x] reale ‘ava € ne a a 2 —4 . . . 
velocity relative to the laboratory is Be, {y = (1— B?)-*. The physical signi- 
ficance of the C-system and of y, is clear when the interaction is between like 
particles as in the collision of an incident proton with a hydrogen nucleus 
in the paraffin target and, among our Interactions, some are probably of this 
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type. When the interactions are more complex the significance, if any, of y 
depends upon the phenomenological model that is used to describe the inter- 
actions. Here we describe only the procedure for finding y, and defer the dis- 
cussion of its meaning until later (Section 8). 

The main problem in determining y, is to choose from an infinity of possi- 
bilities some criterion of symmetry that must be satisfied in the C-system. 
Tf all the momenta in the laboratory system are accurately measured the choice 
of y, can be made so that the total momentum of the charged particles in the 
C-system is zero or 

Ns 
(yes DD cos Gu 0, 

it 
where p* and 0* are the momentum and co-latitude of a particle after trans- 
formation. In our experiment this criterion of symmetry is not usable because 
there are, in nearly every interaction, one or more particles with momenta 
to high to measure. 

We adopted the criterion 


ns 


Ss if 
cos 0* = — Y cos? = 0, 
Ns 1 


— 
bo 
— 


as our condition for symmetry. It has several properties that are worth noting. 
First, the value of cos 0* is a smoothly varying function of y and the value y, 
that satisfies condition (2) can be found graphically by interpolation after 
making only two or three trial transformations. 

Second, there are cireumstances in which the application of criterion (1) 
would lead to results that are physically quite unreasonable and where cri- 
terion (2) would be more satisfactory. For example, if a neutral particle is 
emitted with a high momentum from an interaction the application of (1) to 
the charged particles only may lead to an unrealistic value of y,: For reason- 
ably high multiplicities condition (2) is less affected by the presence of 
high-energy neutral secondaries. 

Third, condition (2) has the advantage that the statistical error in y, may 
easily be found for any assumed angular distribution of the secondary particles 
in the C-system. Specifically, if the angular distribution is of the form 
N(0*) d0*xccos"0* sin 0* d0* the average value of cos 0* will be zero with 


an R.M.S. deviation of V(m+1)/(n.(m+3)) so that the condition (2) be- 
comes 


<="? mHt 1 
9 FR 
(2a) cos 6 0 + be 73) 
or 
(2b) = 0 + (E for m — 0 (isotropy) 
Ns 
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hes 
(26) =0 +} — for large m . 
Tes 


In its practical application, the procedure was as follows. A guess at the 
value of y, was made and the Lorentz transformation was calculated using 
the measured values of p; and 0; and assuming all the fast charged particles 
to be z-mesons. The value of cos 0* for the interaction was found. This was 
repeated for two or three further values of y and the results were plotted in 
the form of a graph of cos 0* against y. The intersection of the graph with 
the axis of y gave the value y, and the values of y corresponding to cos 0* = 
— + V1/3n, gave limits to y. at one standard deviation. 

Fig. 3, curve 1, shows a typical result for interaction XS 29 (Plate I) where 


the value of y, derived is 9.3*3} and these numbers give some idea of the 


08 
RA ARE IT, 


0.64 
J Event XS 29 
=c0s 0* el 
n Ose Fema paenc 


Fig. 3. — Diagram showing the determination of 7. and of the anisotropy in an inter- 
action. In a) are shown the variation of cos 0* and |cos 0*| as functions of the y of 
the Lorentz transformation for event XS 29 (see Table II and plate 1). In b) the 
same variations are shown for an artificial event with y= 10; consisting of 10 r-meson 
isotropically distributed and each having a p* of 350 MeV/c. The points to notice 
in these curves are the large errors in y, and the slow variation of |cos 6*| with y. 
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TABLE II. — Details of 41 analysed interactions. 


i 

No. Origin | np My Ne n_ ty | YetAyo | AEAA > La ; 
x1) 826 P 10 0 3 3 4 DA 0.04 4012, 2.50 
XO 826 P 7 0 0 0 7 ISS 0.475220, 14a red 
XO 693 M 9 0 3 3 3 5.625) 1 0.69 220,13 1.05 
KO 857 M 8 0 4 2 2 4.6712 | 0.55 40.16 | 2.22 
STE] M 6 0 1 1 4 20 059010210) 092 
Om 74 P 6 0 3 2 1 6.6200 1 070.7 lass 
XO 968 M 8 0 2 3 3 Gee a 0.68 0-01 2046180 
XO 1229 M 6 0 2 2 2 5.9534) 0,66 20.18 NL 
XO 1568 M 4 0 2 3 2. 1 10 2750.66) 0.19) 2079 
XO 1696 2 8 0 3 2 3 8.973. | 0.805+0.16 | 1.01 
XO 1928 M 4) 0 3 1 3° 110.3252 10:56 0.16 2075 
XP 505 2 5 0 1 0 4 5.625% | 0.51) 20.16) 1.08 
XP 506 M 6 1 3 0 2 10.722 0575 420.38) 0.52 
NON P 12 0 5 1 6 5.8525 00.88 020.1400250 
NOA? L 7 0 5 2 N22 10:63 0:13 0:02 
XQ 185 M 5 0 2 0 SSA 03320240 Mo 
XR 109 P 8 1 3 1 3 4.5232 (| 0.76" 2201610200 
DAS OOH P 9 1 4 2 2 6,255 1.0.80 LO 40 1904 
REI M 6 0 2 0 4 ASSO 4000 5a ULI: 
XR 880 È 5 0 2 1 PMN ees AE ER TES ALAN 
XR 971 M 7 1 2 2 2 4460 0.62 40.33 | 0.94 
XR 1019 M 6 0 1 1 PI SE 2 0.62 20.32 0.37 
XS 29 10 9 0 4 1 4 93253 00,65. +02] 1.08 
XS 190 P 7 0 3 0 Ae ee) 0.49) £0.37 Ord? 
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TABLE II (continued). 


No. | Origin Nm ma Na Nn No VetAye A+AA eT 
—|- | 
XS 250 | M 11 0 3 3 5 5.27% | 0.65+0.14 | 2.62 
FU 7 | M 8 1 3 2 2 552 | 0.56+0.18 | 1.55 
UL Sal Moni 11 0 6 2 3 042210068201 MI 
GY UF pe ae 10 1 3 2 Pee ie; +35 0.48+0.14 | 0.82 
RTE 18 1 2 5 10 | 12.5435 | 0.5640.09 | 1.48 
| 

XT 466| P 7 1 0 0 6 | 23 +2” | 0.82+0.26 | 0.27 
XT sai | M 11 0 3 2 6 9.5723 | 0.48+0.12 1.30 
XT 1071 | M 6 0 3 1 2 39.6558.) 0.40 SO SE PO ES 
XT 1186 | M 7 0 3 3 1. | 1535757 150. 2.0.14 


XT 1785 | M i 0 8 2 7 14.5732 1) 0.650 120126 
XT 2053 | M 11 0 2 6 3 4.0 22+) 0722 0.16) 3.06 
XT 2156 | P 8 1 4 2 1 5.6722 | (06432018152 
SRL M 20 1 9 5 5 <3.5 = = 
XS id | M 7 1 2 3 1 2 Ba OOO == 
XS 408 | mM | 90 | 1 | 4 | 6 | 9 | eet? toeo0 comme 
XT 1627 | P 10 4 2 1 3 227% = = 
XR 1004 | P 208 OMIS LE = 


This table lists the details of 41 interactions that had ns > 5 and could be analysed by the 
method of Fig. 3 to give a value of y,. The origin P is in the paraffin target, M is doubtful- 
copper, aiuminium or pareffin. A is the mean vajue of |cos 0*| for the tracks from each inter- 
action after applying the Lorentz transformation Vou 

The errors in y, and A were calculated assuming that the angular distribution in the C,M.S. 
was anisotropic with a value 0.65 for A. The mean value of A found, in fact is 0.615-+ 0.025. 
This means that the errors have been over-estimated. Moreover, the error in 4 due to the 
error in y, has been taken as symmetrical whereas, in fact, there is a slight bias that tends 
to underestimate the value of A (see Fig. 2). Therefore the 4A are generally rather larger than 
they should be and the overall mean value (4=0.61) is probably slightly too low. 

The five interactions below the interrupted line are excluded from the later analysis because, 
for the first four there are very good reasons to think that they are very complex interactions 


and for event XR 1004 the resolution of the separate tracks on the photographs does not allow 
accurate measurement. i 
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PLATE 1. — Event XS 29. This event is listed in Table II and the determination of » 
is shown in Fig. 3. 
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PLATE 2. — Event XT 280. This event is listed in Table II and Table IIT. The tracks 
A and B come from a neutral V-decay while 0 and D are probably due to an electron 
pair produced in the gas of the cloud chamber. 
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errors in the determination when #,= 9. Curve 3 in Fig. 3 shows «a similar 
determination for an artificial event that was constructed by assuming ten 
secondary 7-mesons to be emitted isotropically with a constant p* in a system 
for which y, was 10 (*). 

The method of finding y. outlined above was used on 41 interactions with 
ns>5 listed in Table II. For the first 36 events the results were straight- 
forward. For four events the deduced values of y, were extremely low. It was 
impossible to interpret these four events as nucleon-nucleon collisions since the 
visible energy of the interaction was greater than that corresponding to a N-N 
collision with a C-system velocity of B, (E. = (2y;—1)m,). For the last 
interaction the high multiplicity made measurements difficult and a reliable 
value of y. could not be obtained. 

The spectrum of the 36 values of y, is shown in Fig. 4. 


36 events 


LIMI WE et LI 


2 gs 


tC a 


Fig. 4. — The distribution of values of y, found by the criterion DI cos 0*—9 for the 
first 36 interactions listed in Table II. For a simple nucleon-nucleon interaction 
the laboratory energy is related to y, by the equation Fay = (27 — 1) Myc. Thus a 
value of 100 GeV for Ex would correspond to y,=7.3 and 200 GeV to y,=10.4. 


4. — The distribution of transverse momenta. 


The 36 interactions for which a value of y. could be found without diffi- 
culty contain 273 charged secondary particles. Of these 102 were positive, 
64 negative and 107 had momenta too great for their signs to be determined. 

Before considering the distributions of momenta and angles of emission 
in the C-system we describe the distribution of transverse momenta 
(ps = p;sin0; = pî sin 0°). This is invariant under the transformation from 
L-system to C-system and is unaffected by errors in y,. 

For the 166 particles whose momenta were measured the distribution of p, 
is given in-Fig. 5. This shows clearly that the vast majority of the measured 
particles (that is, at least 50% of all the charged secondaries) have transverse 


(*) In Fig. 3 the errors have been found assuming an isotropic distribution of the 
secondary particles. The observed anisotropy in our events led us to assume a value 
m=1 in (2a). The error in y, for event XS 29 is then increased to the figure in Table II. 
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momenta less than 400 MeV/c. The number of particles with measured values 
of p, greater than 400 MeV/c is very small. 

Since only 60% of the charged secondaries are measured, it is important 
to consider whether the remaining 40% of particles could all have very dif- 
ferent values of p.. In 
general, the unmeasured mo- 
50 menta belong to particles 

in two classes; those giving 

long straight tracks at very 
166 particles small angles to the axis and 
those at large angles to the 
axis where the visible track 
is very short. Those wide- 
angle secondaries whose 
tracks lie nearly in the plane 


n 
oe 


Number of particles 


oO 


D 0202000 IT IO 
Pr (GeV/c) —+ measurable and Fig. 5 in- 
Fig. 5. — The distribution of transverse momenta cludes values from many 


for 166 measured secondary particles. This distri- such tracks. It is unreason- 
bution is biased since, in the same interactions there able to imagine that the 
are 107 particles whose momenta were too high to wide-angle secondaries giving 

measure. This bias is discussed in the text. 


short tracks have different 
momenta from those giving 
long tracks and we can therefore ignore the bias due to the exclusion of wide- 
angle secondaries from the distribution of Fig. 5. 

The p,-values from the high-momentum (> 10 GeV/c) particles inclined at 
small angles to the axis are a more serious difficulty since quite small changes 
in the assumed axis direction give large changes in p,. The choice of axis 
—described in Section 2—essentially leads to small values of p, (*). Without 
observing the direction of the primary particle it is very difficult to assess the 
error involved in our method of choosing the axis. What we have done is to 
investigate, in particular cases, the effect of altering the momenta assumed for 
very energetic particles. The resulting change in the axis direction is generally 
small (~ 2°) and the corresponding alteration in the values of p, is of the 
order of 200 MeV/c for a particle with a momentum of 5 GeV/c. 

Our conclusion from this discussion is that at least half of the secondary 
charged particles have transverse momenta less than 400 MeV/c. We have 
no reason to think that our results are strongly biased by the exclusion of the 
particles whose momenta could not be measured. We therefore expect about 


() Considering the p,-values as vectors in the plane of the «target diagram » it 
is clear that our choice of axis makes YP: =0 and makes the sum XV pi a nies 


2384 


A CLOUD CHAMBER STUDY OF NUCLEAR INTERACTIONS ETC, ity 


80% of all the secondary particles to have p, below 400 MeV/c. The errors 
in the determination of the axis direction make it impossible to say anything 
about the distribution of the p,-values below 400 MeV/c. 


5. — The momentum spectra in the C-system. 


The momentum of a particle in the C-system is obtained from the labo- 
ratory measurements by the Lorentz transformation 


(3) pi = D; sm (Lee: (cos he 5 JE 


where £., y, are the constants of the transformation and £, is the particle’s 
velocity in the laboratory system. The above transformation was carried out 
for all the 166 measured secondary particles assuming, for the calculation of 
BP:, that every particle was a 7-meson. 

The resulting momentum spectra for positive and negative secondary par- 
ticles are shown in Fig. 6 (a and db). Both spectra show that the majority of 
the measured secondaries have p* less than 1 GeV/c; only 2 out of 64 nega- 
tive particles and 9 out 102 positive particles have momenta above this limit. 

The criterion of measurability used in selecting the particles whose momenta 
are plotted in Fig. 6 may well lead to a serious bias in the distribution. This 
is because both the classes of unmeasured particles referred to in Section 4 
may consist of particles whose momenta in the C-system were high. For very 
high energy particles near to the axis of the interaction (cos 0; =,= 1) 
equation (3) reduces to p; = p;/2y. and, with our values of y,, we obtain 
large values of p* when 7p; is greater than 20 or 30 GeV/c. For the wide-angle 
particles of low momentum the quantity (cos 0;— f./f:) in (3) is negative and 
can easily be large enough to give a high value for p*. In fact, of the 11 mea- 
sured particles with high values of p*, 8 are slow, wide-angle particles of 
this type. 

As far as the measured particles are concerned, it is clear from (3) that the 
large error (— 25%) in y. will be transferred nearly directly to give similar 
errors in the values of p*. Moreover, even if the hypothesis that all the nega- 
tive particles are m-mesons is reasonable, among the positive particles there 
are certainly many protons. The errors in p* due to a wrong assignment of 
mass can be found from (3). 

The only difference in the two momentum spectra of Fig. 6 is the suggestion 
of a larger proportion of positive particles with p* above 1 GeV/c. If the nine 
positive particles with high p* are assumed to be protons and f; is reduced 


12 - 7] Nuove Cimento. 


2385 


178 I. MONTANET, J. A. NEWTH, G. PETRUCCI, R. A. SALMERON and A. ZICHICHI 


accordingly, it is found that, in eight cases, the altered values of p* are even 
greater than the original values. 
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Fig. 6. — The values of p* for 64 negative particles (a) and 102 positive particles (b) 

from the first 36 interactions of Table II. These values have been calculated using 

the best values of y, and on the assumption that all the particles are z-mesons. The 

error due to this assumption and the effect of the 107 unmeasurable particles are 
discussed in the text. 


From the above discussion we conclude first that more than 50% of the 
charged secondary particles have momenta in the O-system of less than 1 GeV/e. 
Secondly, we attribute large errors to the momenta shown in the spectra of 
Fig. 6 and fear that the unmeasured secondaries could alter their Shapes sig- 
nificantly. Finally, we observe some indication of a positive excess among 
the particles with high p*; if these particles are protons, their momenta in the 
C-system are very high (> 5 GeV/c). 
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6. — The angular distribution in the C-system. 


The angle of emission (0*) of a particle in the C-system is related to the 
quantities we have determined by 


4 to 62 — sin 0; 
63. g oi = ve (cos 0;— B.JB:) 


This involves the laboratory momentum and the mass of the secondary par- 
ticle only through the value of f;. Since p; does not appear explicitly, it is 
possible to find values of 0* even for unmeasured secondaries by assuming 
that f, lies between unity and some lower limit set by the maximum detectable 
momentum. In nearly all cases the uncertainty in 0* is only a few degrees 
although in rare cases it can be 30° or more. 

The values of 0* calculated from (4) were, of course, used to determine y, 
from the condition Y cos 0% — 0. When this condition is satisfied some mea- 
sure of forward-backward symmetry of the angular distribution is demanded. 
The anisotropy of the distribution can be conveniently specified by the quan- 
tity 1/n, Y|cos0}|= A. The expected value of A for a distribution propor- 
tional to cos” 6* is: 


m + jee m +1 
4ns(m re yp 


which, for isotropic emission (m = 0) reduces to 
(6) A=0.5+1/V12n, 


In our experiment the value of A was calculated for each event during 
the determination of y.. The results are given in Table II and the distri 
bution of values of A is shown in Fig. 8. The statistical errors are so large 
that the value of A obtained from a single interaction rarely differs significantly 
from 0.5 but the distribution of Fig. 8 shows a clear tendency for A to be 
greater than 0.5 (the mean value is 0.61, corresponding to m = 1). Further 
evidence for the existence of anysotropy is show in Fig. 7 where the angular 
distribution of all 273 secondary particles is shown. 

There are two sources of error in the determination of A that should be 
considered. The first is the error due to the large uncertainty in y, In Fig. 3 
(curves 2 and 4) the value of A is plotted as a function of y and itis clear that 
the variation is extremely slow. There is, moreover, a systematic tendency 
for A to have a minimum value when y= y,: Thus there is a small syste- 
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matie error which tends to underestimate the anisotropy of the angular distri- 
bution. 


273 particles 
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Fig. 7. — The angular distribution of the secondary particles in the first 36 inter 


actions of Table II. There is no significant forward-backward asymmetry but the 
anisotropy of the distribution is clear. 


The second source of error is important for interactions of low multiplicity 
and is associated with the determination of the axis direction. It can most 
easily be visualized by imagining an interaction with only two fast secondary 


36 events 


— 
(©) 
1 


Number of events 
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Fig. 8. — The distribution of the anisotropy A (=|cos 0*|) for the first 36 showers of 
Table IT. A=0.5 is expected for an isotropic distribution. 


particles (p> 10 GeV/c) which are inclined at an angle 20 to one another. A 
strict application of our analysis procedure to such an interaction would lead 


to the results: y, = cosec0, and A—0. This arises from our inability to 
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measure momenta above 10 GeV/c — the axis is therefore chosen to make 
the angles 0; equal. Any other choice of axis would lead to a smaller value 
of y, and, of course, a larger value of A. It is difficult to estimate the size of 
this bias in interactions with higher multiplicities but we excluded all inter- 
actions with multiplicity less than 5 from'our analysis. The residual bias is 
probably small and, in any case, it again leads to an underestimate of the 
anisotropy. 

Taking into account the statistical error in A (using m — 1 in equa- 
tion (5)) and the error arising from the uncertainty in y, we obtain for all our 
interactions a weighted mean value of A=0.61+0.025. We consider this to 
be proof of an anisotropic angular distribution in the C-system. 


7. The production of strange particles. 


Among the 36 well analysed interactions of Table II there are only two 
(XT 123 and XT 280) that show the decay of a strange particle. To obtain 
better statistical information about the frequency of production of strange 
particles, all the 188 interaction (Table I) were searched for V-particle decays. 
A further nine were found and the interactions were measured. The results 
for the eleven interactions are given in Table III. 


Taste III. — Interactions in which V-particles are observed. 
. Nature of 
Event no. Origin ge Tr Ny 
XT 123 M wo 11 0 
XT 232 M VE | 11 1 
XT 280 E not A° 10 1 
XW 64 IR vic 11 0 
XW 660 1P AE In | 0 | 
XX 392 IP We 5 | 0 | 
XX, 399 IP Vo ~ 45 1 | 
XX 731 IR vo 11 Il | 
XX 1425 1 Ao 10 | 0 
XY 1370 ? vo? | 5 0 | 
XZ 448 M vy 4 0 | 


| 

| 

The interactions listed here include all those observed in the 9.153 photographs referred | 

to in Table I. The origin P indicates the paraffin target, M indicates a doubtful origin that | 
could be in aluminium, copper or paraffin. «V°» means a neutral V-decay that cannot be 
separated as A° or 6°-decay. For event XY 1370 the origin is not coplanar with the V°-decay. 
This could be an example of anomalous, 3-body decay or could be due to the presence of a 


second, unseen, interaction. 
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Little comment is required. There is a tendency for the interactions to 
have higher multiplicities than the average but the difference is barely sig- 
nificant statistically. A high proportion of the interactions occur in the pa- 
raffin target; this could be due, in part, to the higher probability of observing 
a V°-decay when the V-particle is produced close to the sensitive volume of 
the cloud chamber. 

If we confine attention to the interactions with between 5 and 15 secondary 
particles we find the following results. There are 129 interactions with a total 
of 953 fast charged secondary particles. The observed number of V°decays 
is 8. Making approximate allowance for the particles decaying outside the 
cloud-chamber and for those decaying into neutral secondaries, we estimate 
that about 40 V°-particles must have been produced in the interactions. This 
corresponds to about 4% of the charged secondary particles — a similar pro- 
portion to that found in many experiments concerned with interactions of 
much lower energy (*). 


8. — The nature of the interactions. 


81. The significance of y.. — The physical meaning of y, depends upon 
the nature of the interacting particles and on the model that is adopted for the 
description of high-energy nuclear interactions (4). 

In a collision between two protons the secondary particles must, statisti- 
cally, be symmetrically emitted in the centre-of-mass frame of reference. 
If this symmetry, which is demanded by the principle of special relativity for 
an ensemble of interactions, is assumed to exist in every individual interaction 
then the C-system that we have considered is the centre-of-mass system of 
reference. The relation between y, and the energy (Z,) of the primary protons 
in the laboratory system is then 


— 
=] 
<= 


E, = (2y, —1)M, » 27°M,, 


where M, is the rest-energy of the proton. 

If a primary z-meson of high energy strikes a proton there is no a priori 
reason to assume forward-backward symmetry of the secondary products in 
any frame of reference. However, symmetry of the produced particles in the 
centre-of-mass frame of reference is a not unreasonable hypothesis and for an 
interaction with large multiplicity our C-system may again be identified with 


(?) €. FRANZINETTI and G. MorpuRGO: Suppl. Nuovo Cimento, 6, 469 (1957). 


(4) Article by D. H. PERKINS in Progress in Elementary Particle and Cosmic- Ray 
Physics, vol. 5 (Amsterdam). 
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the centre-of-mass system of reference. The laboratory energy of the 7-meson 
is again given by 2y:M, for values of y. such that y > 1. 

The third type of collision we consider is that between a single primary 
proton and a number of nucleons inside a target nucleus. If the interaction 
can be separated into a number of distinct nucleon-nucleon and r-nueleon 
reactions occurring successively in a cascade process, no single frame of re- 
ference can be used to describe the complete interaction. If, instead, the 
incident proton can be assumed to interact simultaneously with a group of 
n nucleons in the target, a centre-of-mass frame of reference exists. If » is 
small compared with the multiplicity of secondary particles and if the produced 
particles are symmetrically emitted in the centre-of-mass system there will 
again be an approximate correspondence between this centre-of-mass system 
and our O-system. 

If all these conditions are satisfied the relation between y, and the labo- 
ratory energy of the proton is E, & 2ny°M, if }>1 

The above considerations show how important it is to decide on the nature 
of the interacting particles. 


8°2. The nature of the primary particles. - Whatever the exact significance 
of y. there is no doubt that the primary particles responsible for our inter- 
actions have energies greater than 50 GeV—mostly in the range from 50 to 
200 GeV. In principle, the counter tray A (Fig. 1) ensures that the primary 
particles are charged and not neutral although discharges in A may some- 
times be due to an upward-moving secondary particle from an interaction 
produced by a neutral primary. 

The only other direct information we have is that the primary particles 
are isolated in the sense that they are not accompanied by associated charged 
particles that traverse our cloud chamber. This suggests strongly that if our 
« primaries » are « secondaries » from atmospheric nuclear interactions the latter 
generally occur more than 50m above our apparatus. 

Indirectly, we may use the results of earlier experiments (°) on the «local 
penetrating showers» which showed that about half the primaries of such 
showers were charged and half neutral. From this fact it was inferred that 
they were high-energy protons and neutrons which, at this depth in the atmos- 
phere, are expected in roughly equal numbers. 

Our conclusion is, then, that the majority of our interactions are produced 
by high-energy protons and that these protons have travelled more than 50m 
from the interactions in which they originate. A small proportion (~ 10%) 
of the interactions may be produced by z-mesons or other unstable particles. 


(5) W. D. WALKER, S. P. WALKER and K. GRIESEN: Phys. Rev., 80, 546 (1950). 
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83. The number of target nucleons involved. — For interactions that occur 
in paraffin with a geometrical cross-section about 40% will be interactions 
with free protons and 60% with carbon nuclei. Of the latter a considerable: 
fraction (perhaps one third) should be « peripheral » and involve only a single 
nucleon in the carbon nucleus—apart from the transfer of some small ex- 
citation energy to the residual nucleus. The remaining interactions will be 
complex, involving several target nucleons intimately. Several different theories 
of complex interactions have been put forward and they predict very different 
values for the number of nucleons that take part in the interaction (°). 

Ideally, the positive excess among the secondaries of each interaction gives 
directly the number of protons involved in that interaction. Since all our 
interactions show secondary particles whose charge cannot be determined we 
cannot use the positive excess in this way. Some idea of the average positive 
excess is given by the figures in Table IV. Despite the uncertainty due to the 
large number of unmeasured particles the figures suggest that the average 
number of target nucleons involved is not much greater than one. 


TapLe IV. — The positive excess in the analysable interactions. 
Number of | Mean + ve 
interactions ry n, TE Ny excess 
(5<np<15) | (deduced) 
129 48 366 261 326 1.6 


The interactions with ny between 5 and 15 are collected together in order to find the 
mean positive excess. ny is the number of slow, heavily ionizing, particles observed. n. and n_ 
refer to the fast positively and negatively charged particles and n to the particles whose sign 
could not be determined. 

The deduced positive excess is found by assuming first that 90% of the heavily ionizing 
particles are protons and, second, that the ratio of positive to negative particles among Ne is 
the same as the ratio ni: n_. 


It is clear that, assuming uniform detection efficiency, it is reasonable to 
suppose that more than half the interactions we observe in paraffin are nucleon- 
nucleon interactions. In the remainder, the number of target nucleons in- 
volved seems to be small. 


84. The laboratory energy of the primary particles. — As a first approximation 
we assume that our interactions are simple nucleon-nucleon interactions and 
use (7) to calculate the primary energies from the values of y-: The results 
range from 30 GeV to 2000 GeV with the majority of values between 70 GeV 
and 250 GeV. 


(9) For example, T. Gozani and K. Sirre: Nuovo Cimento, 11, 26 (1959). 
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In this range of energies the cosmic-ray protons have a steeply falling 
energy spectrum of the form N(F)dEcx E-‘dE with s~2.5. If the inter- 
action cross-section is independent of energy we can use (7) to find the expected 
spectrum of values of y. We find N(y.)dy.o y “dy. which gives a power 
law with exponent —4 if s — 2.5. 

The observed distribution of values of y, (Fig. 4) is certainly much more 
widely spread than the expected distribution. There are two effects that could 
account for the discrepaney. In the first place, the selection of our events 
discriminates against low-energy interactions but does not impose a sharp 
cut-off. This undoubtedly reduces the number of events with low values of y,. 

Secondly, the large error in the determination of y, coupled with the steep 
form of the proton energy spectrum leads to a systematic tendency to over- 
estimate the energy of interactions. The number of low-energy interactions 
is always greater than the number of high-energy interactions. As a result of 
fluctuations there will always be more low-energy interactions that appear to 
have a high energy than the reverse. 

Some idea of the magnitude of this effect can be obtained in the following 
way. First, we note that a set of interactions with fixed true value of y will 
give a distribution of values of y. that is symmetrical about the true value 
on a logarithmic scale. We therefore change to the quantities /=Iny and 
l.=Iny,. We make the assumption that the distribution of /, about 7 can 
be represented by a Gaussian with standard deviation o. The true frequency 
distribution of / is obtained from the expected spectrum of y and is N(/) dl = 
= exp [21(1—s)]dl. The likelihood that an interaction with true value (J, dl) 
should occur and lead to an observed value (2, dl.) is thus 


exp [2(1 — s)] dl exp 


The logarithm of the likelihood function for is 


se (l Er L.)? 


(8) log LZ = 211 —s) 


and this has a maximum corresponding to 


= exp [2(s — 1)o*] 
y 

or, for 8 = 2.5, 

(9) 5 — exp [3°]. 
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For the majority of our events where the statistical error in y, is of the order 
of 25% the corresponding systematic overestimation of y, is less than 20%. 
However (9) shows that the bias rises rapidly with increasing error and it 
would be quite reasonable to take the few values of y, that we find to be 
ereater than 20 as being spurious. 


8°5. The inelasticity of the interactions. — For a nucleon-nucleon interaction 
the kinetic energy available for particle production is, in the center-of-mass 
system 2(y,—1)M,. The energy radiated as m-mesons is n_E*, where n, 
includes both charged and neutral r-mesons and E* is the average total energy 
of a 7-meson in the center-of-mass system. If we neglect the nucleons and 
heavy mesons among the secondary particles from an interaction we may 
express the inelasticity (A) in the form 


1.5, HÉ 


(10 e E, 
fie MO DENT. 


Experimentally, H* does not seem to vary much from one interaction to 
another. The quantity n,/(y.—1) should therefore be closely related to K. 
The distribution of this 

quantity for our 36 inter- 

een actions is shown in Fig. 9. 

If a reasonable value of 

E* (~ 500 MeV) is inser- 

ted in (10) we find that a 

| | | mu Jo of unity for 1,/(y.—1) 
ZO SO 40 corresponds to A~0.4. 

Yo! The distribution of 

Fig. 9. — Distribution of the quantity n,/(y, — 1) for the Fig. 9 must not be given 
first 36 interactions of Table II. n, is the number of too great an importance 
fast secondary particles from each interaction. The because of the very large 


relation between n,/(y, — 1) and the inelasticity is di- 
scussed in the text. 


| 


oO 
re 


Number of events, 


(=) 


(=) 
= 
(=) 


errors in y,: In particular, 
the low values of ”,/(y.—1) 
correspond to the events 
that give very large values of y, and, as explained above, there is reason to 
think that these may be spurious. 

The conclusion that we draw is that the inelasticity in the interactions we 
have studied is generally significantly below unity (~ 0.25 --0.7 5) but 
can say nothing about the form of the distribution. 


we 
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9. — Discussion and conclusions. 


The first conclusion we would draw from our work is that a cosmic ray 
study of proton-proton interactions at energies of about 100 GeV is perfectly 
feasible. At an altitude of 3500 m, using apparatus with an acceptance of 
40 em’ < sr and a target thickness of 0.1 interaction length about one useful 
interaction with energy above 70 GeV was observed per day. 

Secondly, the type of study we have made is severely limited by the large 
uncertainty involved in the indirect estimation of the energy of the interactions. 
The statistical error in y. amounts to about 25% in our method of determi- 
nation; for low multiplicities there is a tendency to choose the axis direction 
in such a way as to overestimate y,; finally, due to the shape of the spectrum 
of cosmic ray protons there is a systematic tendency to overestimate y, for 
all interactions. 

With these errors, it is clear that no quantity that depends strongly on y, 
can be well determined. The inelasticity of the interactions is the most im- 
portant and interesting such quantity. The quantities that are not dependent 
on y, or vary very slowly with the assumed value, are the transverse momenta 
and the anisotropy in the angular distribution of the secondary particles. For 
the transverse momenta, our results agree with many others in finding pre- 
dominantly low values (7). For the angular distribution of the secondary par- 
ticles we find a definite anisotropy in contrast to the results of DULLER and 
WALKER (8) and of LOHRMANN (?) who found no sign of anisotropy in inter. 
actions of comparable energy. We can find no satisfactory explanation of this 
discrepancy. 

Concerning the various models that have been proposed to describe high- 
energy interactions there is little that we can say since we have made no 
attempt to obtain an unbiassed sample of interactions. In particular no varia- 
tion of multiplicity with energy is found; our criteria of selection were such 
that interactions with low multiplicity were rejected or discriminated against. 

We have, however, two indications that the inelasticity of the interactions 
is considerably less than unity. The first is the presence of positive particles 
(presumably protons) with high momenta in the C-system. The second is the 
distribution of »,/(y.—1) discussed in Section 8.5. These indications are evi- 
dence against a model of the type proposed by FERMI (1) in which the entire 

(7) O. MINAKAWA, Y. NisHImuRA, M. Tsuzuxi, H. YAMANOUCHI, H. Aizu, H. HASE- 
GAWA, Y. Isun, S- TOKUNAGA, Y. Fusimoro, 8. Hasegawa, J. NisHimura, K. NIv, 
K. Nisnikawa, K. IMAEDA and M. Kazuno: Suppl. Nuovo Cimento, 11, 125 (1959). 

(8) N. M. DuLLer and W. D. WALKER: Phys. Rev., 93, 215 (1954). 

(*) E. LOHRMANN: Nuovo Cimento, 5, 1074 (1957). 

(19) E. FERMI: Progr. Theor. Phys., 5, 570 (1950). 
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energy of the interacting particles is statistically distributed among the se- 


condaries. 


In the initial stages of the experiment we had the valued collaboration of 
Dr. H. Frrsurn amd Mr. 8. O. Larson in CERN and much help from 
Mr. A. H. CHAPMAN of Imperial College, London. 

At the Jungfraujoch Research Station Mr. HANS WIEDERKEHR gave us 
every assistance. We are indebted to him and to the Administration of the: 
Research Station for providing us with facilities for ous work. 

While considering the results of our experiment we had an extremely va 
luable discussion with the groups at the Ecole Polytechnique, Paris, and Im- 
perial College, London, who are performing an experiment of a similar type. 
In particular this discussion made clear to us the importance of the argument 
summarized in Section 8'4. 

Finally, we wish to record our debt to the late Professor C. J. BAKKER, 
the Director-General of CERN, and to Dr. L. KowARSKI who encouraged our 
work and supported it whole-heartedly. 


RIASSUNTO (*) 


Si riferisce un esperimento con raggi cosmici per studiare interazioni nucleari di 
100 GeV con una camera a nebbia in campo magnetico. Si descrive l’analisi delle 
interazioni e si discutono le ipotesi necessarie per interpretare le misure. Si danno i 
risultati di un primo campione di 41 interazioni. Essi dimostrano che: 1) gli impulsi 
trasversali delle particelle secondarie sono generalmente bassi (< 0.5 GeV/c) e 2) nel 
centro del sistema di riferimento degli impulsi le particelle emesse per la maggior parte 
hanno impulsi bassi e la loro distribuzione angolare è anisotropa. Ci sono anche prove 
che in una interazione nucleone-nucleo di questa intensità di energia il numero di 
nucleoni interessato è piccolo. La proporzione di particelle strane fra le particelle secon- 
darie è analoga a quella che si riscontra in interazioni di energia molto inferiore. 


(*) Traduzione a cura della Redazione. 


2396 


IL NUOVO CIMENTO Vou. XVII, IN. 2 16 Luglio 1960 


Internal Conversion Coefficients of #2 Transitions. 
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(ricevuto il 18 Marzo 1960) 


Summary. — A survey and analysis of all available data on internal 
conversion coefficients of pure £2 (2 + —0 +) transitions is presented. 
The results indicate the possibility of the internal conversion coefficients 
depending on the deformation of the nucleus. 


1. — Introduction. 


It is well established (1-2) that proper charge and current distributions of 
the nucleus in its finite dimensions should be taken into account in all cal- 
culations of internal conversion coefficients. Further theoretical work (25) has 
shown that internal conversion coefficients also depend upon the detailed 
nuclear dynamics. However, all these effects are considered to be negligible 
for #2 transitions. For this reason, sometimes theoretical conversion coeffi- 
cients of #2 transitions are taken as a standard in order to evaluate conversion 
coefficients. But a few reports (*7) indicate that some conversion coefficients 
of 52 transitions show (10--20)% deviation from their corresponding theoretical 
value. All attempts (+%) at understanding these deviations in terms of nuclear 


(1) L. A. Surv and M. LISTENGARTEN: Zurn. Eksp. Teor. Fiz., 22, 29 (1952); 
L. A. Sriv and I. M. BAND: Coefficients of Internal Conversion (Moscow, 1956). 


) 
1, 326 (1956). 
(3) E. L. Caurcx and J. WENESER: Phys. Rev., 104, 1382 (1956). 
(4) T. A. GREEN and M. E. Rose: Phys. Rev., 110, 105 (1958). 
(5) A. REINER: Nucl. Phys., 5, 544 (1958); S. G. Nizsson: UCRL-3803 (1957) 
) 


F. K. McGowan and P. H. StELSON: Phys. Rev., 107, 1674 (1957); M. C. JosHI, 
B. N. SugB4A RAo_and B. V. THosaR: Nuovo Cimento, 10, 775 (1958). 


(7) G. J. NryGx and A. H. WapstRA: Nucl. Phys., 9, 545 (1959). 
(8) L. §. KISsLINGER: Phys. Rev., 114, 292 (1959). 
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structure effects have failed, because these Æ2 transitions are enhanced by 
large factors. These considerations prompted a survey of all experimental data 
on pure #2 transitions and an analysis of it to be made. 


2. — Analysis of data and results. 


Table I shows most of the available data (*) on Æ-shell internal conversion 
coefficients (x) of pure 2+ > 04 #2 transitions. In the first column, the 
nucleus and the energy in (keV) of the transition are shown. In the second 
column, experimental values (xÉ ) are shown together with the reference. In 
the third column the mean of all observed values weighted according to the 
inverse of errors are shown. In doing this, a 12% error was assigned wherever 
errors were uncertain or were not available, since this appears to be the normal 
range of experimental error. The fourth column gives the theoretical values 
obtained by numerical log-log interpolation of the values tabulated by Ross (7°), 
and SLIV and BAND (!. In almost all cases, these tabulated values do not 
differ by more than a few per cent. The fifth column gives the ratio R(«*) 
of the experimental to theoretical values. The last column shows the defor- 
mation parameter obtained from #2 transition probabilities (1). In Fig. 1 
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Fig. 1. — Variation of R(o*) with mass number. 


(°) Data were obtained in many cases from: D. STROMINGER, J. M. HOLLANDER 
and G. T. SsABORG: Rev. Mod. Phys., 30, 585 (1958); B. S. DzeLEPov and L. K. PEKER: 
Schemes of Radioactive nuclei (Moscow, 1958). 

(!°) M. E. Rose: Internal Conversion Coefficients (Amsterdam, 1958). 

(1) K. ALDER, A. Bonr, T. Huus, B. Morrenson and A. WINTHER: Rev. Mod. 
Phys., 28, 432 (1956). 


2398 


INTERNAL CONVERSION COEFFICIENTS OF H2 TRANSITIONS 191 
TABLE I. — Conversion coefficients. 
Nucleus AS E ron 
x KA E Deformation 
(E keV) “Theor. Jet) parameter, 8 | 
Observed Ref Adopted | 
ue 5.2-10-3 (1) 
(564) (5.2-+0.4)-10-3 (?) (9*£2250-6) = 10-2))0:3  10-* | 1.020511 = 
(6.6+1.7)-10-3 (È) 
dle 3.4-10=3 (6) 7 vw 
(603) SMIL (5) (3.8+0.5) -10-3 | 4.15-10-3 | 0.92+0.12 | — 
4.3-10-3 (8) 
126 Xe (1.6 +0.2)-10-? (7) a 65 .2)-10 
| .6: 0.2): —2 | 65: —2 LE € = 
Eee ped eg |(1-65+02) | 1.65-10-2 | 1.00+0.12 
130X6 10 (0) È | 
TE = .10-3 | 
(528) oe (20) (6.5+0.8) -10 6.6 -10 0.98 +0.12 = 
132X0 TH; i ¥ 1 | 
EEE a 10.39): 103 103 8 I ne 
(673) 2.7-10 () (2.7+0.32)-10 HO. lO) 0.75-+0.09 
154Ba | (4.7-40:3)-10-2 (11) 
(605) | (5.7-+0.7)-10-8 (23) e e Lera 
| . . 2 me | . bi ae . 20), aan 
(4.6+0.4)-10- | (22)¢, (19) (5.1 +0.8) 1.00 +0.16 
(6.1 +0.7)-10-% (ESE 
14006 | 5 i 
.3-10-4 59 È i o —4 È . 4 ale = 
(1597) 6.3:10 (15) (6.3+0.76)-10 6.9 -10 0.91 +0.11 
152Sm 0.80 (15) | 
MAD .65 28 | 
(122) val (85) DÉCO 0.65 1.08 +0.16 0.28 | 
DIE BAPE LO ca (3-25+0.3)-107 0.032 1.02+0.1 0.2 
(344) | (3.3+0.1)-10-2 | (35) na Sa) Wiis. eso Gal 
“aca f È 7 | l 
(123) (0.54-+0.13) (10) 0.54+0.13 0.64 0.84 +0.21 0.27 
160Dy 1.65 0.2 (2) 
(87) 1.6 (18) 5 diro ? 
tr (19) 1.75 -+0.20 1555 1.14+0.12 0.35 
20.20 a) | 


(a) These were obtained from «, using (K/L+M= 6.7). 


(b) This is derived, assuming the theoretical value for the «* of the 344 keV transition in !*Gd. 
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Nucleus Xexp. orco Ro) Deformation 
(H keV) = È parameter, 6 
Observed Ref Adopted 
| 164]) 
a (2.70--0.5) (21) Dap SO 2.19 1.23+0.22 | 0.41 
164Kr 
an 1.9 +0.2 (21) 1.9 +0.20 1.60 1.18+0.12 | 0.33 
| Er 1.9 +0.3 (27) 
| (22 LAZIO So 1.88+0.18 1.65 1.1540.11| 00.38 
ae © 3 . , : È 
1.85-£0.13 (22) 
10Y 152202 (17) 
(84) 1.6 +0.15 (25) 
| 1.56+0.15 (28) 
Le ace (27) 1.6 +0.15 1.36 |1.18+0.12 | 0.30 
1.69 (28) 
1.70 £0.20 (48) 
1.65 £0.12 (22) 
176 Hf 1.3 (17) 
(89) n. (22) L31201] 1.12 IHISSZO 0.29 
180]f 
(93) 1.30 +0.4 (859) 1.30+0.4 1.04 1.25+0.38 | 0.20 
186)g 0.37 +0.08 (29) 
(137) 0.35+0.10 (20) A co) 
Ne @) 39 +0. 43 0.90+0.14 | 0.20 
0.36 (32) 
186W 
(123) 0.45 +0.10 (38) 0.45 +0.10 0.57 0.79+0.17 | 0.20 
18808 0.40 +0.05 (31) 
(155) 0.32 40.05 (3) 0.36 +0.05 0.33 1.09+0.16 | 0.18 
Vest n) ùl 
192()$ 
(208) | 216 (49) | 0.16 0.02 0.152 | 1.05-40.13 | 0.18 
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TABLE I (continued). 
Nucleus Xexp. al R(o*) Deformation 
E ke es = : por arameter, B 
| ae Observed Ref. Adopted i i 
192Pt 
0.054 (34) 0.054 0.054 1.00+0.12 | 0.2 
(317) 
194py 
(328) 0.056 +0.006 (42) 0.056 +0.006 0.051 1.09+0.12 | 0.2 
196 Pt, 0.042 +0.003 (35) 
(354) | 0.059 (3) 0.046 +0.006 0.042 1.09+0.15 | 0.2 
0.05 (Es) 
Hg | 0.037 (5) 
(411) | 0.03 (54) 
0.028 +0.007 (38) 
0.026 (55) 0.0281+0.0005 | 0.0302 | 0.93+0.02 | 0.2 
0.025 (29) 
0.029 (65) 
0.028 1 +0.0005 (&5) 
202 
He | 0.034 (4) | 0.034 0.027 | 1.26 = 
(439) 
202 57 
Fe ii ©) 0.0058+0.0007 | 0.0058 | 1.00+0.11 ae 
(961) | 0.006 (8) 
206Ph 
9 -10-3 (49) (QT SET) 10-24)" 8-8-10-2) 1:08-.0.12 | eee 
(803) 
210Pp 
| 4.8-10-3 (5°) (4.8+0.6)-10-3 | 4.4-10-3 | 1.09+0.13 aes 
(1181) 
ES ie (8) | 0.13 0.111 1.17-40.14 | 0.2 
(24)1 
222 
tao, (52) = = n QE 
(111) 
2 51 
noie) MN Ge 0.19 +£0.03 0.180 1.06-+0.16 | 0.28 
(187) | 0.22 +0.02 ©) | 
— ——— |—— —— — = ila A 
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TABLE I (continued). 
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R(«*) is plotted against the nucleon number A, in order to bring out the pos- 
sible dependence of R(«*) on nuclear structure, in a good approximation. In 
this figure, the values of the deformation parameter, f, are given within brack- 
ets against the nucleus concerned. 

It is seen from this figure and Table I that in the region of nearly spherical 
nuclei (in the neighbourhood of A — 140, 208), R(«“) is approximately unity 
with much scatter on both sides and there is a trend of x, being slightly lower 
than on. This is in particular supported by the well studied (7) 412 keV 
transition in 1*Hg, as also by the transition in Ce, But the more impres- 
sive trend is the systematically higher values of most of the af, in the highly 
deformed nuclei. Many of them have been determined by different methods 
and different groups, as shown in Table I. Thus the possibility of a constant 
error entering the experimental values seems to be reduced. It was for this 
purpose that in the analysis of data the weighted mean procedure was adopted. 
Moreover the deviation of x&, of the 122 keV transition in 18m is larger than 
in the 344 keV transition in 1?Gd. Similarly, the deviation in the transition 
of ***Dy is higher than that in the transition of the less deformed '*4Hr. The 
general situation shows that experimental errors are rather large to permit 
a quantitative evaluation of the dependence on deformation. It may be em- 
phasized that apart from the cases which had been specifically reported (57) 
to be deviating from theoretical values, many other experimental data are 
in good agreement with the general trends of the large positive deviations in 
the highly deformed nuclei and a small and sometimes negative deviation 
in the near-closed shell nuclei. 

Most of the (24 > 0+) transitions in the deformed actinide nuclei fall 
below the threshold for A-conversion, so that not much information can be 
obtained in this region of nuclei. The dependence of conversion coefficients 
on the nuclear structure is expected to be prominent in heavy nuclei. In ad- 
dition, due to the high energy of such transitions in light nuclei and the low 
conversion in these light nuclei, not many cases have been studied accurately. 
For this reason, the survey and analysis have been limited to the rare earth 
and neighbouring nuclei. 

Similar analysis of conversion ratios X/Z and (£,+ L:;)/L1 showed that 
no systematic deviations are clearly observable. 


3. — Discussion. 


There are a few cases in which the generalization that conversion coeffi- 
cients of transitions in deformed nuclei is higher, appears to fail. These are 
the 123 keV transitions in !54Gd and 18%°W. No particular comments seem to 
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be relevant, since both these transitions have been studied in a rather approx- 
imate way and only one measurement for each of them is available at present. 
It is interesting that both of them are on the borders of the region of large 
deformation. The conversion coefficients of transitions in ?°*Hg and 27 Xe, 
similarly, shew large deviations compared to others in that region. Again 
only one measurement is available in these cases. Thus, while more accurate 
determinations are indicated by this survey for almost all transitions, these 
singular cases deserve particular attention. 

Experimental errors in conversion coefficients depend to a large extent on 
the particular experimental procedure adopted. There are many direct and 
indirect procedures (#2) of evaluating conversion coefficients. But the most 
favourable method which can lead to very small errors is to be decided for 
the particular case. Thus, as pointed out above, while the magnetic spectro- 
metric measurement of conversion line intensities and the corresponding f-group 
intensities has led to accurate values (7) in many instances, for highly con- 
verted and low energy transitions this method is very unfavourable. McGowan 
and STELSON (5), who first reported about the systematic deviations in a few 
nuclei, pointed out that systematic errors in their method could be a reason 
for such deviations. As has been stressed earlier, the consideration of values 
obtained by different methods obviates this doubt to a large extent. 

From such considerations, it is seen that possibly the conversion coefficients 
are nearly normal with a tendency to being slightly lower in the closed shell 
nuclear region, and are mostly above the present theoretical values in the highly 
deformed region. 

Most of #2 transition probabilities in the deformed nuclei are enhanced by 
large factors due to collective motion of nucleons. This makes the dynamic 
nuclear structure contribution (4) very small. The same conclusions were ar- 
rived at from considerations of nuclear configuration mixing (8). Thus, nu- 
clear dynamics may not be contributing appreciably to such deviations. How- 
ever, effects of static structure due to nuclear deformation can present themselves 
through electron wave functions entering into calculations of conversion co- 
efficients (*). 

It should be quite interesting, therefore, to understand the cause of these 
apparent anomalies from theoretical and more accurate experimental inves- 
tigations. 


(12) D. ALBURGER: Handb. d. Phys., 42, 26 (1957). 

(*) Recent investigations of the influence of nuclear deformation on electron wave- 
functions by P. Varner et al.: Zurn. Eksp. Teor. Fiz., 38, 870 (1960) and A. S. MeLIGY: 
Nucl. Phys., 14, 248 (1959) adduce further support for such a proposal. 
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I am very thankful to Dr. B. V. THosAR for many helpful discussions. 
My thanks are also due to Drs. M. E. ROSE, E. L. CHURCH and J. M. HOLLANDER 
for helpful communications. . 


RIASSUNTO (*) 


Si presenta una rassegna analitica di tutti i dati disponibili sui coefficienti interni 
di conversione delle pure transizioni #2 (2 + 0 +). I risultati mostrano la possi- 
bilità che i coefficienti interni di conversione dipendano dalla deformazione del nucleo. 


(*) Traduzione a cura della Redazione. 
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Nuclear Magnetic Relaxation Times 7, and 7, in Some Liquids. 


G. BonERA, L. CHIODI, G. LANZI and A. RIGAMONTI 


Istituto di Fisica dell’Università - Pavia 


(ricevuto il 30 Marzo 1960) 


Summary. Nuclear relaxation times 7, and 7, of a set of pure liquids 
have been measured. The results we have obtained give in general for 
low viscosity liquids, according to theory, T, —T, within the experimental 
error. We then examine the possibility that in some liquids having more 
than one chemical type of protons, different longitudinal relaxation times 
are present. 


1. — Introduction. 


Tt is known that nuclear relaxation times have a fundamental importance 
in describing the nuclear magnetic resonance in condensed matter. 

The reciprocal of the longitudinal relaxation time 7, gives a valuation of 
the interaction between the nuclear spin system and the lattice; the reciprocal 
of the transversal relaxation time 7, gives a valuation of the spin-spin inter- 
action. 

The nuclear relaxation theory in liquids has been developed by various 
authors (1°) and foresees 7, = 7, for pure liquids of low viscosity in which 


(1) N. BLOEMBERGEN, E. M. PurcELL and R. V. Pounp: Phys. Rev., 73, 679 (1948). 
() R. K. Wanasness and F. BLocH: Phys. Rev., 89, 728 (1953). 

(*) A. ABRAGAM and KR. V. Pounp: Phys. Rev., 92, 953 (1953). 

(*) R. Kuso and K. Tomrra: Journ. Phys. Soc. Japan, 9, 888 (1954). 

(5) I. SOLOMON: Phys. Rev., 99, 559 (1955). 

(5) F. BLocH: Phys. Rev., 102, 104 (1956). 

(*) N. BLOEMBERGEN: Phys. Rev., 104, 1542 (1956). 

(8) E. M. PurcELL: Suppl. Nuovo Cimento, 3, 961 (1957). 

(9) K. Tomira: Progr. Theor. Phys., 19, 541 (1958). 
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the relaxation is due only to a pure dipole-dipole interaction between the 
protons. The main aim of this work is to verify this expectation of the theory (*). 
A comparison between the two relaxation times in liquids has been already made 
by some authors (1-15), However the results till now published are not in good 
agreement with one another. Therefore we think that the present work can 
bring a contribution to the question. 

In this work we have then examined some liquids that have protons in 
non equivalent positions in order to verify the possible presence of more than 
one longitudinal relaxation time. 


2. — Experimental apparatus. 


The measurement ot the relaxation times has been carried out with the 
methods we have described in previous works (17:14). 

The measurement of the longitudinal relaxation time 7, is based on the 
recording of a sequence of nuclear signals in phase with the rotating magnetic 
field H,, due to the complete reversal of nuclear magnetization. This reversal 
is obtained performing adiabatic fast passages through resonance. 

The transversal relaxation time 7, is measured observing the nuclear mag- 
netization decay when it precesses together with the rotating field in a plane 
normal to the constant magnetic field H,. The magnetization M is set parallel 
with H, stopping an adiabatic fast passage in resonance. 

The measurements of 7, and 7, are made with a magnetic field A, of 
1650 gauss. 

The precision we have in comparing the two relaxation times in a liquid 
is determined in practice by the precision obtained in the measurement of T,. 


(*) If a sealar interaction with another spin is present and if the correlation time 
relative to this interaction is great enough 7, can be greater than 7,. Furthermore 
in this case the value of 7, can depend in a very remarkable way on the magnitude 
of the rotating magnetic field. [J. M. Winter: Compt. Rend., 249, 1346 (1959); 
I. SoLomon: Compt. Rend., 249, 1631 (1959)]. 

(10) J. G. PowLEes and D. Cutter: Nature, 180, 1344 (1957). 

(11) S. MerBoom and D. GILL: Rev. Sci. Instr., 29, 688 (1958). 

(2) J. G. PowLes and D. Curter: Arch. des Sci., 12, fasc. spéc., 135 (1959). 

(13) G. Hocusrrasser: Arch. des Sci., 12, fase. spéc., 132 (1959). 

(4) G. Boxera, L. CHiopt, L. GruLorto and G. LANZI: Nuovo Cimento, 14, 119 (1959). 
(15) J. G. Powzes and D. CuTtLER: Nature, 184, 1123 (1959). 

(16) I. SoLomon: Journ. Phys. Rad., 20, 768 (1959). 

(17) G. CxraroTri, G. CRISTIANI, L. GIULOTTO and G. Lanzi: Nuovo Cimento, 12, 
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Indeed the measurement of the transversal relaxation time presents greater 
experimental difficulties mainly if this time has a magnitude of several seconds. 

One could expect that in liquids having protons in non-equivalent posi- 
tions more than one longitudinal relaxation time and more than one trans- 
versal relaxation time are present. We have performed a research of this type 
for longitudinal relaxation times only, because of the greater precision we can 
obtain in this measurement. 

The investigation of more than one longitudinal relaxation time in liquids 
has been performed observing the return to the equilibrium value of the nuclear 
magnetization M after it has been reversed by means of a fast adiabatic pas- 
sage. If many longitudinal relaxation times 7, are present, M reaches the 
equilibrium value M, according to the law: 


M(t) = M, > A,(1 — 2 exp [—1?/T,)) , 


t being the time elapsed after the reversal and A; the fraction of protons having 
longitudinal relaxation time 7,. 

From the behaviour of M as function of tit is possible in practice to obtain 

the values of relaxation times only if they are two or three different ones. To 

this object it is convenient to set 
Ho In #(1 — M(t)/M,) versus t. A suffi- 
Î ciently sure experimental datum is 
then given by the slope of the tangent 
at the curve for t= 0. This slope is 
equal to > A,/T;. 

We have measured the rate 
| M(t)/M, relative to different values 
| of t modulating the constant magnetic 

field H, throughout the resonance 

(b) value H* as is shown in Fig. 1-a. One 
| obtains in such a way two fast adia- 
batic passages elapsed in time of #. 

Fig. 1. — a) Magnetic field vs. time; The height of the first signal 
b) record of nuclear signal. (Fig. 1-b) is proportional to M, and 

that of the second one to — M(t). 

We have obtained a modulation like that of Fig. 1 by means of the electro- 
lytic bridge described in a previous work (a) 

The central electrode of the electrolytic bridge is stopped outside the re- 
sonance during a quite long time in order that the nuclear magnetization reach 
the equilibrium value. Then by means of a particular device the central electrode 
makes a rotation of 360° with a suitable angular speed. 


| 
| 
| 
| 
| 
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3. — Experimental results. 


We plot in Table I the values of the longitudinal and transversal relaxation 
times in some pure liquids we have measured in order to verify the equality 
between this two times foreseen by theory. In the same table some of the 
data obtained by other authors are plotted. 


TARE 
Literature values Literature values | 
T, at | T, at |temp. | 
Liquid 1650 G|1650 G | TS 
4 a È È 4 (20) | Tha | field Itemp.| T, field | temp. 
à | x ON | A | Si | on 
| (s) | (G) (°C) | () [Mateo 
pera es - im ee ve Sa i e a 
Water EN RES ZERI NON 1650 23° | 3.1 a 1650 | 20° 
| 3.2 b 2000 and 7000 23° | 2.7 el 0.5 20 
0 mG Zo Oran 5 OO0OM Zr e227 Gn Owns 95° 
| oe 2.65d | 6500 | 24° 


d 6500 | 24° 


Cyclopentane | 12.8 | 10.5 | 20° | 


Cyclohexane | 6.2 fies |) Paes” (if @) 
| | 5.5 oc | 250 25% 13.5 el O75.) 2252 
| 
lec 5000 25° | 
Cycloheptane | 4.4 4.4 | 21° | 
Methyl-alcohol | 7.5 6.5 | 22° 8 ie 1650 (=) 
| (2.5-8.5)g | 250 25° | 2.5 g| 0.75 | 25° 
Ethyl-alcohol 3 aE 225 o Wii 1650 (19) 
| Whe Wit 250 25° | 2.3 g | (007500258 
| | 2.2(**)h | 6800 and 1100 | 20° | 
| = | | | A 
Propyl-alcohol | 1.1 | 1.1 | DID 7 250 25° | 1.9 g| 0.75 | 25° 
Butyl-alcohol | 0.6 | 0.6 | 19° | 14 gl 250 25° | 1.5 g| 0.75 | 25° 


References for literature values in Table I: 
a—C*) andy G4) G==(F2)'s c= (35); d= (1°); e=(); LC): TC h = (2). 
(*) These results, obtained some years ago with the same experimental apparatus, are rela- 
tive to a temperature that had not been exactly measured, but which was probably quite greater 


than 20 °C. | 
(**) This measurement has been accomplished in presence of oxygen in solution. 


(18) L. GruLorto, G. LANZI and L. Tosca: Arch. des Science, 10, fase. spéc., 250 (1957). 
(19) J. H. Simpson and H. Y. Carr: Phys. Rev., 111, 5 (1958). 
(29) G. CararorTri, G. CrISTIANI and L. GiruLorTo: Nuevo Cimento, 1, 863 (1955). 
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We can remark in some cases remarkable differences between the results 
relative to the same liquids. Only in part these differences can be due to dif- 
ferences of temperature or to the presence of oxigen in solution. 

If we restrict our considerations only to the results we have obtained at 
present, we can observe that in general the equality between the two nuclear 
relaxation times is well verified for liquids whose relaxation times are less or 
equal to about 4s. For 7, and 7, greater than 48 we can remark that our 
measurements should give for 7, a value slightly less than 7,. However it 
is probable that such a result can be due to a systematic error. These differ- 
ences could at least be partially due to surface effects. A systematic study of 
surface effects in nuclear magnetic relaxation is now in progress. 

We can remark besides that for some liquids plotted in the table (methyl- 
alcohol, ethyl-alcohol, propyl-alcohol and butyl-alcohol) many chemical types 
of protons are present. However, the return of nuclear magnetization to equi- 
librium value we have examined for these liquids with the method previously 
described, looks still exponential within the experimental error (see Fig. 2). 


Ti 5 
0 0 Ime tins 20 


Methyl Alcohol at 22 °C 


=3 | 


Fig. 2. - In3[1— M(t)/M,] vs. t for Methyl- Alcohol. 


This result seems to prove that for the primary alcohols we have studied, the 
longitudinal relaxation times of the various groups of protons in non-equi- 
valent positions, are not very different from one another. However this con- 
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clusion does not seem in good agreement with the results recently obtained by 
POWLES and CUTLER (1°). 

On the contrary an analogous research we have performed on toluene has 
shown in a rather evident way the presence of two longitudinal relaxation 
times as one sees in Fig. 3. The dotted line is the tangent at the origin of the 


Time t ins 
0 10 20 30 40 


Toluene at 24°C 


Fig. 3. - In {1 — M(t)/ My] vs. t for Toluene. 


curve which represents In 3(1— M(t)/M,) as a function of t; the continuous line 
is the theoretical plot for two longitudinal relaxation times (Ton, = 15 8 and 
(T.)cn, =88. This line is quite well adherent to the experimental results. 
We can observe that these values of the two longitudinal relaxation times of 
toluene are in good agreement with the results obtained by NEDERBRAGT and 
REILLY (2!) with a high resolution device ((Ti)o4,=168; (Tan, =98 at 
25 °O). 
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RIASSUNTO 


Sono stati misurati i tempi di rilassamento nucleare 7, e 7, per una serie di liquidi 
puri. I risultati da noi ottenuti danno in generale per liquidi puri a bassa viscosità, 
in accordo con la teoria, 7,=7, nei limiti dell’errore sperimentale. È stata inoltre 
esaminata la possibilità che in alcuni liquidi aventi protoni in posizioni non equiva- 
lenti siano presenti più tempi di rilassamento longitudinale. 
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The Characteristic-Functional Method in Cascade Theory. 


J. W. GARDNER 


Atomic Power Division, English Electric Co., Ltd. - Whetstone, Leicester 


(ricevuto il 30 Marzo 1960) 


Summary. — In a previous paper the last-collision equation treatment of 
the general Markoff cascade was developed and applied to examples 
drawn from cosmic radiation and nuclear reactor theory. The present 
paper gives an alternative approach via the characteristic functional, 
the latter being introduced as a logical extension of the concept of a 
moment-generating function. Potentially the C.F. embodies complete 
information about all the probability relations of the stochastic process 
in question; in practice the extraction of information about any par- 
ticular distribution depends on a suitable choice of the arbitrary auxiliary 
function occurring in the definition of the C.F. The power and versatility 
of the method are illustrated by reference to both point-source and 
extended-source cascades; it is particularly suited to the treatment of 
cascades from line sources, examples of which would be, in descending 
order of energy: a high energy nuclear «jet»; the track of an ionising 
particle; the source representation in the Feinberg-Galanin method of 
nuclear reactor theory. 


1. — Introduction. 


In a previous paper (1) the present author has given a general mathematical 
solution for the probability distributions of a Markoff cascade, the latter being 
defined as any multiplying system (physical, biological, etc.) whose growth can 
be characterized by parameters defining a Markoff process. An example chosen 
for illustration was the longitudinal development of the nucleon cascade with 
general source terms and certain simplifying, albeit inessential assumptions. 


(1) J. W. GARDNER: Nuovo Cimento, 16, 977 (1960). This paper will be designated 
hereinafter as I. 
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The method adopted was that of the so-called «last collision equation », which 
is itself a particular form of the Chapman-Kolmogoroff equation, relating the 
state of the system at #-+-At to the state at t. (Here t is the parameter with 
respect to which the process develops; conventionally it is referred to as 
«time» in the mathematical treatment although in physical applications it 
may be depth, energy, etc.) 

General solutions were obtained for the distribution functions p, and 4, 
whose definitions are here recapitulated for convenience: p,(11,+--, Mn 3 7) A1 +. An 
is defined as the differential probability that the postulated source distribution 
give rise at depth x to x particles with energies in the ranges 7, 4,+d, with 
k=1,..., nin any order; 9,(77, ... Nn3 ©) dm -. dm, is defined as the differen- 
tial probability that the postulated source distribution give rise at depth x 
to n particles with energies specified as above, plus any number of particles 
with arbitrary energies, including if we wish the same energies specified above. 

In addition to p, and q, a further function of considerable practical interest 
was defined, to wit g,(n; x) expressing the probability that the postulated 
source distribution give rise at depth x to n particles with energies greater 
than 7 and any number of particles with energies less than 7. Although it was 
pointed out that the moments of y, are simply related to the q, no equation 
was established for g, itself. This was because the Chapman-Kolmogoroff 
equations in question apply only to Markoff processes and the function g,(7; 7) 
does not define a Markoff process since the information that there are n par- 
ticles above energy 7 at depth x is insufficient to predict the state of the 
system at æ+Ax. 

However, @,(7; 7) is still a stochastic variable and it will be shown in what 
follows that it can be extracted from another quantity, the « characteristic 
functional », which obeys the Chapman-Kolmogoroff equations (*) and which 
contains in portmanteau form as it were, all the probability relationships asso- 
ciated with the state of the population of particles at a particular value of the 
parameter ¢ (or x). Whether these probability relationships can be extracted 
in explicit analytical form from the characteristic functional will depend on 
the particular application envisaged; in general one must not expect to avoid 
considerable numerical work with this method, as with other methods of solving 
the cascade problem, in any but a few mathematically favourable cases. This 
necessity for extensive computational work appears to be an ineluctable con- 
comitant of the physical nature of the problem. What is claimed, rather, for 
the characteristic functional is that it provides a convenient means of handling 
in one package operation any or all of the probability distributions of the 


(*) In the present paper we shall use only the first-collision equation satisfied by 


the characteristic functional; a subsequent publication will extend the discussion to 
other equations. 
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cascade, some of which are extremely intractable to handle by other means. 
From this it also follows that hitherto unsuspected, or at least unformulated, 
relationships between the various distributions may well be disclosed by the 
characteristic functional method. 

Rather than proceed immediately to thé details of the method, it may be 
to some heuristic and epistemological advantage to review certain other 
methods of cascade theory which have been used with varying success in their 
own particular fields. This we do in the next section. 


2. — Previous methods. 


Each of the distribution functions q,, Pn and y, represents a historically 
different approach to the cascade problem: the product density q, is associated 
primarily with the work of BHABHA (?) and of RAMAKRISHNAN (3), and the 
Ph and y, with the work of JANossy (15), although subsequent developments 
have of course been made by many authors, see for example the bibliography 
given by GARDNER et al. (6). 

The Janossy function p,, could it be determined, would provide complete 
information about the longitudinal development of the cascade since it spe- 
cifies precisely the energies of all particles at depth x, arising from a given 
source distribution. However, just because it contains potentially so much 
information the function p, is in general enormously complex, and its diffu- 
Sion equations can be solved directly in only a very few ideal cases. This led 
Janossy to develop and apply to cascade theory the method of regeneration 
points, already a familiar tool in other contexts to statisticians, see for example 
the review by BARTLETT et al. (7). The method uses the following argument: 
at some point (the «regeneration » point) located in the interval d£ at depth 
x —& a primary particle undergoes a collision, giving rise to two secondaries (*) 


(?) H. J. BHABHA: Proc. Roy. Soc., A 202, 301 (1950). 

(5) A. RAMAKRISHNAN: Proc. Camb. Phil. Soc., 46, 595 (1950); 48, 451 (1952). 
(4) L. JAnossy: Proc. Phys. Soc., A 63, 241 (1950). 

(5) L. JANossy and H. MEessEL: Proc. Roy. Irish Acad., A 54, 245 (1951). 

(6) J. W. GARDNER, H. GELLMAN and H. MesseL: Nuovo Cimento, 2, 58 (1955). 
(7) 


(*) Note: (i) As in I, we use the word « collision » only to fix our ideas: the mul- 
tiplication could occur for example by spontaneous fission, or some other process not 
involving a collision in the conventional sense. 

(ii) One of the « secondaries » is the original primary. 
(iii) We neglect ternary and higher order collisions which could however be 
accommodated by an obvious extension of the notation. 
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which become independent primaries of stochastic processes; the first-gene- 
ration secondaries of these new processes become primaries of further pro- 
cesses, and so on. This argument leads to the following well known first-colli- 
sion equation for a point-source cascade: 


ei 1 
cio, fon [ (@ = E)]pi(nr; É)Drlims EVO (Mz Nm) DN Am BE , 


kyim=n . 
O 7 


where (1.5 Hm) dx dj, AE denotes the appropriately normalized (*) probability 
for a primary particle of unit energy to have in distance dé a collision re- 
sulting in two secondaries with energies in the intervals dy, at 4, and dy» 
at Nm. 

Equation (2.1) may be formally solved either directly by means of multiple 
Mellin-Laplace transforms (8), or indirectly using the generating functions 
introduced by JANossy (45). 

Despite the fact that y, contains only the one energy variable 7, contrasted 
with the n energy variables 7,,..., 7, of p,, its practical evaluation is still in 
general a formidable computational task even with modern automatic equip- 
ment. Means of approximating to it have consequently been considered with 
varying success by a number of authors, a favoured method being an expan- 
sion of y, in terms of its moments which may be, at least in principle, com- 
puted from the product densities q,. The limitations of this method have 
been discussed elsewhere by the present autor (*): its usefulness is for practical 
purposes restricted to rather small energy ranges and numbers of particles 
(n not more than 2 or 3, say). 

Summarizing then the results of earlier work on q,, p, and y, it may be 
said that while formal solutions have been obtained for all three distributions 
the numerical evaluation of these solutions is practically limited, with existing 
computing facilities, to a rather small range of cases. This range may be 
extended somewhat—usually at a certain cost to accuracy—by employing an 
expansion of g, in terms of its moments. 

Naturally the quantities q,, p, and g, are not independent and for com- 
pleteness we recapitulate below the well known relations between them, the 
most widely useful of which is (2.4) expressing the moments of y, in terms 


(*) The normalisation in (2.1) is such that the © defined in (4.10) below is unity. 
We shall later remove this restriction. 

(8) J. W. GARDNER and H. MessEL: Phys. Rev., 86, 808 (1952). 

(9) J. W. GARDNER: Nuovo Cimento, 5, 1368 (1957). 


2416 


THE CHARACTERISTIC-FUNCTIONAL METHOD IN CASCADE THEORY 209 


of the qn, 


œ 


| ET 
(2.2) PRE 2 fe fr, don nr ne Mami ©) Ad 


do th 

(2 3) 5 Pesi L 

2. CAVA) at DM) MA) 

with 

(2.3a) P(n, x) JE Seni) dano 
0 0 

(2.4) Lia) so. Ae) en, 
n n 


where 7,(7;) is the n-th factorial moment of g,(7;x) defined by: 


Ga) Tals) = a a). 


The above are general mathematical results which are not restricted to any 
particular cascade model or source distribution; however, in physical appli- 
cations it is standard practice to express energies as fractions of the highest 
primary energy, so that the upper integration limits would then be unity rather 
than infinity. 

The remainder of this paper will be devoted to consideration of the cha- 
racteristic functional as an alternative description of the cascade. 


3. — Characteristic functional - general definitions. 


The concept of the characteristic functional of a distribution is a logical 
extension of the ideas of the moment-generating function and the characteristic 
function, familiar from standard statistical theory (7°). We recall first of all 
the definition of the moment-generating function of some stochastic variable ¢, 
having a distribution F(¢). In terms of some auxiliary variable v the moment- 
generating function M(v) say, is defined: 


(3.1) M(v) = E {exp [ré]}, 


(19) See for example: M. G. KENDALL: Advanced Theory of Statistics (1955). 


14 - Il Nuovo Cimento. 
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where H{z} means the expectation value of ¢. The function M(v) always 
exists; moreover when » is imaginary, say iy, the moment-generating function 
is then a function e(u) of the real variable u. This characteristic function c(u) 
is related to the distribution function F(f) by the wellknown formula (7): 


L 


BD AF(¢) = lim eer 
where 
(3.3) Ay(î) = y + h) — y(È) 


for any function y, and where F(¢) is assumed continuous at & and ¢+h. The 
n-th derivative (if it exists) of c(u) at u=0 gives the quantity iH{l"}. 

It is possible to generalize c(u) by replacing the argument w by the pro- 
duct-sum: 


(3.4) gun. 


where the R; are values of any random variable extracted from the stochastic 
process C(t) and the J’; are corresponding values of some arbitrary function It). 
In practical applications /Y" would be chosen to facilitate the evaluation of 
whatever parameters of the distribution one is particularly interested in; more- 
over both ¢ and J’ may be—and in practice often are—continuous variables, 
in which case the product-sum is replaced by the integral [Rte] ar (2). 
The characteristic functional, denoted by €, is then defined as c(u) with y re- 
placed by the above integral; it may be considered as an alternative theoretical 
specification of the stochastic process. While having a less direct interpre- 
tation, it is potentially useful as a concise form of definition which moreover 
lends itself to mathematical manipulations, especially linear transformations 
or operations on the random function C(t). BARTLETT et al. in references (7) 
and (11) have discussed the application of characteristic functionals to a variety 
of stochastic processes, including such disparate examples as the turbulent 
flow of a fluid and the ageing of a biological population. 

We conclude the present section by remarking that the complete charac- 
teristic functional with respect to the distribution function F(¢) automatically 
contains all the factorial moments. For the case in which F(f) is a step 
function, say b for ¢ between w and w+h and 0 otherwise, the verification of 
this statement follows immediately from expanding the exponential in the 


(11) See for example: M. S. BARTLETT: Stochastic Processes (Cambridge, 1956). 
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definition 


(3.4) OCU t= FE [exp AGES) 5 
From this result and (2.4) of the preceding section it is clear that the cha- 
racteristic functional also contains all the product densities q,, although not 
necessarily in such a way that they can be explicitly extracted in closed form. 
We shall revert to this point in the next section. 


4. — Application to point-source cascades. 


Although the definitions given above permit one immediately to write 
down the formal expressions for the characteristic functionals of any arbitrary 
cascade, a better insight into their properties and general usefulness can be 
obtained by considering first the characteristic functional of a point-source 
cascade. Let us consider therefore the cascade initiated by a primary par- 
ticle of unit energy, at depth «=0, and having the multiplication cross- 
section w defined in Section 2 above. 

The quantity À with respect to which the characteristic functional is taken 
can in principle, as we saw in Section 3 above, be any random variable of the 
cascade. For simplicity and definiteness we shall in the present context take 
it to be the number N(7, x) of particles having energies >7 at depth x, arising 
fom the specified primary. Then if /"(7) is an arbitrary real function of 7 the 
characteristic functional is given by: 


il 


(4.1) | CLL yy, #) = E |exr fr) AN (7, all R 


0 


We may note at once from the definition (4.1) a relation between © and the 
Pn defined in Section 1, viz: 


1 Mn N 


(4.2) OT, 7), 7] = n DI Ni Ja [i DICO 1: Pris +) Nn} 0) An + di - 


As in the case of the relation between € and g, mentioned above, the value of 
(4.2) as a means of determining p, must depend on the distribution law; both 
relations simplify considerably for the case of a Poisson law. (See for example 
BARTLETT, reference (11), chapter 3). 

From C can also be obtained the probability generating function and mo- 
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ment generating function ot the cascade. Thus the probability generating 
function for the number of particles of energy >7 is: 


(4.3) Ge, n, 2) = > 2" y(n, «) = C[L,(u), ©], 
n=0 
where 
= In e when >, 
(4.4) Pays) 
0 when u< 7. 


The moments of this distribution are given by: 


= ' Où 

(4.5) Um N", x)= |- 5 M(s, 1, 2) a 
where 
(4.6) Misi, ay Genie: 
with 

—8 when «>, 
(4.7) 1T(u) = 

0 when u< 7. 


Yet another property of the cascade which can be obtained from C is the 
total energy distribution (7, 7), whose Laplace inverse L(«, 7) is given by 


DD 


(4.8) La, #) = fexp [— an] o(y, +) dn = CI:(m), ©], 
with 
(4.9) il (7) = — an. 


The above examples illustrate the usefulness and versatility of the cha- 
racteristic functional in yielding any distribution of interest connected with 
the cascade by an appropriate choice of the arbitrary function I’. 

We now proceed to show that the characteristic functional satisfies the 
first collision equation. Since the auxiliary function J’ does not enter into 
the ensuing argument, which is quite general, we shall for the remainder of 
this, section use the abbreviated notation C(y, x) for the functional of a cascade 
initiated by a primary of energy 7. Let us consider the change in O(#) in the 
initial small interval (0, Ax). In this interval there is a probability 
O(Mmr Hr) Am Amx Av for the primary to suffer a collision giving rise to two 
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secondaries having energies in the intervals dm, at 7, and dy, at %, with 
Mm + Mx =1. The probability for the primary to suffer no collision in Ax is 
(1— 2) Aw where, as in I, 


1 dl 
(4.10) 2 =| (nx, Mm) Ay, + foun, Te) UY are 
0 0 


L if 


== 2fom, Nin) Ai, = 2 Jolin, a) Ole 


0 0 


Since any secondaries produced in the interval Ax may be regarded as the 
primaries (with energies 7x, 7) of two independent cascades, we have: 


1 
(4.11) O(1, x) = (1— Q Ax) C1, æ — Ar) + 245(0%, x — Ax) ©(mi2) dn, » 


0 


where v has been written for 1—7,—y7,. On dividing throughout by Ax and 
proceeding to the limit in the usual way we obtain the first collision equation 
for_0, ve: 

2 1 
(4.12) (7. + 2) Chay = afew, œ) O(N Ho) An, « 


0 


The solution of (4.12) will of course depend on the functional form of the 
cross-section ©, which is assumed to be specified either empirically, or by a 
suitable mathematical model for the cascade under consideration. It will be 
convenient to postpone further discussion of the solution of (4.12) until we 
have set up the corresponding equation for the general-source cascade, and 
can then discuss the solution of both equations together. 


5. — Application to generai-source cascades. 


As in I we postulate a source function S(7, x) representing the probability 
of a primary particle’s being «emitted » with energy 7 at depth x. Thus, for 
the case considered in Section 4 of a primary having unit energy at zero depth, 
S(y, x) obviously has the form 0(j—1)¢(x), where the deltas are Dirac func- 
tions. For the cascade having more general source terms there are two alter- 
native methods of setting up the characteristic functional equations. The first 
and more obvious method is simply a generalization of that given in Section 4, 
with the auxiliary function redefined to secure consistency with the new boun- 


(mi 
nN 
= 
nN 
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dary conditions. The second method utilizes the fact that in the present one- 
dimensional treatment the most general source can consist only of point-sources 
and line sources, and that the latter can be considered as generated by recti- 
linear motion of the former. ‘We shall consider each of these approaches in 
turn and discuss their equivalence. ; 

Let us write Y[y(y), €] for the characteristic functional of the cascade 
arising from the source distribution S(7, 4); here y(7) is a new auxiliary function 
replacing the (7) of the previous sections: it must obviously satisfy the 
boundary conditions imposed by #(7, «) but is other wise arbitrary. The fami- 
liar arguments of the regeneration-point method may now be used to express 
the change in Ÿ in the interval Ax. It must be borne in mind of course, in 
applying these arguments to the present case that a change in Y can arise 


not only from the primary’s undergoing a collision in the interval Ax but 
Aa 

from the emission of | (7, x) dx new «primaries » from the source S. How- 
0 


ever, this latter possibility can be accomodated, formally at least, by a re- 
definition of the collision cross-section © to include particles emitted by S. 
We shall denote by / this redefined cross-section which may be regarded as 
an empirically, if not analytically known quantity. Thus f(7%:#7wm)d7:d7m dx 
denotes the probability for two «secondaries » with energies in the intervals 
dy, at 4, and dy, at 7, to be generated in the depth interval de at x, regard- 
less of whether these arise from collision of the primary or emission from the 
source S. By analogy with (4.10) we may also define a total cross-section B 
thus 


(5.1), B= | B(x, Mm) x + | B( ey Ym) dm = 2 fi B(1x3 Nm) Any, = 2 | B(x) Nm) Mm - 
0 0 - 0 0 


The upper limit of integration, formally taken as infinity above, would in 
practice be determined by the highest energy of particles emitted by S. 

The equation satisfied by Y is then, by regeneration-point arguments, 
found to be: 


UD li) 
(5.2) [x(); No) ©] -| exp [— Br] P[x(n); m, ©—T1B(m, Ne dp dr . 


0 0 


If we now differentiate this equation with respect to x, at the same time 
suppressing the auxiliary function and defining our energy scale so that the 
maximum energy 7, is taken as unity, equation (5.2) becomes 


; | 1 
(5.3) (È "E à) (1, €) = Bf, ©) (M3 M2) Am , 
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which compares with equation (4.12) for the point-source cascade, showing 
that the characteristic functional satisfies the same type of first-collision equa- 
tion in either case. 
| We now proceed to the alternative approach in which the extended source S 
is supposed to be generated by the motion of a point source, Q, say. This 
latter must have a certain probability, P, say, to emit in a small interval 
Ax a primary particle of given energy which will then originate a point source 
cascade, described by a characteristics functional © as in Section 4 above. 
In general P will be a function of the energy 7 of the emitted primary, and 
of the depth x; in fact it is immediately apparent that it is only necessary 
to identify P(7, 7) with S(7, x) to ensure that the moving point-source Q fully 
reproduces the cascade of the stationary extended-source S. (Things would 
not be so simple of course if we were concerned with the temporal as well as 
the spatial development of the cascade, but that is outside the scope of the 
present paper.) 
As before, we consider the variation of Y[#] in the interval Ax. There is 
a probability 1— W(x) Ax that Q emit no primary in the interval Ax, where 


(5.4) W (a) =| Stu, x) dn ; 

0 
and there is a probability S(7, ©) dn:Ax that it emit a primary in the energy 
interval dy at 7, in the detph interval Ax. Such a primary can then initiate 
an independent point-source cascade whose characteristic functional at depth 
x is C[n, x]. Hence 


(5.5) x) = [1 — W(a) Ax] Y [a — Av] + 


+ W(x) Ar ¥ [a — Ari] Con x — Ax) S(n, ©) dy . 
0 


The upper limit on the energy integral is replaced by unity if we use the same 
energy scale as in (5.3), and it then follows from (5.5) that W[a] satisfies the 
first-collision equation 

al 


Wie tu) ¢ [n, €] S(y, €) dy . 


2 + Wo) 


(5.6) # 


Equations (5.3) and (5.6) will be equivalent if the auxiliary function y(7) 
in W{y(n), æ] is appropriately chosen; in fact these equations may be regarded 
as definining a relation between y and © for the cascade in question. 

The solution of (4.12), (5.3) and (5.6) will be discussed in the next section, 
and some applications considered. 


cn 
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6. — Solutions and discussion. 


As has already been mentioned the methods to be used for solving integro- 
differential equations such as (4.12), (5.3) and (5.6) will be selected in the 
light of the analytical or empirical forms of the functions © and $, which are 
assumed to be known. In view ot this it would scarcely be profitable to 
proceed further with our (hitherto) completely general discussion, leaving © 
and S unspecified. We shall first of all therefore suppose « and S to be such 
that known functions K(y, x) and H(y, x) can be found to satisfy 


n 
(6.1) K (7, &) -{¢ [n — u, #] w(u, 7) du, 


0 
and 


No 
H(i, «) = [vtr —u, &]P(u, 4) du. 
0 


Thus, for example, in the ionization cascade discussed in I and elsewhere (12-14) 
if C[y, æ] is known for 7 < kI,, where k is an integer and J, is the ionization 
potential, X(n,) for given x is known as a function of n in the range 
kI, << (k+1)I,, and hence by iteration in the range (k+1)1, << (k+2)1,, 
and so on. 

Assuming then that K(7, 7) and H(y, 2) are available in analytic or tabular 
form the solutions of (4.12) and (5.3) for the specified initial conditions can 
be written respectively: , 


(6.3) O[y, æ]= exp[i/(M) — Q(M)x] + 2) [exp [— Q(y) |e — tr] | K(m, t) dr; 
(6.4) ln, «]= exp [iy(n) — Bin)xl + Bin) [exp [— B(n) [# —T]] H(n, 7) dr . 


0 


To express the solution of (5.6) it is convenient to define a function Y(#, «) 
as follows: 


N 


(6.5) Mia) = [Ot e] S(u, x) du. 
0 
(2) J. W. GARDNER: Phys. Rev., 94, 764 (1954). 
(1°) J. W. GARDNER: Proc. Phys. Soc., A 68, 1072 (1955). 


(14) J. E. Morar: Nuel. Phys., 1, 180 (1956). 
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Then, for the same initial conditions as applied above to (6.4), the solution 
Of (5.6) is 


(6.6) Pin, “|= exp weft Y(n, t) — 1] ar 7 


0 


The use of expression (6.6) to compute Ÿ requires Y, and hence 0, to be 
known. In other words, one must in this case solve the point source cascade 
first and then feed the solution into the general-source cascade equations. 
This would seem a tiresome procedure if one were interested only in the general- 
source cascade and one might try to avoid the computation of © by using the 
expression (6.4) for Ÿ. Naturally both (6.4) and (6.6) must give the same 
ultimate answer and since the latter depends on © so, somehow, must the 
former. As we saw in the previous section this dependence enters in (6.4) 
via the choice of the auxiliary function y. In principle y can always be ob- 
tained if we know €; in practice y can often be much more expeditiously ob- 
tained by applying simple trial and error procedure to a few likely-looking 
functions. If y cannot be obtained by trial and error then the computation 
of C cannot be avoided, and if one has C computed it is probabiy better to 
proceed via formula (6.6) for the computation of 4. 


7. — Conclusion. 


We conclude with a few remarks about other forms in which the character- 
istic-functional method might be applied to cascade theory. In the definition 


of Section 3 above the R; occurring in the product sum > R,/ are values 
i 


of any random variable extracted from the stochastic process C(t). For the 
applications considered in this paper we have taken À to be continuous and 
have in fact replaced it by the average-number distribution. It would how- 
ever have been perfectly legitimate mathematically to replace it by any other 
distribution, continuous or discrete, extracted from the random variable ¢(t); 
one could for example define the characteristic functional with respect to any 
of the distributions p,, 4,, 9, introduced in Section 1. It might well be asked 
of course, in view of the complexity of these functions themselves, whether 
a characteristic functional so constructed would be of any value for the prac- 
tical computation of quantities of physical interest. It has already been re- 
marked that p,, could it be determined, would give all the physical information 
about a cascade, but that this very fact means that p, must be an enormously 
complicated function. Would not this same limitation apply also to the cha- 
‘acteristic functional constructed from p,? Not necessarily: the virtue of the 
characteristic functional method is that one can often choose the auxiliary 
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function Jin such a way as to suppress parameters not of immediate interest, 
leaving a relatively simple, or at least computable expression in the others. 

Tt should also be remarked that in the present paper we have concentrated 
on the first-collision equations satisfied by the characteristic functional. These 
equations correspond to the Chapman-Kolmogoroff equations relating the state 
of the system at time At to its state at time zero. For the applications dis- 
cussed in Sections 4-6 above this seemed the most direct and useful approach; 
for other applications some other forms of the equation satisfied by the cha- 
racteristic functional might be preferable. 

It is hoped to publish in a subsequent paper a more detailed development 
of the alternative forms of the characteristic-functional method adumbrated 
in this section. 


The author is indebted to Dr. J. BUNBURY for valuable discussion, and 
to The English Electric Co., Ltd., for permission to publish. 


RIASSUNTO (’) 


In uno scritto precedente si sviluppava e si applicava ad esempi tratti dalla radia- 
zione cosmica e dalla teoria dei reattori nucleari, il trattamento dell’equazione della 
ultima collisione. Il presente scritto da, in alternativa al precedente, un altro metodo 
di affrontare il problema a mezzo del funzionale caratteristico, il quale ultimo viene 
introdotto come un’estensione logica del concetto di funzione generatrice di impulso. 
Potenzialmente il funzionale caratteristico contiene informazioni complete su tutte le 
relazioni di probabilità del processo stocastico in questione, in pratica il ricavare infor- 
mazioni su una distribuzione particolare dipende da una opportuna scelta della fun- 
zione ausiliaria arbitraria presente nella definizione del funzionale caratteristico. Le 
possibilità e la versatilità del metodo sono illustrate facendo riferimento a cascate sia 
con sorgente puntiforme sia con sorgente estesa; esso è particolarmente adatto al trat- 
tamento di cascate con sorgenti lineari, come per esempio, in ordine decrescente di 
energia: un getto nucleare ad alta energia; la traccia di una particella ionizzante; la 


rappresentazione di sorgente nel metodo Feinberg-Galanin per la teoria dei reattori 
nucleari. 


(*) Traduzione a cura della Redazione. 
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Production of Bremsstrahlung 
by Longitudinally Polarized Electron-Proton Collisions. 


SASABINDU SARKAR 


Department of Theoretical Physics, 
Indian Association for the Cultivation of Science - Jadavpur, Calcutta 


(ricevuto il 6 Aprile 1960) 


Summary. — The purpose of this paper is to calculate the differential 
cross section for the production of bremsstrahlung in the high energy 
collisions of longitudinally polarized electrons and protons considering 
the effect of the anomalous magnetic moment and the recoil of the proton. 
We have neglected the matrix element for the bremsstrahlung of the 
proton and so our formula will be valid when the angle between the 
directions of the electron and the photon will be of the order of mc/p 
(m and p being the mass and momentum of the electron respectively). 


1. — Introduction. 


In this paper we intend to calculate the cross-section of radiative collisions 
between a longitudinally polarized electron and a nucleon polarized in the 
same or opposite direction. Several authors, e.g. BINCER (1), FORD, MULLIN (?) 
and ScoFIELD (*) have investigated different cases of the differential cross- 
section of scattering of two polarized fermions. The effects of polarization in 
the case of bremsstrahlung production becomes important when the kinetic 
energy of the recoil proton is comparable to its rest energy. RZEWUSKI (*), 
DRELL (5), BERG and LINDNER (6 have studied high energy electron proton 
bremsstrahlung. In our case there will be two matrix elements one of which 


AS Me BINCER Pivjs. Revs, 107, 14679 (1957)- 
G. W. Forp and C. I. MuLLINn: Phys. Rev., 108, 477 (1957); 110, 1485 (E) (1958). 
J 


J. RzEwuski: Acta Phys. Pol., 9, 121 (1948). 
SD DRELL: Phys. Rev, 87.753 (1952). 


R. A. Bere and C. N. LINDNER: Phys. Rev., 112, 2072 (1958). 
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corresponds to the radiation emitted by the electron (M) and the other by 
the proton (M,,). There is a sharp rise in the value of the differential cross- 
section of ultra relativistie energy of the electron when the bremsstrahlung 
photon emerges nearly parallel to either the incident or the final electron. 
Tt can be shown that when the angle of emergence is of the order of me/p 
(m and p being electron mass and momentum respectively), and the electron 
energy is of the order of the rest energy of the proton (Mc?), the contribution to 
the differential cross-section from |M;,|? and that from the product of M, and 
M, is smaller than that from |M,,|? by factors of the order of (m/M)? and 
m|M respectively as observed by RZEWUSKI (1). We have considered only M, 
and so our formula for the differential cross-section is valid when the emerging 
photon makes an angle of the order of me/p or the kinetic energy acquired by 
the proton is small compared to its rest energy. The formula for the diffe- 
rential cross-section can be verified by coincidence experiments. The formula 
for the total cross-section obtained after integrating the differential cross- 
section over the momenta of the proton and the photon will be fairly correct 
because BERG and LINDNER (5) have estimated that in their case the product 
of M, and M,, contributes at most 3% to the final integrated cross-section 
when the maximum of the momentum transferred to the proton is of the order 
of 500 MeV. In our calculation we have taken into account the effect of the 
anomalous magnetic moment of the proton which is taken to be g in units 
of nuclear magneton. 


2. — Calculation. 


Let p,(8;), p(£), k(k) and q=p,—p—k(#,) momenta (the correspond- 
ing energies) of the incident electron, final electron, the photon and the recoil 
proton respectively. For convenience of calculation we introduce unit vectors 
e., e, and e, in the direction of ps, pox q and (piX g)X ps. 

Taking o, = 6° Po/|Po|, the spin projection operator for the initial state of 
the electron is $(1+-e,0,.) where e,=1 or —1 according as the electron is polarized 
io: ene forward or backward direction. The corresponding operator for the 
initial state of the proton at rest is taken to be 4(1+e,0,). 

The required differential cross-section d®(e,, ¢,) considering the matrix ele- 


ment M, of BERG and LINDNER (*) is evaluated in the usual way and is given 
Dy “(putting fi ¢ = 1). 


Ga Ki, M E, dE AQ dQ, 
(La I RITA 
era 
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where 
shi Ey KR) (He M} = 9? 4K H+ (#,— BE)? + — @ 
Pi = Der a 2A? > oe 9 “2 cal US, ok 
4 Put Pr, BE + À) + (A, M} — 
AA, 2A? ) 


Pi 2PaPi , Por A ee 
ne 2 » 2 | p 2, 
$2 = (A ee (5 110 a a È | + à 


2 D DENT: ae 2 ; 3 
ee Os és PoP A2 Ey k) 4 (Hy re ee (PS i Ps 


A2 12 IAE PAR Ao 
ky + 2Piea (E, — k) A — q?/2 kE — kh? Ay (dz — 26:19 di 
k E è e = f | y 0 0 x CUA x 
i gg ( 3 A, | DO A, A 
1 Fe +2enPa +2en(KX e. hr E IAT 1 
ois 24, FH) 4 
ikkyo + Porkkyde + Cyik? Aoda 
I Cyr A? ? 


SI 2PaPi | Pu , (kxexi | Pi i) + 
0 


an DAI I TR) ope ot OY A 
ke2—kA (pes Por ar Por Eur y + KA ni — à 3) Daga 
AE A A; T E fea | 
CxrPollikzA + hide) — Pordilk® — hz) | 1 2PorP1 — Bede — 2P re 92 — bdo 
| : Lp ee De rat 24 a 


Po1— Porea(2k, + q: ) + Pek, — Co Cyt hy Gn — kKAC 


k 7 3 


where the suffixes a, y, z and / denote the components of the quantities in the 
direction of e,, e,, e. and the polarization vector of the photon respectively. 
If we do not want to observe the polarization of the photon we have to sum 
the above expression over all its states of polarization. In the eq. (1) 

= (kH —k-p)/k, A, = (kE,— k-p,)/k, dQ and dQ, are the elements of solid 
angles about the directions of p and È. 

The expression (1) clearly shows that the effect of polarization increases 
with q/M. 

The above expression can be simplified by taking p,, p, k and the pola- 
rization vector of the photon to be observed to be in the same plane, then 
terms involving ky, @y:Pi, CyiPory Cyr€xt) CR enlkxe:); will all vanish. 
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When the angle between p and k is of the order of mp, the expression (1) 
can be put in the following simple form 


ro 


MOR DI E, — M 5 2) 9 
(2) CORRE era Il g aM Jus + k)+ (E, M) q )/ + 


E,—M €183 | E, dE 49490, 


a i ar (1+ 9) (GE +(1+9(E. Micol: i aa 


Further for E,= M, k — 1 MeV, the angle between p and k is of the order 
of m/p, pp, and 1+ 9 = 2.793 4 + 0.0003 


do, 
db, 


= 0.0315 using eq. (2) 


where suffixes p and a correspond to cases where the spins of the proton and 
electron in the initial state are parallel and antiparallel respectively. 

When q is small compared to M, kK<1 MeV and the angle between pho- 
ton and the electron (both initial and final) is not necessarily small, expres- 
sion (1) can also be put in the following simple form 


a, 2 E,— M 
6) ADE i teh BV (+e ai) (+ E224 


© 800 pg’ 4 AA, at 


i E> 
DI eee no 


| E,d4E AQ dQ, 
2M 


i kE,—k-q 


- 


(1 + 9) (GE. + (1 + g)(Es — M)q°/2) 


Finally the formulae (1), (2) and (3) should all be multiplied by the square 
of the form factor of the proton which is same for both its charge and mag- 
netic moment distribution. 

In the experiment of high energy electron-proton scattering some soft 
photons which escape detection are emitted by electrons which are counted 
in the experiment as having been elastically scattered, though, actually, they are 
not. For this case SCHIFF (7) has given a correction term which is obtained by 
integrating the differential cross-section for the production of bremsstrahlung 
over the momenta of soft photons and recoil protons. A similar correction to 
the experiment of high energy longitudinally polarized electron-proton scatte- 
ring can be obtained by integrating eqs. (2) and (3). 

We find that the radiative and non-radiative scattering in the forward 
direction is nearly the same for the two cases of the spins being parallel or anti- 


(7) L. I. ScHIrr: Phys. Rev., 87, 750 (1952). 
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parallel, in these cases the influence of the interaction is not much felt. 
However, the difference between them is most prominent when the influence 
of the interaction is most felt, i.e. the value of the cross-section of the scat- 
tering of the electron in the backward direction for the case when the spins 
are antiparallel is very much larger than when the spins are parallel. From 
this it appears that the interaction between an electron and a proton is much 
stronger when the spins are antiparallel than when they are parallel. 


The author is greatly indebted to Professor D. Basu for his constant help 
and guidance throughout the progress of the work. 


RIASSUNTO (*) 


Lo scopo di questo lavoro è di calcolare la sezione trasversale differenziale per 
la produzione di bremsstrahlung nelle collisioni di alta energia di elettroni e protoni 
polarizzati longitudinalmente, tenendo conto dell’effetto del momento magnetico ano- 
malo e del rinculo del protone. Abbiamo trascurato l’elemento di matrice per la brems- 
strahlung del protone e quindi la nostra formula sarà valida quando l’angolo fra le 
direzioni dell’elettrone e del protone è dell’ordine di mc/p (m e p sono rispettivamente 
la massa e l’impulso dell’elettrone). 


(*) Traduzione a cura della Redazione. 
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CERN, Theoretical Study Division - Geneva 


(ricevuto 18 Aprile 1960) 


Summary. — Combining the assumption of a Mandelstam representation 
with approximate unitarity one can get the pion-photon scattering ampli- 
tude in terms of the pion-pion scattering amplitude in even states and 
of the matrix element forp hotoproduction of pions on pions. In solving 
this problem attention must be paid to the fact that in the low energy 
limit this amplitude should become equal to the classical amplitude. 
With some simplifying assumptions we obtain an expression which contains 
four quantities: the two pion S-phase shifts, the total cross-section for 
photoproduction of pions on pions, and the electromagnetic coupling 
constant. Numerical estimates will be given later. The interest of studying 
this process les in its possible influence on the nucleon-photon scattering. 


1. — Introduction. 


In a preceding paper (1) the authors have studied the reaction VTT HT, 
using the Mandelstam representation and assuming a pion-pion resonance in 
the T=1, J=1 pion-pion state. From this reaction amplitude and the pion- 
pion scattering amplitude in even angular momentum states one can get the 
amplitude for pion-photon scattering, using again the Mandelstam repre- 
sentation, and neglecting, in the application of unitarity, all the processes other 
than two-body two-body. 


(1) M. GOURDIN and A. MARTIN: Nuovo Cimento 16, 78 (1960). See also: H. S. Won: 
Bull. Am. Phys. Soc., 4, 407 (1959). 
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The pion-photon scattering amplitude cannot be directly confronted with 
experiment, but it plays a role in the nucleon-photon scattering, when the 
photon is scattered by the meson cloud of the nucleon. More precisely, when 
one tries to apply the Mandelstam technique to the process YEN > y+, 
one has to know the imaginary part of the reaction amplitude in the channel 
y+y = V+; the lowest intermediate state is, apart from the x, giving rise 
to a pole, the two-pion state. This is the main motivation of the present work. 
On the other hand, the present problem constitutes a good illustration of the 
Mandelstam technique where an explicit solution of the equations can be 
given. 

In the body of the paper several simplifying assumptions have been made, 
but we show in the appendices how one can solve the problems in the general 
case. First, we only retain the J=1 part of the y+r >r+x amplitude, 
neglecting contributions from J >3; second, we neglect the D waves and 
higher waves in pion-pion scattering, because this is essentially a low energy 
treatment, and, due to the range of the forces, these waves are expected to 
be very small. Then it is possible to express the y-x scattering amplitude in 
terms of three quantities: 

— the electromagnetic coupling constant which enters as a multiplying 
factor in the low energy limit of the charged scattering amplitude (this limit 
has to coincide with the Thompson amplitude); i 

— the total cross-section o, for the process Y+T>T+7r; 

— the two S phases of pion-pion scattering in 7=0 and T= 2 isospin states. 

In the present paper we make no attempt to replace these quantities by 
numbers. We feel that for the moment it is not really worthwhile. One could 
in principle use the figures proposed by FRAZER and FuLco (?) in the 7=1 
J=1 pion-pion state, which gives, up to a factor, o,, and try to estimate the 
T=0.2, J=0 pion-pion amplitudes by inserting the T—1, J=1 data in the 
Chew-Mandelstam (*) or Cini and Fubini (4) equations for the pion-pion scat- 
tering. This will be done later. 


2. — Structure of the matrix elements. - 


We define p; p, as the incoming four momenta of the pions, k, k, as the 
coming momenta of the photons with polarizations e, and e,. The relevant 


(2) W. Frazer and J. FuLco: UCRL 8880, to be published. 
(3) G. F. CHew and S. MANDELSTAM: UCRL 8728, to be published. 
(4) M. Cini and S. FuBINI: to be published in Ann. of Phys. 


15 - Il Nuovo Cimento. 
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channels are indicated in Fig. 1. ‘We introduce the three invariants: 


2 
Pr Channel IL“ py [s=-(k+p) =— (he + po) ; 
a oe (1) Se 
Channel I | a inl les)? ns PRET 
Ki k2 with, on the mass shell: 
e; e) 
Fig. 1. (2) 81 + 82 +t = 2u?, 


where 4 is the pion mass. 
Due to energy momentum conservation all the quantities can be expressed 
in terms of three four-vectors: k,, k, and A= p, — Pi. 
The matrix element e,, ¢,, M,, (k, k,, 4) must fulfil the gauge invari- 
ance condition which can be expressed as well known by the following neces- 
sary requirements: 


KM e) and k 


If we retain only the terms giving non vanishing contributions we can write 
the matrix element ® as 


(3) P(e, €33 kx, ky, A) = D, (81, 82, t)L, + D, (81, 82, t) Ig, 


where the two invariants I, and I, are given by: 


(er ke) (€.° ky) 
ky ° ke 4 


Tu = (4° 62) 


(4) Ig = (e1°4)(e2: A) (Ky: h,) — (61° he) (€2° A) (Ke: A) Fi 


(61° ky) (€3° ky) (ky A)(k,* A) 


— (€,° A) (€:h,)(k* A) + (K+ Kp) 


With the above notations the S-matrix element has the form 


I 1 


2 


(5) SID = 


and we obtain for the differential cross-section y+n—>y-+7 in the centre of 
mass system: 


(6) 4 
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We notice that I, and J, are invariant in the exchange p, p,, i.e. 
A—— A. The matrix element ® must be even in A by crossing property. 
It follows that D, and D, are symmetrical functions of s, and sy. Then it is 
easy to see that in channel III only even angular momentum states of the 
pion-pion system are permissible; therefore ‘the isotopic T—1 state is for- 
bidden. Hence the two photons currents appearing in the matrix element 
are either both isoscalars or both isovectors. 

Let us call ®, and ®, the matrix elements for the Compton effects on charged 
and neutral pions (by charge conjugation one sees that there is no difference 
between x+ and x). The pions being in the T=0 and T—2 isotopic states, 
the corresponding reactions amplitudes are: 


1 
(7) toa 
D,=V8(D,— D,) = V2 F,. 
Conversely ©, =2(F,—F,) and ®, =4(2F,+ F,). 
Let us now compute the matrix element ©, in the Born approximation. 


The three contributions indicated on Fig. 2 cannot be treated separately be- 
cause of gauge invariance. 


We find for the two func- 
tions D, and D, the follo- ee SR SR 
wing contributions: Big. 2. 


2e? 


(Si — fl?) (S2 — 27) 


(8) Dy = 2¢, Di(S1, 82) = 


1 


where e is the rationalized electromagnetie coupling constant. 


3. — Kinematics in the various channels. 


Because of the crossing properties one has only to consider two independent 
channels (see Fig. 1). 


3°1. Channel I. — If k and k' are momenta of the pions in the centre of 
mass system, the three invariants are given by: 


| SEEN (k + wz) ; 
(9) ! 8, = pu? — 2kw, — 2k? così, 
| t=—2k*(1— così), 


where 0 is the scattering angle. 


2435 


228 M. GOURDIN and A. MARTIN 


Then 
2 D (Si — ue) 
i Pe ISD ©: 
(10) 28,8 
I a = = DICI 
cos 0 = 1 + PTE 


If we choose the gauge to eliminate longitudinal photons we get the fol- 
lowing expressions for I, and I, 


9 


a 


I, = ey re (e k,)(es-k), 


Tp ——= (8 = u)(er k')(ey- k). 


(11) | 
| 


+| bo 


In spite of the presence of t in the denominators the invariants I, and I, 
are bounded in the limit t >0. In this limit k= k'; so e, becomes ortho- 
gonal to k' and e, to k. In addition, it is easy to see that in the limit of zero 
energy, only I, is non vanishing while I, behaves like k. 


32. Channel III. — p and — p are the photon momenta, q and — q the 
momenta; then the invariants take the form: 


| 81 = — (¢ + wy + 240, COS @) , 
(12) | 8 = — (qd + 0? — 240, 0089), 
| t =4p = 40? . 


The pion-photon angle, y, is given by: 


13 rs: 
(13) DE es AY 


and I, and I, can be written for transverse photons: 


| L, si €;'e2) 


(14) 
I, = 2t(e,-q)(e,-q). 
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4. — Mandelstam representation. 


We assume for D (5,8) and D;(S:;S,t) a Mandelstam representation: 


Di 


(15) DSi, 8, 1) = Di + — lle Q:($1, 2) ds; ds, + 
(eo — 3) )(s2 — Sy) 


L 5:(81't’) i(S2, t ! c 
x 3; ; ss 
| Alle — $,)(t'— 1) dii Al Le t) Genet i=a,f. 


Let us first consider the case of two isoscalar photons. It can be easily 
seen that the variables of integrations s, and s, or s’ and ¢’ cannot reach 
simultaneously their lower limit. More precisely if (24) <s!<(4u) one has 
8, > (4u)?, and conversely; the same result holds for s’ and #. So the Cini 
and Fubini technique can be applied (4), and we introduce unidimensional 
representations. 

The situation is not the same for isovector photons. One can just say that 
the integration domains are limited by 


! 2 2 


ii 20 


2 


s, > (84), 


The weak dependence on s, and s, will be completely neglected. As will become 
clear later it is impossible to take this dependence into account as long as 
nothing is known on the reaction Y4T >T+r+4". 

The amplitude ®, involves both isoscalar and isovector photons but ©, 
corresponds only to isovector photons. 

Under these assumptions one can write generally 


AP 
(S1 — 4°)? 
+ Sym (s; — 82) + Di [t(t — ser Pas| 


(16) D (81, 82, =D" 2 (Si 3)? Pa È Fa | a” (ss) + 


Ch 


nia bit 
aa 0); 


where P, is the ordinary Legendre polynomial. 

This corresponds to a partial wave expansion in the various channels (see 
formulas (10) and (13). 

We assume for a™ and 5b” a spectral representation 


po (n) (/ / 
a” (8) i i Ma for the 7=0 amplitude, 
ug æ —S$ 
(2)? 
(16') af) = 0 for the T=2 amplitude, 
58%) = I TOY 
TT up a ae 


(24)? 
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The convergence of the integrals will be studied later. The low energy 
limits of D, and D, must, in addition, satisfy certain requirements (5) which 
will be taken into account later. 


5. — Unitarity in Channel I. 


In this Section we shall only be concerned with the PM, amplitude. In 
channel I, F, is given by: 


(17) F = (ex €2)Dzo 4 i (e: k')(es- hy) [Dro + 4K°(k + ©x) Dyo] - 


k°(1 — così) 
Application of the unitarity condition enables us to determine the imaginary 
part of this amplitude: 


(18) Im <ya|P | = > GalTin)(a]T*|y2) . 


(n) 


The only intermediate state we retain is the two-pion state, because we 
are interested into an energy domain where s, is less than (34)?. The 
y+r 77 amplitude has been calculated in a preceding paper (1). The 
structure of this matrix element is 


(19) My = 


1 è 2 
an ne De ae as 


where k is the photon momentum and q the outgoing pion momentum in the 
centre of mass system. If one assumes that the reaction mainly occurs in the 
T=1, J=1 state, H does not depend on the angles (this approximation was 
done in (1)). Under this assumption the imaginary part of the yn —> YT 
amplitude can be related to the total cross-section o, for yx +17: 


2(k + wx) 


20 Il = 
(20) m #, 


(e,x k):(esx k')or, 
with 


1° ka 
= — =. eee »2\ 12 
Or 487 (k ae Ox)? | H(k ) ae 


(9) M. GeLL-MANN and M. L. GOLDBERGER: Phys. Rev., 96, 1433 (1954); M. GeLt- 
MANN and V. GLASER: private communications. 
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where 20,=k+©,. By comparison of formulae (17) and (20) one gets: 
Im D, = 2k(k + x) cos6a,, 

21 

(21) SD 


Im Dgo = — Shee na Ra ‘ 


These formulae could have been obtained by use of the « optical » theorem. 

The more general case where H possesses an angular dependence is con- 
sidered in Appendix I. 

If we return to formula (16) for functions D, the knowledge of (21) gives 
the values of the weight functions o. Then one can write D, and D, as: 


Di(si82t) = Di (8182) + O,(8,t) + C,(S2t) + 
SIONE 


= Dt ) È ‘= LA n 


Ta: 
a. sep soa 
(22) | Gus) —[2st4 (s_u2)] = e 
(2)? 
sa 1. al O7(x) da 
g0(8) TT (x — u?)(x — 8) 


6. — Unitarity in Channel III. 


In this channel the unitarity condition connects the reaction amplitude 
Y+y >©+r to pion-pion scattering. Only even angular momenta are present, 
and we shall provisionally assume that the D-waves (and higher waves) of 
the pion-pion scattering can be neglected (the general case will be considered 
in Appendix II). The first invariant J, is isotropic but J, contains both S and 
D parts. Therefore, if Im D, is non vanishing in channel III the J#0 am- 
plitudes have a non vanishing imaginary part in contradiction with the above 
assumption. It follows that all the functions appearing in D; are zero and 
that in D, the only non vanishing weight function, 7° is proportional to the 
J=0 pion-pion scattering amplitude. We now extract the S part of F both 
in T=0 and T=2 states: 


di 


(23) (e. e») Fo.8 (4?) ere 


frites e»; q) dQ, . 
Dieneatote == 0 on 


(24) Im Fr(4°)= ty (40%) = FF (4?) hy (02). 
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hi (q?) = exp [ids] sin ds, is the pion-pion scattering amplitude for J=0 and 


T=0 or 2. 


7. — Low energy requirements. 


From the above equations we finally get the following representations: 


5 IL 67(x) da 

(25) Dyo(815 823 t) = DAS Sy) )= [2 Site (Si a) (a — u 2)(# — 8) Î 
Bre 
Lf osede 4 (Fi) 
+ sit + (a nl oe. =| 4(u? + gt) —t > 
(2)? 

te Rae he es eee 1 a p(x) dx il | On(x) dal 
28) Detti) = Dili [ata Late s) 


. 4 F“(0 ) hs" 2)d 2 
7) Das; 825 t) = Dz(813 82) # "i See : ? 


(28) D,,(51; S25 t) = D'(s1; S2) - 


It can be easily seen that in channel I (Compton scattering) lim 7 MESE Ds 
eed 
cos 0 fixed 


is a finite quantity. 

It is known that in the zero energy limit, for any charged system scattering 
photons (5), the cross-section becomes equal to the Thompson cross-section. 
Moreover the derivative with respect to k, at zero energy, is proportional to 
the magnetic moment of the system; here the magnetic moment is zero since 
we are dealing with spin zero particles; for the neutral case, the cross-section 
and its derivative should vanish. In fact the explicit proof was only given 
for spin 4 particles but everybody believes that this result is independent of 
the spin of the particle. We shall not give here a justification of this state- 
ment but we shall take it for granted. Then one gets the following condition: 


TD, + I,Dz— (I, Dz + 1,D5) ~ ke (8, — pt)? ~ (8. — pr) 1 


for s,— pw? +0, così fixed. 
Since I, behaves like (sj —?)?, the expressions (26) and (28) for D, and 
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D,, are unaffected by this condition, but in D,, and D_, a subtraction is needed: 


Re ; 1 Or(x) da 

29 Dio = Dz + [268 (8 ae Mier 

(29) 0 [ (s pu?) = @=u)@—= 8) "È 

(24)? 
more if. aS ) da Fo(q?) ho" (@?) dq? 7 

vz Jp (e—s2) a) (+ pt) [492 + ut) — 1)” 
RA poeta 

30 lp = 7 


These expressions still satisfy the unitarity requirements in channel ITI. 
The very fact that subtractions are feasible in equations (25), (27) natu- 
rally leads to ask the question: why other subtractions should not be made? 
Here, we have done the minimum number of necessary subtractions. It 
would look quite unnatural that the pion-pion interaction does not influence 
the k dependence of the y scattering cross-section, as viewed, for instance, 
in a perturbation approach. 


8. — Equations for Fî and F and their solutions. 


We take the average of the F, over angle g in channel IIT as was done 
in Section 6 (see formula (23)). It is more convenient to write equations for 
MESSA or wand fy E /0î: 


+1 
12 Wine 12 )d r2 e2 2 d cos 
81) fl LU Les al p 


—q— ie wi | © — gq? cos? p 


) 


+1 


g)dq" e d cos 
(32) fi) Re Liane al ET, 


a — ie wi i — q? COS? @ 
+1 œ a 
2 a(@ ) da 
side | Sym (1 =2)| - 
+2] cosa fi, Wooo NES) 
RI ~ (2u)? J 


These equations are obtained after some manipulations using the expressions: 
of g? and sin?g in terms of s,, s,, and t. These equations are of Omnès- 
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Mushkelishvili type (9); they could be solved directly but it is perhaps better 
to exhibit first the analytical structure of the inhomogeneous terms. By the 


change of variable Z =— 2/sin? y, one shows: 
aril — nè 
d cos y de 
(33) ZO 2 sa = = =" 
Wa — 9? COs? @ a, (q2 — 2) Vele + 2) 


One makes also use of the relation: 


— (@+u°)°/4x 
(34) pece fier = ele; 5 (see Appendix ITT). 
E que vw — So 2(q? a 2) Va(z + u® + IZ) 2) 
Then the equations become: 
s r fs dg ewe fk de 
(35) 2(q?) = TOI 2 | 
> g"— q— ie Wa (g— 2)Ve(2 + u?) 
figa ui d 
» 12 r2 e2u? 2 
(36) fi) et ort oe 
g—p—ie où 2)Ve(e 
co © — (23° +o)? 4a 
1 pees Hal: Pees de 
I e—u ‘an FE Sr sa 
IR n 


Actually the two last terms in the last equation could be combined to show 
that it is still meaningful if o,(7) increases like w'~® at infinity. 
Under the following assumptions: 


[ 350) = di(co) =0, 


One gets the solutions of equations (35), (36) (see Appendix ITI) 


F3(q?) = exp [o2(q?) + à dî(g2)] ey? 
1° 


| ot ade È [— 08()] | exp [— os(2)] — exp [— o3(— u?)] 


L Vale + yt) 


S I ! 2 ? 


Gi 2% z+ 4 


(9) R. OMNES: Nuovo Cimento, 8, 316 (1958). 
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where 


" 38(q ay dq'° 
S( 2 
03(q?) = 21 dE 


and 


Fo(q?) = exp [a0(q?) + i d6(99)]- 


— 3 
[ew de pero , exp [= 96(2)] exp [oil SE, 
I Vaetp)\ 2 24 we | 
— (a+ u?)?/4x 


_ 20% (eee fine. exp [— of (2)] de i] 
a x — p? 2) Ve(e + u?) i 


Qu)? — © 


(38) 


with a similar definition of 9. 

The advantage of putting the solution in this form, is that, apart o(q?) 
there is no principal value integral to evaluate, which is a great advantage 
from a computational point of view. It can be seen that F exists even if 
0,(r) behaves like «'~* at infinity: the assumptions on 6 show that o(z) behaves 
like (log 2/2) at infinity; hence for x + co 


(e+ u®)?/ È 


| exp [— go(@) | de 4) AO Le 
| J (@—-aviertrun ® 


The assumptions on the phase shift behaviour at zero energy and at infinity 
permit to guarantee the uniqueness of the solutions of equations (35), (36) 
(see Appendix ITT). 

By inserting F and F in the subtracted D,, and D,, (equations (29), (30)) 
and substituting these quantities in equation (17) one gets the Compton scat- 
tering amplitude. 


9. — Concluding remarks. 


We have given explicit expressions of Fj and F) in terms of o,, dî, dî 
which permit to obtain the y-x scattering amplitude or the m+m—yty¥ 
amplitude. However one cannot consider this problem as completely solved 
since we have not made any attempt of numerical evaluation of these ampli- 
tudes, because of the lack of information on pion-pion scattering. Nevertheless 
we felt it might be useful to give this algebraic solution. Of course it would 
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be very interesting to give numerical estimates of the deviations from Born 
approximation as the energy of the photon increases. This will be done in 
the future. 

Tt is important to notice that the equations we solve have a unique solution 
if the S-phases of pion-pion scattering are such that 6(0) = d(00)=0, but 
that the equations are not really unique, because an arbitrary number of sub- 
tractions can be made before solving them. Here, we have done one sub- 
traction, which is the minimum number of subtraction needed to ensure that 
the low energy cross-section for y-7 scattering is given by the Thompson for- 
mula. On the other hand one cannot see any mathematical argument to forbid 
a larger number of subtractions at t=0; nevertheless this would give the 
result that the derivative of the amplitude F,, related to the difference between 
the charged and neutral amplitude, with respect to k, at zero energy, would 
not be influenced by the pion-pion interaction, which looks rather unsatis- 
factory. So, on the basis of this argument and also for simplicity reasons we 
felt that one should not make more than one subtraction at zero energy. 


** OK 


The authors wish to thank Professor 8. FUBINI for useful suggestions and 
comments. They are indebted to Professor GELL-MANN and Professor GLASER 
for clarifying discussions about the low energy limits. One of us (M.G.) would 
like to thank Professor FrERZ and Professor BAKKER for hospitality at CERN 
during the early stage of the present work. 


APPENDIX I 


Unitarity in Channel I. 


The matrix element for the process y + x y + x has been defined in the 
centre of mass system as: 


(I.1) kh’, es | T|k,e) = ; ; P(e,, e,, k, k'), 


and for the process ytn—>x+x 


il 


È Il 
Je ‘q|T|k, e) = 
{ ) q | | 7) €, (277)8 4(kox)t 0, We, k, q) . 


be 
i 


cv 
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to 
w 
=] 


The unitarity condition has the simple form: 


ik Ò AU Le == k 
(2 zt)? 802 € q (2 Wg Ox jie 
yy (e,, k, q)¥ (es, k', q), 


(1.3) Im (e,, es; k, k') = 


and after integration over the modulus of q 


ig 
(1.4) Im Per, €>; k, k') = (82)? ” - fre, k, q) Le k', q) do, 


with 20,=k+a,. 
From gauge and Lorentz invariance conditions, we know that ‘Ye, k, q), 
has the form 


P(e, k, q) = (e, k, q)HK, (k-q)] 
It follows that the integral over the q angle can be written as a linear com- 


bination of two invariants directly related to I, and Z,. By definition of 
UW, ERO0Ì 186, 


(1.5) fre, kev gh (e,,. kg) dQ; = dg q@? {(e,-e,)(K-k')U,(k*, cos? 0) — 
Sa *k)(e.-k)U,(k?, cos? 0)}. 


By comparing this expression of Im ® and the definition 
Im D = I,ImD, — I, ImEDeE 


we easily deduce by identification: 


2 


ide GPE 


siete 0U,(k°, così 
247 k + x eae ik, ); 


Im D, = 
(1.6) 


Tr nere ni [cos (Us Su 
> 962 (k +) i 


At this stage, it is interesting to relate U, to the total cross-section for the 
reaction y+n—>7-+7:. In order to obtain the latter we must calculate the 


integral 
4 
| \We, k, g)|°dQ, = RU, 1), 


and the total cross-section is given by 


‘ 1 ka? 2 | 
= 1 ie 
KT) or = 487 (hb Lan) EE SUR 
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In the particular case considered in this paper, where H does not. depend of 
the angle between k and q, we simply have 


(1.8) U,(k?, cos? 0) = U,(k?, cos? 0) = | H(k*) |?. 


From the general form given for D, and D, we can expand their imaginary 
parts in Legendre polynomials of cos 6 


Im D; = > [2k(k + @;)}?" P.(c08 0) o (81), 


and the unitarity condition permits to know the weight functions of?(s,): 


cal 
i 2n +1 qk? L EURE 
oe (8) = gn = où [2k(k + a | RI Aa 
È \ —1 
(1.9) e 
È n+1 9 1 
op(S1) = 1927 (k + op [2k(k + 0)" fr —U,)—U;]P,(u)du . 
-1 


APPENDIX II 


For the channel III, is given by: 


(II.1) p = (€1° 2) Di — 8a4(e1: g)(€2°q)D, . 
Taking the incident momentum p as ¢ axis and the bissector of e, and es 
as x axis, we have: 


2 
(€1°q)(€2°q) = A sin? p[e,:e, + cos 20], 


where 9 and © are the polar angles of the final momentum q. Consequently 
the amplitude takes the form: 


p= (€1° €2){D,— È (qu) [1 — P,(cos p)]} D.} — (qa) P3(cos p) cos 2wD,, 


where P, and P; are respectively the Legendre polynomial and the Legendre 
associated function of order 2. 

The first part of g contains only cylindrical symmetrical waves and the 
second part non, cylindrical symmetrical waves. For one given value of the 
orbital angular momentum | we can define two partial waves (except for 
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at +1 
Faq) = —— 


+1 
fo Ego)? D, (1 — P,(cos g))]P..(cosg)dcosg, 
(IT.2) di 


+1 


for (q?) = [D.Pi(eos p) Psi(cos yp) d cos@ , 


—1 


corresponding to the following expansion for y: 


p = (e:°e2) > F31(9°) Pa(cos @) —i(q0,) > fu(q?) Pa (cos g) cos 2. 
1=0 


v=1 


Il 


The imaginary parts of F,,(q?) and f,,(q?) are easily related to the imagi- 
nary parts of D, and D,. With the definitions given in the text 


Im Dis = > (440)? Pa;(608 p) 01°2(40%) , 


I 


It follows that angular integrations over g can be performed and we define 
for the cylindrical part: 


+1 


Xn =fa- u?) Pon()Ponto(U) du ’ 


+1 


BS = fa ae U*)[Pon(u) du, 


-1 


and the imaginary part of F,,(q2) appears as a linear combination of the un- 
known weight functions 0: 


(11:3) Tm (97) = (4¢0-)?° 


AI le, dl me 
{eon — STI [ear 


2 


(409) + (440)? 8102 (404) + (Ago) ang dî) ; 


In the same way, we obtain for the non cylindrical functions: 


(11.4)  Imfp(4°) = (go) *[yua02 (409) + (490) 71 02 (409) + 


+ (Go) Yu 03 (4w?)] > 
where 


+1 


Vim -fa = u?) Pau) Pom(W) du. 


—1 
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It is easy to apply the unitary condition for each partial wave. For a 
given isospin, we have: 


Im P,:(q°) = Fa(0°) 21(°) ; 
Im f2:(q?) = fei (q?) 21(q°) o 


where 4,,(9?) is the pion-pion scattering amplitude for the angular momentum 21. 

By inspection of equations (II.3) and (II.4) we see that oŸ appears in 
Tine (gq?) and~o, in im nni ee 
assume the pion-pion scattering amplitude to be zero for the angular momenta 
higher than 2L, we can conclude: o) =0 for 1> L and 0)’ = 0 forl>L—1. 
The non vanishing weight functions are solutions of a linear system of (21 +1) 
equations (corresponding to (2L-+1) unitarity conditions) for (2Z-+1) functions 
(oY for 1=0, L and of for 1=0, L—1). In general we have a unique 
solution. | 

The particular case considered in the text corresponds to L —0 and we 
have only one weight function of given by the unitarity conditions corre- 
sponding to the S-wave. 


APPENDIX III 


We first transform the integral (34). A general technique has been given 
by the authors in ref. (1), but here a trick due to S. FUBINI gives the answer 
in a straightforward way: 


i fie 
© — 8, [y ag) 


Since s, is a linear function of cos y the last integral can be performed and 
gives a rational expression: 


foe) 


fa 
J 2400 + (+ 29)/49)"]” 


which, by a trivial change of variable gives: 


—(æ+ u°)/49 


dz 


(ITI.1) I 
2(q?— 2)Ve(e + ut) 


_ o 


2:40 
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Then we try to solve equations (35) and (36). Let us first solve the general 
equation 


| | MY fa") R*(g?) dq” — M< io] de 
(III.2) f(q?) = al pepe +f 29 oe. 
0 


— D 


where h(q?) = exp [10(q?)] sin 6(q?) and where p(z) may include a finite number 
of 6-functions and behaves at most like 1/2 at infinity. 
Equation (III.2) defines an analytic function f(Z) with the following prop- 
erties: 
— it is such that Jim Î(Z)= ©, if an oscillator behaviour of f(q?) h*(q?) at 
infinity is excluded; 


— it is holomorphic everywhere except along two cuts from — oo to 
— M and from 0 to +c; 


— the discontinuity across the right hand cut is such that 


f(q? + te)/f(q? — ie) = exp [27 6(q?)] 
If 


| 6(0) = 6(co) =0 and 


(IIT.3) 0(q*) ae q as q i 0 ? 


il 
Asa as q> 09, 


one can define 


1 Ô 12 d 12 : F 
Wi | DUR + ie) = MF. 


Then the quantity 


(IIT.4) exp [U(Z)] o : i SL Sa eee de] = exp [U(Z)]g9(Z), 


& — 


—_ œ 


fulfils all the above requirements. To check that this is a solution one can 
insert it in (IIT.2). Then one has to compute the integral 


1 exp [U(a + ie)] g(a) h*(x) d 1 pees [U(x +ie)] — exp [U(e — ie)]) g(w) de 
a i 7 Dis ESA 


16 - Il Nuovo Cimento. 
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This integral can be computed 
by following ©, (Fig. 3) (the con- 
tribution from the circle surroun- 
ding the origin vanishes). Then 


= Residue (Z,) X271 — 


Fig. 3. | il J 
1 exp [u(e + ie)]R*(æ)g(æ) dæ — I z 1 exp [u(Z)]g(Z)dZ 0 
dei; = exp ll TI à. 
0 d= 
Now, if g(Z) is regular at Z=— M (which is the case here) 
PI dll e . 
È exp[u(Z)g(Z)}AZ _ (exp [o(e)] [ge + te) — g(e — ie)] Tee 
(LIT.5) | SA TRY 
C_ — 
SU 
> 2ri] p(e) de i 
2z- Za 


— Oo 


Then it can be seen that the equation is satisfied. The multiplicity of the 
solutions has already been discussed elsewhere (7). The homogeneous equation 
has no solution when the hypothesis (III.3) are made on the phase shift. 
This is shown by inserting in the homogeneous equation the possible solutions, 
i.e. Z" exp[U(Z)], m integer (°), and using contour integration. 

The equations we have to solve in the present paper can be formally written, 
in the form (III.2), but the singularities of the weight functions make the 
integral divergent: one writes the term arising from Born approximation as 


_ nè 


È eu? de 
(ITT.6) —— | a 
D+) (g? — 2) Vale +p?) 


foo) 


— n° 


il dz dz 
= COTE | a | | = ! | = 
PHI (+ 2)Vee+ n) J (a 2e + pe) Vele + ur) 


If we disregard the divergences and apply formula (ITI.4) we finally get for 
the solutions convergent expressions. Though this procedure is not well justified 
we can see by contour integration that the result is correct. As we said the 
only trouble comes from the term (IIT.6) which appears both in the equation 
for fo and in the equation for fi. The contribution of (IIT.6) to the solution 


(7) M. GOURDIN and A. Martin: Nuovo Cimento, 8, 699 (1958). 
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is found to be 


fs(4°) = exp [u(q? + ie)]: 


= pe 


ee J exp [— o(2)] , exp [— o(e)]—exp[—o(— p)]] dz 
ET +e gd | {RA pò Vele + na) 
Now 
SAC C aC GP 
(i= Fate) — Sri 
DI » = Be 
| exp [u(Z)] dZ f[exp[—o(e)] exp [—o(e)—exp [—o( PE ni 
AIA Zn] Zoe | RU Vele+u?) * 
SUVA, sE] dZ I _ de De dZ 
È Hg (Z vai Zo)(Z ur Lu?) (Z SLA 2) Ve(e + u?) 2701 (Z PU Zo)(Z + pu?) 
os DI | — © Ole 
| exp [= ele] — exp [— u(2)] | exp [— el = exp [= et a) de 
Ze + pe Vale + i?) 


In the first part of the integral the path O_ can be obviously replaced by a 
contour surrounding Z,. In the second part of the integral the integrand has 
no discontinuity across the left hand cut and the only possible singularity 
inside C_ lies at Z=— u?, but the residue is zero, so, finally 


RI 


E (q'°)h gi eu? i dz 
: = f}(Zo Le .0.d. 
if Ge et f2(Z0) dari (Zi — 2) Vale + u?) 2 


RIFATSISIUINERTOR$) 


Accoppiando i presupposti di una rappresentazione di Mandelstam con una unita- 
rietà approssimata si può ottenere l’ampiezza di scattering pione-fotone in termini 
dell’ampiezza di scattering pione-pione in stati pari e degli elementi di matrice per 
la fotoproduzione di pioni su pioni. Nel risolvere questo problema si deve fare atten- 
zione al fatto che nel limite della bassa energia questa ampiezza potrebbe divenire 
uguale all’ampiezza classica. Con alcune posizioni semplificatrici otteniamo un’espres- 
sione che contiene quattro quantità: gli spostamenti della fase S di due pioni, la sezione 
trasversale totale per la fotoproduzione di pioni su pioni, e la costante elettromagnetica 
di accoppiamento. Stime numeriche saranno date in un secondo tempo. L'interesse 
dello studio di questo problema sta nella sua possibile influenza sullo scattering 


nucleone-fotone. 


(*) Traduzione a cura della Redazione. 
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Formulas for the Matrix Elements of Connected Diagrams 
in Field Theories. 


F. R. HALPERN (*) 


Scuola di Perfezionamento in Fisica Teorica e Nucleare dell Universita - Napoli 


(ricevuto il 9 Aprile 1960) 


Summary. — A concise formula for an element of the S-matrix has been 
given by CAIANIELLO in terms of the real particles appearing in the 
process. This formula contains in addition to the contributions of the 
diagrams that describe mutual interactions of all the interacting par- 
ticles with each other, the contributions of disjoint processes occurring 
simultaneously. In the present note a method of factoring out unwanted 
disconnected pieces is given. As an example the removal of vacuum- 
vacuum transitions in quantum electrodynamics is carried out. 


1. — Introduction. 


A compact formula for S-matrix elements in quantum-electrodynamics has 
been given by CAIANIELLO (1). This expression contains the contributions of 
all diagrams with a prescribed number of external lines. Among these are 
found not only those that describe interactions of all the particles but also 
those that describe the simultaneous interaction of two or more disjoint sets 
of particles, i.e. disjoint graphs. The most common example of this situation 
is the occurrence of vacuum-vacuum transitions in the conventional expres- 
sions for the S-matrix. In low order calculations where only a small number 
of graphs are encountered it is a simple matter to discard these irrelevant 
diagrams. However if a serious attempt is to be made to present the renor- 


(*) Fulbright Fellow on leave from Princeton University, Princeton, N. J. 
The research reported in this document has been sponsored in part by the Office 


Chief of Research and Development, U.S. Department of Army, through its European 
Office under contract no. DA-91-591-EUC-1096. 


(*) E. R. CAIANIELLO: Nuovo Cimento, 11, 492 (1954). 
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malized perturbation series it may be desireable to start by discovering expres- 
sions for an S-matrix element that is free of vacuum fluctuations. 

It has been shown (!) that several algebric algorithms are sufficient to 
express the general S-matrix element. In the present development a new 
expansion is introduced for determinants and with the help of this it will be 
possible to start at any point in a matrix element, internal or external and to 
determine which other points are connected to it. A formula will be given 
for the S-matrix element that is free of vacuum-vacuum transitions. There 
is no unique formulae for this since it is possible to begin at any point or 
group of points and find all the connected points, and for each such choice a 
new formula results. 

A typical matrix element M,, will be the product of the vacuum-vacuum 
transition V and M,,, the vacuum free matrix element. If each of these 
three quantities is expanded in a power series in e, the n-th order matrix ele- 
ment M}? is given by 


re [MI | N 
n Wry —28) 17 (28) 
(1) I Nago == > Myr J . 


s=0 \2s 


Equation (1) is a typical convolution type sum and may be solved for M® 
in terms of the other quantities quite easily. This solution although apparently 
more compat than the one to be presented is a typical combinatorial one in 
which a particular piece is added in and subtracted out many times to achieve 
a small net result, and thus will not be particularly useful for computational 
purposes. However if the matrix element M,, can be displayed in the form 
of (1) the identification of M® will be possible. The occurrence of an 
equation like (4) is one of the reasons that makes the separation of vacuum 
pieces much simpler than the separation of pieces connected by several lines 
self-energy and vertex parts. 


2. — An expansion theorem for determinants. 


The formula for the S-matrix in quantum-electrodynamics is (*) 


© 


D o (__qp\n 
) VES An > CS far faz anne 


5 Ai 
Vol... Toln=o N° 
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The following expansion theorem for hafnians is contained in (1), 


n 
0 0 
(3) [Se del DD (0 ie RASE 


Sat OnE 
The first sum is over the ("_') ways of separating «,,,... 2, into two groups. 
One containing n — s elements and the other s —r. This expansion requires 
that the points 7, x, are connected to the points 7,...7,. This follows 
since in the expansion of the hafnian 


(4) [e; oj] = [a æ;] = (x; è;] = D,(&; — Zi), LT: ;] — 0: 


A somewhat similar expansion is required for determinants. However the 
nature of the connection generated by the electron-positron line is somewhat 
different than that for the photon line since it proceeds either straight through 
the graph or is a closed loop. The set of all points connected to a given point 
by an electron line is the set of points on the line or loop containing the given 
point. The considerations lead to the following classification of the terms in 
the expansion of a determinant. Consider the determinant 


As, G2, An 


My, A, An 


In the expansion of this determinant the n! products are to be arranged in 
cycles. All of the cycles are to be divided up into two groups, those that 
contain the elements a,,..., a,, and those that don’t. The products are to 
be arranged according to the number of elements a’s there are in the cycles 
of the second type. There is a minimum of none, for example all the elements 
are in one cycle, and a maximum of n —r. A particular set of cycles of the 
first kind does not occur only once but will occur in fact in (n — s)! pro- 
ducts if there are n— s a’s not included in the cycles of the first type. The 
cycles of the second type are exactly those occurring in the expansion of the 
determinant composed of the » — s elements in the cycles of the second type. 
In a similar way a particular set of elements does not give rise to just one set 
of cycles of the first type but to several. In this case the determinant of this 
group of elements is not formed since at least one of the first r a°$ must occur 
in each cycle of the first type. However all the terms in the expansion of the 
determinant of this group of elements consistent with this restriction does 


occur. For this reason a determinant with a modified expansion rule is intro- 
duced as follows. 
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The determinant 


Oy A Ur o) Arti , ds 


Operon elite Os 


is to be expanded according to the conventional rules for determinants and 
then any term in the expansion that has at least one cycle composed entirely 
of the variables a,.,...ds is to be discarded. 

Thus the above discussion leads to the expansion of a conventional deter- 
minant as a sum of products consisting of one conventional determinant (cycles 
of the second kind) and one with the new rule for development (cycles of the 
first kind); more formally we have, 


Chose Dia n y see Ur. Gir, 4 U; 


Ri is “e 


I 
im 
M 


CAGE MO TETE TA ID aa 


EAN 


The inner sum is over the (7-7) ways of dividing the a,,...a, into 
two groups. A more formal proof may be given by counting terms on both 
sides. There are r(s—1)! in the barred determinant, (n —s)! in the un- 
barred and 


No signs occur because the structure of the cycles is left undisturbed in this 
expansion. The virtue of this expansion is that the points under the bar must 
be connected with the unbarred parts. 


3. — Separation of the vacuum-vacuum amplitude in electrodynamics. 


A systematic application of the two lemmas (3) and (5) to the S-matrix (2) 
permits it to be cast into the form (1). The first step is the expansion of the 
determinant. Only the determinant and hafnian are exhibited in what fol- 
lows since no alteration is made in the other quantities. Advantage will be 
taken of the #°s being dummy variables and the symmetry of the expressions to 
always replace the inner sums in (3) and (5) by a binomial coefficient. Ex- 
pansion of the determinant using the external variable as the unbarred lines 
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TS VIS Wii cnn HbR 
(6) ARSA A 
DN ao OO HP cao Hee 


Ra UD UTO LV | Baty one Va ae 
= Sla) [eb Co Om a iene 
0 


OP? eo NUS 


Because of Furry’s theorem the second determinant must be of even order. 
Tf it is of odd order there is at least one closed cycle with an odd number of 
elements which vanishes. Thus n — a = 2b and b may be used as a variable 
of summation. To simplify the notation the variables 7, ...7, are replaced 
by the subscripts 1.... With these modifications (6) becomes: 


ea | n UD ...U%? À... n — 2b, 
m (n) Gi 
(7) “LN, Py — > 
L= ORO D AU ED 2 D 
n— 2b, + 1...n 
feti re Len 
n— 2b, +1... 


The points 1...n —2ò, are now connected to the external lines. To sepa- 
rate the vacuum parts it is next necessary to expand the hafnian according 
to (3), using the external photon lines 2, 47° and those internal points 
Wy e Cn, that are connected to the external electron lines as the first group 
of terms. With this expansion (7) becomes 


na iz 6, | N DD NU ue Te ne 202 
(8) NA Pi > 


%=0%=0 \2h,/\2b,/\vD oo 1... n — 2b,/ \n — 2b, 4-1... n 
- [g® Pea Al n—2b,, n 2 + 1 n 2 JD 2b 1 
ean: Lou 19 n—- 2b, 3 cn N — 2b,][n — 2b, +1, «.. n]. 


The lower hmit for the b, sum is more precisely maximum (0, 2b, — (P, +n) /2). 
However the second limit imposed by the vanishing of the hafnian if there 
are two many circled terms will be ignored. 

In (8) the term b, = b, contains an obvious vacuum part. The remaining 
terms may or may not contain vacuum parts but all the points 1... n —- 2b, 
are connected to the external lines. To discover the remaining vacuum parts 
it is necessary to repeat the two expansions until all terms are found to be 
connected to the external lines. In practice this may be done by dropping 
the term b, = b, from the sum: note that when this is done the sum on bi 
begins at 1 not at 0, the remainder is then expanded again. It is not difficult 
to guess and then prove by induction the following result for the repeated 
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bo 
a 
(e) 


application of this process. 


(2/2) bi-1 dy gl Danes 


Me oe a ee Yay 


b,=c ba=c-1 bg=c-1 bi=c-2 b,,_1=1 


einen (26/26) Qboo-1\ (UV... uS® 1...n— 20, 
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n_2boaij+1...Nn_-2ba n-2ba+1...h_-2Dga4 


n_2ba3+1..h_—-2bg hn_-2boagt+1...n_-2Doa 
0 0 
[D 2%) 1 n—2b, n—2b,+1 n—2b,]- 


.[n—2h3 +1 n— 2bog4] VOM. 


Tf this expansion is continued until ¢ = [n/2]+1 the first term will vanish 
and only the second will remain in which all of the vacuum pieces are expli- 
citly exhibited. It may appear that the terms for which b,,,, = b,, have con- 
cealed vacuum parts but because of the hafnians these terms vanish; use will 
be made of this to achieve a slight simplification. To bring the matrix ele- 
ment into the form (1) it is necessary to move the sum on b,,,, to the extreme 
left as b,,,, is the index of the vacuum part. If we call the product of binomial 
coefficients, determinants and hafnians under the summation (BDH) after this, 


transposition we obtain, 
[n/2] [n/21—bsa41 [n/21 b,-1 bsqg_1—1 


(10) Lee sree > ae > PDH) 


Doat1=0 d=0 di=bog41+4 ba=b2941+4--1 boa=b2a41 


The binomial coefficients must be adjusted to bring out a factor (2b,, ,). 
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The required identity is that 


n 2b, 2004 n (p= DD Wea a, 20501 
= Sa = no 
2b,/ \2b, 2bsat1 Dear] \201 — 2beat1, 205q — 2bsata 


After making this substitution for B it is clear that the term b.4,,—0 
al VA ( D 0 RO 
is M© and is given by, 


AA (— 1)r7 
(11) NM ee a D PROS 
V5;! SO 
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Sey | DH. 


0 2b, 2b, 2054 


Where the terms with b,,:, — b,, have been eliminated and the exapnsion of 
this expression gives all diagrams without vacuum-vacuum parts. 


4. — Alternative expansions and applications to other processes. 


The previous expression was derived by introducing the external electron 
and-photon lines in a particular way. There is no particular reason for this 
method, it is possible for example to start at a single photon line and find all 
points connected to it. Since the expressions encountered in expanding the 
determinants and hafnians will be less symmetric if the approach is followed 
to some extent the inner sums in (3) and (5) can no longer be lumped into 
binomial coefficients this will make this work more complicated although this 
program may be followed to find connected diagrams. 

The extraction of self energies by this process does not seem very hopeful. 
It is possible to find self energies attached to a particular pair of points but 
this is a very cumbersome method of removing all self energies. There is one 
class of self energies that is obvious in the present formalism. When 
be, = bss. —1 there is a part joined by only two photon lines. 


RIASSUNTO (*) 


Una formula concisa per un elemento della matrice S è stata data da CAIANIELLO 
in termini delle particelle reali che appaiono nel processo. Questa formula contiene 
in aggiunta ai contributi dei diagrammi che descrivono le interazioni mutue di tutte 
le particelle che interagiscono fra loro, i contributi di processi disgiunti che avvengono 
simultaneamente. Nella presente nota si dà un metodo per escludere pezzi sconnessi 


non desiderati. Come esempio si esegue la rimozione di transizioni vuoto-vuoto nella 
elettrodinamica quantistica. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 12 Aprile 1960) 


Summary. — The elastic scattering of electrons by a proton with an 
intrinsic electric dipole moment (e.d.m.) is caleulated in first Born approx- 
imation. The cross section is quite sensitive to the e.d.m. interaction 
at high energies and large angles. It therefore seems likely that high 
energy electron-proton scattering may give valuable information of the 
value or upper limit of the proton e.d.m. coupling strength. Certain 
spin correlation effects between the recoil proton and the initial electron 
arising from the e.d.m. interaction are also discussed. In the last section 
the calculation of the ionization power of a neutrino due to its intrinsic 
e.d.m. is briefly outlined and compared with the experimental results 
of Cowan and Reines. 


1. — Introduction. 


It has been shown (2) that if space reflection (P) or time reversal (7) 


invariance holds, the existence of electric dipole moments (e.d.m.) is impos- 
sible. Therefore, detection of an e.d.m. would mean the non-invarianee of 
physical laws under the P and T operations. 


However, RAMSEY (*) has pointed out that if magnetic monopoles exist 


electric dipole moments are possible without violating the invariance under 
the combined inversion PC, where C represents electric charge conjugation. 


) This work was supported by the U.S. Air Force Office of Scientific Research. 
) L. LanpauU: Nucl. Phys., 3, 127 (1957). 

2) T. D. Lee and C. N. YANG: BNL 448(T-91) (1957). 

) N. F. Ramsey: Phys. Rev., 109, 225 (1958). 
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Such possibility would imply an additional degeneracy for elementary par- 
ticles, and because of the Pauli principle would change the form of the periodic 
table and the nuclear shell structure. It is therefore very unlikely that the 
electron, proton and the neutron have such a degeneracy. 

In recent years upper limits for the intrinsic electric dipole moment ot 
several of the elementary fermions have been determined by different techni- 
ques. The lowest upper limit for any e.d.m. has been determined for the 
neutron. SMITH, PURCELL and RAMSEY (4) obtained by a magnetic resonance 
method a value of 5-10-2e-cm, where e is the charge of the electron. Very 
recently, BERLEY and GIDAL (5) have found that the e.d.m. of the muon must 
be less than 2-10-!8e-cm. This has been accomplished by deflecting an ini- 
tially longitudinally polarized muon beam by a strong magnetic field and putting 
an upper limit on the polarization component perpendicular to the final mo- 
mentum. 

Upper limits for the e.d.m. of the electron have been founa by various 
methods. FEINBERG (°), SALPETER (7) and SACHS () have obtained an upper 
limit of 1-10-1*e-cm for the e.d.m. of the electron by calculating the effect 
on the Lamb shift of hydrogen-like atoms due to the hypothetical e.d.m. 
interaction of the electron and comparing it with known experimental results. 
Until now, the precession technique (*!°) has yielded the lowest upper limit 
for the e.d.m. of the electron: 3:10-%e-em. By a quite different technique, 
viz. high energy Coulomb scattering of electrons by spinless targets (1), an 
upper limit for the e.d.m. of the electron of 8-10- e-em has been obtained (1). 
It should by pointed out, however, that the e.d.m. of the electron may be 
considered momentum-transfer dependent. Therefore, as the last mentioned 
experiment is performed at a much higher energy, it does not necessarily yield 
the same information as the two other ones. 

STERNHEIMER (1°) has discussed the effects of a proton electric dipole moment 
on the energy levels of the hydrogen atom and obtained that the e.d.m. of 
the proton must be less than 1.3-10-!* e-cm. 

In Section 2 of the present paper, we discuss in first Born approximation 


(4) J. H. SMITH, E. M. PurcELL and N. F. Ramsey: Phys. Rev., 108, 120 (1957). 
(5) D. BerLEY and G. GiDAL: Phys. Rev. in press. | 

(5) G. FEINBERG: Phys. Rev., 112, 1637 (1958). 
(7) E. E. SALPETER: Phys. Rev., 112, 1642 (1958). 
(8) M. SacHs: Ann. Phys., 6, 244 (1959). 

(°) 

( 0 


R. L. GARWIN and L. M. LEDERMAN: Nuovo Cimento, 11, 776 (1959). 


(11) B. MARGOLIS, S. ROSENDORFF and A. SIRLIN: Phys. Rev., 114, 1530 (1959). 
(G2) Tel, E KENDALL and BURLESON: Kiev Conference on High Energy Physics (1959). 
(18) R. M. STERNHEIMER: Phys. Rev., 113, 828 (1959). 
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the elastic scattering of electrons by a proton target with an intrinsic e.d.m. 
It seems likely that high energy electron-proton scattering may give valuable 
and independent information on the value or upper limit of the proton e.d.m. 
interaction. As in the electron case, the proton e.d.m. coupling strength is 
to be considered momentum-transfer dependent. Therefore, such a high energy 
experiment would shed light also on this particular point. Certain spin cor- 
relation effects between the recoil proton and the initial electron arising 
from the e.d.m. interaction are discussed also in this section. In Section 3 
we calculate the ionization power of a neutrino with an intrinsic e.d.m. and 
compare it with known experimental results. 


2. — The e.d.m. of the proton. 


Let us consider the elastic scattering of an electron by a proton having a 
total magnetic moment (1-+K)e/2M, and an e.d.m. 2e/2M. K is known to 
be equal to 1.79. The matrix element for the first Born approximation in 
momentum space may be written in the form: 


oer ie 
(1) M = — An ÿo(da; 41) Tr ACE 
where the proton charge-current density is equal to 


) lao, do) = eT (ae) fy" Fula?) + ggg (EP) + ABQ? )y5) ofa] Ua) - 


bo 


( 


U(q,) and U(q;) are free particle spinors describing the initial and final proton 
states, ¢, = (d2.—), =—(P2—P.), is the four-component momentum-energy 
transfer (14), and q? = 4H? sin? (0/2)/(1+(2E,/M) sin? (0/2)), where 0 is the scat- 
tering angle in the lab. system, and £, the initial electron energy, assuming 
always E,> m,. The three functions Y,(4*) (à = 1, 2, 3) are form factors asso- 
ciated with the charge, the anomalous magnetic moment, and the e.d.m. of 
the proton, respectively. The electron charge-current density is, as usually, 
given by 


(3) Jeu (Pa; Di) = €U(Ps)y, Pi); 


(14) We have set A=c=1. Furthermore, we use the definitions 


y Ù 1 
NN A i 
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where u(p,) and u(p,) are the appropriate initial and final free electron spinors. 
Bq. (1) is invariant under the © operation, but not under the P and T ope- 
rations. Taking advantage of the fact that U(q,) and U(g») are solutions of 
the free Dirac equation, eq. (2) may be reduced to: 


(2') (qe, 1) = e U (qn) [y" F: — DÉF, — Di F,)U (qi), 


whe re 


A : 
Di = (nd) My) NDS (Gate) ys 


Summing over the initial and final spin states of both particles, the diffe- 
rential cross-section in the lab. system is readily obtained: 


(4) 6/Op, = 1 + (AF3)? (q/2.M)?/R(q, 0). 


or, is the well known « Rosenbluth scattering» cross-section by a physical 
proton: 


| 1 
(2) 1 GEM) sin? (0/2) OG 9): 
mes 
(6) R(q, 0) = Fi + (q/2M)? |2(F, + KF,)? tg? 9 RGB) 


where ©, is the Mott scattering cross-section from a spinless point nucleus. 

A priori, there is no reason why the three form factors F should have the 
same functional form. Yet, so far, no significant difference between F, and F, 
has been detected experimentally (1°), and all electron-proton scattering exper- 
iments up to 550 MeV are consistent with the assumption F, a F,. This 
however, could be a misinterpretation of the experimental data, as it is 
based on the a priori assumption of vanishing e.d.m. interaction of the proton. 
Nevertheless, we will, for the time being, stick to F, — F, even for a non- 
vanishing e.d.m. interaction. Later, at the discussion of possible polarization 
effects, it will become clear that this is no more necessary. Of course, nothing 
can be said about F,. With F, = F,, it is a relatively easy task to determine 
for a given value of ¢ an upper limit to the e.d.m. term JZ, (19). All one has 
to do is to compare the values of o/o,, at two different angles with the same 
value of g. As a numerical example, if AF,/F, = 5 (according to the Lamb 


(5) R. HoFsTADTER: Ann. Rev. Nucl. Sci., 7, 231 (1957). 


(16) See also ref. (11), Sect. 2 for a similar discussion about the e.d.m. interaction 
of the electron. à 
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shift calculation of Ref. (18) À < 13) the ratio of eq. (4) at 0= 60° and 0-- 140° 
for q/M — 0.5 would be 1.75. Accurate elastic electron-proton scattering ex- 
periments have already been performed (#) up to energies of 550 MeV and 
angles of 135°. As this example shows, such an experiment seems to provide 
a sensitive method to put an upper limit on the e.d.m. interaction. As this 
hypothetical interaction may be regarded q-dependent, the experiment should 
be performed for different values of q. 

Additional information on the e.d.m. interaction in particular, and non- 
invariance of the electromagnetic interaction under space reflection and time 
1eversal operations in general, could be obtained by detecting certain spin 
correlation effects between the recoil proton and the initial electron. In order 
to see that, we have to calculate the cross-section o(q,¢;¢.,¢,), where &, 
describes the initial electron polarization and €, the final proton polarization. 
This cross-section is readily obtained by standard trace technique, if use is. 
made of the known (1) relativistically covariant spin projection operator: 


PG.) = (1 + tys8,) , 


where 


CONE Pi(P1%) 
m(H, + m) 


with an analogous operator for the proton. The result, assuming again 
E,>wm,, can be written in the form: 


Ou 
1+ (2#,/M) sin? 6/2 


AG CR) R(q, 0) + (4/2 M) (AF)? + 


+ (AF,)(F, + KF,)(d/2M)° sec?=(n,"C.)(n. Xn,:GY+E S(q, 9; Cl 


bo | D 


where n, and n, are units vectors in the direction of the initial and final electron. 
momenta p, and p,, and R(g,0) is defined by eq. (6). As is obvious from 
eq. (7) the term proportional to À arises from the interference between the 
e.d.m. interaction and the Dirac plus Pauli interaction. There is no pola- 
rization dependent 2? term. The function S has not been calculated explicitly, 
but is certainly /-independent. Detection of the term proportional to À in 


(7) H. A. TOLHOEK: Rev. Mod. Phys., 28, 277 (1956). 
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eq. (7) would imply the non-invariance of the electromagnetic interaction 
under P and 7 operations. This is obvious, as this term follows directly from 
the e.d.m. interaction which, as has already been pointed out, violates both 
P and T invariance. 

More generally, any interaction Hamiltonian which violates P as well as 
T will in first Born approximation (i.e. neglecting final state interaction) give 
rise to terms in the cross-section with similar structure as the 2 term in 
eq. (7), namely, terms which involve an even number of spins and an odd 
number of momenta. This follows immediately from a well known theorem 
proved by LEE and YANG (2). Of course, nothing can be said about the energy 
and angular dependence of these terms in general. This is possible only, when 
the interaction Hamiltonian is given explicitly (e.g. the e.d.m. interaction). 
It is easy to show that terms of this kind, and only of this kind, (in first Born 
approximation) contribute to polarization of the recoil proton which is per- 
pendicular to the scattering plane, provided the initial electron beam is lon- 
gitudinally polarized. The simplest representative of this type of terms is 
(q:°¢.xC,); another one is (p,-C,)(piX p.°S,) which appears in eq. (7), etc. On 
the other hand, terms which do not violate P and T invariance are even in 
both momenta and spins, like (€,:¢,) or (p.Xp.):(G.X%,), ete. It is evident 
that these terms cannot contribute to perpendicular proton polarization (for 
initial longitudinal electron polarization). Thus, the function S in eq. (7) 
which ineludes only P and 7 invariant terms, has no influence whatsoever on 
the perpendicular polarization of the proton. 

Now, if the initial electron beam is longitudinally polarized, the polariz- 
ation of the recoil proton perpendicular to the scattering plane follows in an 
obvious way from eq. (7). It is given by the expression: 


_2(4P,)(F + KF,)(q/2M)? tg (6/2) 
R(q, 0) + (q/2M)(AB5)? 


(8) P, 


As a numerical example, consider an initial electron energy of E,/M = 0.25 
and 4F;/F.=1, then for a scattering of 120°, the recoil energy of the proton 
would be 7,/M = 0.07 with a perpendicular polarization of 12% (assuming 
again F,= F,). Thus, a substantial proton polarization could be expected for 
a relatively small value of the e.d.m. coupling strength. There is an addi- 
tional transverse polarization of the proton to be expected which comes from 
the ordinary (P and 7 invariant) proton-electron interaction in higher Born 
approximations. However, this effect, which arises from a spin-orbit coupling, 
is estimate to be of the order of 0.75%. It is therefore negligibly small. The 
transverse polarization of the proton could be analysed by proton-proton or 
proton-nuclei scattering. Of course, such an experiment would involve the 
production of a high energy longitudinally polarized electron beam, which may 
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be a difficult task at the present time. Finally, it should be pointed out, that 
the measurement of the differential cross-section (4) as a function of the scat- 
tering angle 9 (for a constant value of q) plus a measurement of the polari- 
zation (8) provide a unique solution for the three unknown functions F,, F,, 
and ÀF,. 


3. — The e.d.m. of the neutrino. 


As for the e.d.m. of the neutrino, it is most easily determined by consi- 
dering the ionization power due to the interaction of the e.d.m. with the 
electrons. If we take into account only those electrons, whose kinetic energy W 
is much larger than the ionization potential of the atom, it is easy to show 
that, to a good approximation, momentum is conserved in the collision (18). 
Therefore, the electron may be regarded a free particle, and its total energy 
after the collision will be V/ (P2— p.)} +m. It follows, that the matrix ele- 
‘ment responsible for this interaction is given by eq. (1), provided we put there 
50: 

Then, after summing over the spin states of the neutrino and the electron, 
we obtain for the energy distribution of the secondary electrons the expres- 


sion: 
(9) DW)dW= Ax(e2,F,)} (1 —W/H)aw/w, 


where 4, is the e.d.m. of the neutrino and Æ its initial energy. Eq. (9) is 
obtained after taking the limit of vanishing neutrino mass. 

Now, exactly the same expression is obtained for a magnetic dipole inter- 
action. This has been calculated a long time ago by BETHE (1%). CowAN and 
REINES (2°) have performed an experiment in order to determine an upper 
limit for the magnetic dipole of the neutrino. We can directly take over their 
result: 


i PAT IGT 


(18) For a more thorough discussion of this problem, see H. A. BetHE: Handb. d. 
Phys., Bd. 24(1) (Berlin, 1933), p. 491. 

(19) H. A. Berne: Proc. Cam. Phil. Soc., 31, 108 (1935). Actually, our formula 
differs from Bethe’s expression by a factor m/(m+W). This is so, because we use a 
slightly different approximation in the derivation of the density of final states. 

(20) C. L. Cowan and F. REINES: Phys. Rev., 107, 528 (1957). 


17 - Il Nuovo Cimento. 
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It is perhaps worth-while mentioning that a hypothetical e.d.m. due to the 
weak interaction would most probably be smaller than (10729 — 10~**) e-cm. 


The author is much indebted to Prof. H. PRIMAKOFF and Dr. S. TANI for 
many helpful discussions. 


RIASSUNTO (*) 


Si calcola in prima approssimazione di Born lo scattering elastico di elettroni su 
un protone con un momento elettrico dipolare (e.d.m.) intrinseco. La sezione trasver- 
sale è molto sensibile alle interazioni e.d.m. di alte energie ed angoli grandi. Sembra 
quindi verosimile che lo scattering protone-elettrone di alta energia possa dare valide 
informazioni sul valore o sul limite superiore della forza dell’accoppiamento e.d.m. 
del protone. Si discutono anche alcuni effetti di correlazione di spin fra il protone di 
rinculo e l’elettrone iniziale nascenti dalla interazione e.d.m. Nell’ultima sezione si 
delinea brevemente il calcolo del potere ionizzante del neutrino dovuto al suo e.d.m- 
intrinseco e si paragona ai risultati sperimentali di Cowan e Reines. 


(*) Traduzione a cura della Redazione. 
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(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciaia 
dalla Direzione del periodico ai singoli autori) 


Search for Antineutron Annihilation Events in Nuclear Emulsions. 


Tsar-CHÜ, M. Moranp, G. BourLET, C. SIMONIN and D. SCHUNE 


Laboratoire de Physique Enseignement, Faculté des Sciences - Paris 


(ricevuto il 13 Aprile 1960) 


Shortly after the production of anti- 
protons (1) by the Bevatron, antineu- 
trons (?*) were created by charge ex- 
change of antiprotons in matter and 
detected by electronic means. However, 
up to now, no antineutron annihilation 
event has been observed in nuclear 
emulsions. We followed more than 
250 antiproton tracks through the plates 
contained in a stack of K-5 emulsions 
exposed to the enriched antiproton beam 
of 740 MeV/c. About half of the anti- 
protons had been annihilated in flight; 
among them we observed seven « stops ». 
In two cases a small nuclear recoil was 
possibly present. An antiproton stop 
should indicate either a charge exchange 
of the antiproton into an antineutron or 
an annihilation star without any charged 
prongs. Antineutrons can also be created 
in charge exchange stars produced inside 


(1) O. CHAMBERLAIN, E. SEGRE, C. WIE- 
GAND and T. YPSILANTIS: Phys. Rev., 100, 
947 (1955). 

(?) B. Cork, G. R. LAMBERTSON, O. PIc- 


CIONI and W. A. WENZEL: Phys. Rev., 104, 
1193 (1956). 
() J. BUTTON, T. ELIOFF, E. SEGRE, 


H. M. STEINER, R. WEINGART and T. YPSI- 
LANTIS: Phys. Rev., 108, 1557 (1957). 
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the stack or in antiproton interactions 
occurred outside. In the course of 
scanning, we happened to find a star 
produced by a neutral primary, with 
more than two minimum-ionizing tracks 
emerging from the star. We suspected 
it to be an antineutron annihilation 
event. A systematic search was made 
afterwards along the prolongation of the 
stops and their neighborhood and in the 
proximity of the charge exchange stars. 
As a result, four other antineutron 
annihilation events were found. 

Fig. 1 shows a projection drawing of 
one of our antineutron annihilation 
events observed. The star has 17 prongs; 
it involves an associated production of 
K7-K°-mesons, with the K° decaying as 
6° into two r°’s immediately after the 
production and each 7° in turn decaying, 
respectively, into a Dalitz pair and a 
y-ray. In our experiment, the antiproton 
beam arrives at the stack from the 
upper side of Fig. 1 and follows a di- 
rection approximately parallel to that 
of prong 15. After 1.5 cm, the anti- 
proton comes to a stop in the emulsion. 
This occurs at 50 um from the glass of 
the second plate. It has a kinetic energy 
of 240 MeV, and a dip angle of — 4° 20’ 
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in the direction from glass to air. The 
antineutron annihilation star and the 
stop appear in the same plate 1.95 cm 


31 


25 pm 
16 RR) © DRE 
ayy DS 
f 
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apart. The probability for an accidental 
coincidence of the two events is less 
than 10 *. The arrow shown in Ties, 1 
represents the direction of the anti- 
neutron and the dip angle of the anti- 
neutron is estimated at —1° 20’. The 


resultant momentum of the prongs of 
this star shows a major unbalanced 
component in the direction of the anti- 


{18 


14 « 


— The antineutron annihilation star. 


neutron. The two electron pairs of the 
§°-meson are formed by the four min- 
imum-ionizing tracks (prongs I 2, 5, 
and 14): the space angle between the 
first pair is 33°; that between the second 
is 46°. Electrons are identified by ioni- 
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zation and multiple scattering. Distor- 
tion is considered and eliminated for 
steep tracks. Prong 14 is doubtless an 
electron, as it degenerates into a low 
energy electron in the middle of the 
emulsion. The K~-meson (prong 7) is 
emitted with an energy of (124+11) MeV. 
After 1.63 cm across 13 plates, it makes 
an inelastic scattering involving an 
energy loss of 55 MeV and a deflection 
of 22° 30’. 1.3 mm further on, it inter- 
acts again and produces a small star 
containing two protons and a nuclear 
recoil with a total visible energy of 
20 MeV. The mass of the K -meson 
measured by ionization and scattering 
is (944 +90) m, before, and (940 +200) m, 
after, deflection. The three knock-on 
protons (prongs 8, 11, and 16) are 
identified by the increase of grain den- 
sity along the trajectory. Prong 11 
produces a small star of 66 MeV after 
9.6 mm. After 2.45 cm, prong 16 also 
produces a star of two protons: one of 
11 MeV is ejected nearly perpendicular 
to it and the other of 79 MeV at an 
angle of about 10°. Prong 17 is a low- 
energy electron. The other prongs (3, 6, 9, 
10, 12, 13, 15, and 18) are all low-energy 
evaporation protons, except prong 4 
which could be an « particle: all come 
to rest in the emulsion. Assuming a 
binding energy of 8 MeV for a proton 
and 2 MeV for an « particle, we have a 
total energy of 1.44 GeV for the whole 
star (Table I). The total available energy 
of the antineutron and the nucleon 
before annihilation is about 2 GeV 
Therefore the remaining energy, 0.6 GeV, 
is shared among the two y-rays of the 
rs and the knock-on and evaporation 
neutrons. Through charge conservation, 
at least one positive pion is emitted 
together with the K-mesons in the 
annihilation of the antineutron with a 
neutron and at least two positive pions 
are emitted in an antineutron-proton 
annihilation The positive pions are 
absorbed afterwards and produce the 
knock-on and evaporation nucleons. 
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TABLE I. — Nature, direction and energy 
of the prongs of the antineutron annihi- 
lation star. 


Prong Dip Energy 
no. SEMEN | (degree) (MeV) 
1 PERSI 90+ 9 
2 e | 22 67-+10 
3 p — 43 8.7 
+ x | 8 5.3 
5 e | 52 80430 
6 p 59 15.6 
7 ica PEN ZII 
8 p 33 | 106+10 
9 p 2 8.6 
10 p 11 zi 
11 p 45 91+10 
12 p — 32 4.7 
13 p 9 2.0 
14 e 46 73+18 
15 p 30 13.2 
16 p 22 | 168-410 
17 e — a 
18 p 23] 6.4 


It is very unlikely that the anti- 
neutron events mentioned in this paper 
could be produced by cosmic-ray neu- 
trons. High energy protons and neutrons 
in cosmic rays have a sharp zenithal 
distribution in lower altitudes; 94% of 
them are contained in a cone having a 
half opening angle of 45° and an axis 
in the vertical direction. For one day’s 
flight, the number of stars produced by 
fast cosmic ray neutrons outside the 
cone should be 0.01 instead of 5 as in 
our scanned volume. Furthermore, if it 
were a cosmic-ray event, the estimated 
energy for a primary neutron should be 
hardly over the threshold value for the 
production of K-mesons. Again, all the 
secondary prongs of the star should be 
collimated along the direction of the 
primary, instead of the observed distri- 
bution. 

The four electrons of the star cannot 
be interpreted as the decay products 
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from a single 7°. Their resultant mo- 
mentum is p—211 MeV, their total energy 
E=313 MeV, the mass of the x° calcula- 
ted from m,e° = (E? —p%c?)? will be 
230 MeV. This impossible supposition 
is. quite evident from their wide an- 
gular separations (Fig. 1). The next 
alternative is to interprete them as the 
products from two independent 7°’s. 
This is again excluded by the insufficient 
energy available from antineutron anni- 
hilation. 

It is still too early to predict that 
complete similarities exist between anti- 
proton and antineutron annihilation 
events. But, for the moment, it is 
interesting to point out the excess of 
pions of different signs. The study of 
antiprotons has furnished a large number 


TSAI-CHU, M. MORAND, G. BOURLET, C. SIMONIN and D. SCHUNE 


of pp and pn annihilation events; while 
the study of antineutrons will give nn 
and np annihilation events. An excess 
of negative pions is observed in anti- 
proton annihilation events; this follows 
directly from the conservation of charge. 
In contrast with the antiproton annihi- 
lation events, an excess of positive pions 
will be present in antineutron annihila- 
tion events. 


We wish to express our sincere thanks 
to our colleagues in Berkeley and Bristol 
for the exposure and preparation of the 
emulsions. We also want to thank 
Mr. R. CuNyGHAM-BRown, Mrs. M. Cat- 
LIAUD, Mrs. L. van Rossum, and 
Mrs. H. Y. CHANG for technical assistance. 
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Some Remarks about the Spectral Representation 
of the K* —37x Decay Amplitude. 


S. FUBINI 


CERN - Geneva 
Istituto di Pisica dell’ Università - Padova 


R. STROFFOLINI (*) 


CERN - Geneva 


(ricevuto il 15 Giugno 1960) 


The effect of pion-pion interaction on the K+—3z spectrum has been recently 
investigated (1:2). 
The following representation has been proposed for the decay amplitude 
A4(s,, 89, 83): 
1 g(a) dx 1 g(x) da 1 g(x) de 
(1) Alig, ty) = of PO | = | 
ma} x—s —ie ala—-s—ie a 
au au Au? 


where 


(K and x, denote the momenta of the K-particle and of the pions). (*) 

From the representation (1) an integral equation has been derived which can be 
approximately solved by taking into account only the effect of the lowest singu- 
larities. This has enabled some information to be obtained from the K spectrum 
about the s-wave pion-pion interaction. 

Representation (1) can be simply understood in the following manner. Let us 
consider the unphysical case in which the mass of the K-meson is of the same order 
of magnitude as the pion mass; in such a case we have to deal with the scattering 


process 


(2) Kir-r+r. 
(*) On leave from: Istituto di Fisica Teorica dell’Università di Napoli. 

(1) N. N. KHuri and S. B. TREIMAN: Princeton preprint. 

(3) R. F. Sawyer and K. L. Watt: University of Wisconsin preprint. 

(*) In order to simplify the notation we shall disregard here the complications due to the 
charge of the pions. 


* 
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Perturbation theory then suggests the validity of the Mandelstam representation 
for the scattering amplitude 4° of the process (2) 


À res i 
(3) A8(8,, 8, 83) ci na: > 


a(x, Se i 1 o(a, y) de dy 
+=] (a — 8, —te)(y —8—ie) | m2} | (& 8; —ie)(y —8, —ie) © 


It has been shown by M. Cini and one of us (*) that under very general arguments 
on the location of the singularities, the Mandelstam representation can be trans- 
formed into the approximate form, valid for |s;|< 164? 


1 n g(a) da 1 x) de 1 g(x) da 
(4) A*(81, 89: 83) = ={ 2 2 Git i 3 Sar Da ame, 
nu} o—s,—te 2] @—S,—te nf as, —tle 


au? aus 


the spectral function g(x) can be obtained by computing the imaginary part of 45 
for s, > 4u?. Since only the denominator of the first term can vanish, we get 


(4) Im A%(s,, 82, 83) = 9(81) - 


On the other hand, the theorem of final state interaction (5) as applied to the 
process (2) gives 


(5) g(s) = A,(s) exp [— iô(s) sin 6(s) , 


where ô(s) is the s-wave pion-pion phase shift and A,(s) is the s-wave part of the 
amplitude A. Eq. (5) follows from the fact that the K-meson is here connected. 
to pions only through a weak interaction. . 

Now it is easy to observe that the representation for the decay amplitude can 
be written in the form (1) by assuming that eq. (4), together with the expression (5) 
of g(x) valid for the scattering process, can be extended also to the actual value 
of Mx corresponding to the decay process. 

The possibility of justifying such continuation has been studied by G. Bon- 
NEVAY (6) who has investigated in detail the analytical properties of the decay 
amplitude in the mass of the K-meson. 

We point out here that in the representation (1) for the decay process a new 
feature arises. All three denominators in eq. (1) can vanish simultaneously for 
values of s,, ss, 8; in the physical region. | 

If we assume the validity of (1) and (5) for the decay process, then the spectral 
function g(s) must be a complex quantity (7). This can be easily seen in the fol- 
lowing manner. Let us assume g(s) to be real, then eq. (5) tells us that 


Im Aj(s) = g(s) . 


€) M. Cini and S. FuBINI: to be published in Ann. of Phys. 

(5) S. FUBINI, Y. NAMBU and V. WATAGHIN: Phys. Rev., 111, 329 (1958), Appendix II. 

(9) G. BONNEVAY: private comunication. 

(*) This result has also been obtained independently by G. BONNEVAY and A. STANGHELLINI. 
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This result is in contrast with what one would obtain for g real by taking the 
imaginary part of eq. (1) 


Im A°(8,, 89, 83) = 9(51) + 9(52) + 9(83) . 


We see therefore that the representation (1) for the decay process has a completely 
different mathematical structure than the similar equation for scattering. In order 
to understand the physical meaning of the imaginary part of g let us write 


g(s) — g*(s) 
9j ì 


(6) o(s) = 


Let us take the imaginary part of Ad 


oo (co) 


1 xv) dr x) dx à o(x) da 
(7) Tm A = gle) +06) + 06) + 2] | ee +f el dj sei | 


x —s,+ de Z— 8+ die D — 8+ de 
au? au? 4p? 


We note that the term containing o(x) must vanish for Mg 3. 

Now, using the unitarity of the S-matrix, and invariance under time reversal 
and under the Lorentz transformations (5), one obtains the following general ex- 
pression for the imaginary part of the decay amplitude: 


(8) Im A(z, ; 72, 363) = fan. Ge GINS ni, my, my) Aa, n, 13) 
. O(a + u?) O(a + wu?) dn + a?) 0111) 0x2) (og) dia dis dal 


where T is the transition matrix element from the state x,,,, 3 to the state 
M, 7% Te: By comparing eq. (7) with eq. (8) it is possible to recognize the con- 
tributions to eq. (7) as due to the effect of correlations between 2 pions or 3 pions. 

We thus split the matrix element 7 into a connected and a disconnected part 


(9) T= daly Too 
where 
r r 
(10) To = 2008 (10 — TO) tg, 733 Has Hy) + 2 dt 3 My) + 
+ 273.9° (703 — Ts), Ta ; mi; I) 
LT, 7053 1, x) is the pion-pion scattering amplitude, and 7’, (1, 7, 73; Ti; Ty, Ty) 


is the matrix element representing the sum of all those graphs in which all three 


pions are involved. 

Let us substitute eqs. (9) and (10) into eq. (8). It is easy to verify that the 
contribution of the disconnected part 7, (if one keeps only the pion-pion s-waves) 
reproduces exactly the term > g(s;) in the r.h.s. of eq. (7) where g(s) is given by eq. (5). 

i 


The contribution of the connected part 7, can therefore be idertified with the 
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term > J o(x) dx/(a —s,+ ie), such a term, as previously stated, vanishes where Mg 
decreases to 3y. 

We see therefore that the imaginary part o(x) of the spectral function g(x) is 
directly connected with the 3 pion transition amplitude. Thus eq. (1) might be 
a tool for constructing the 3 pion amplitude in terms of the pion-pion phase shifts. 


* OK OK 


We wish to thank Dr. D. Amati, G. BONNEVAY, A. STANGHELLINI and B. VITALE 
for very useful comments and criticism. We are particularly grateful to Dr. G. Bon- 
NEVAY for having shown us his unpublished work; many questions have become 
clearer from comparison with his results. 
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An Improved Version of an Ergodic Theorem in Quantum Mechanies. 


G. M. ProsPERI and A. ScortI 


Istituto di Scienze Fisiche dell'Università - Milano 
Istituto Nazionale di Fisica Nucleare - Sezione di Milano 


(ricevuto il 21 Giugno 1960) 


In this letter we intend giving an improved version of a theorem established 
in a previous paper (1) (this paper will be referred to in the following as A). The 
new version is stronger than the older one and conceptually more satisfactory. 

The fundamental formula proved in A is the following: 


N a oo SN 
(A.10) PLU (la, m0) > 81} < exp|— Vout ++ log], 220.38, 


where #,(a, y,) denotes the time fraction during which the quantity 


(u,(t) — s,/8)? 
A,(t = —————- , 
v(t, Po) He 


(square deviation of the quantum mechanical probability of occupation of the v™ 
cell from the microcanonical value s,/S) remains greater than a and P denotes the 
measure of the set of initial state vectors y, for which even one of the d,(4, yy) 
alone is greater than #. 

Formula (A.10) will here be substituted by the following one, which expresses 
our new version of the theorem: 


CAT 


6 © do SN 
(1) PLO(G, Bo) > 0} <exp|—%V8an0+ % + bog ; 


where now (a, yo) (*) represents the time fraction during which even one of the 
A,(t, Yo) alone is greater than a and P is the measure of the set of initial state 
vectors y, for which # is greater than #. 


(@) G. M. ProsPERI and A. Scorri: Nuovo Cimento, 13, 1007 (1959). 
(*) Notice that the symbol 8 was used in the last section of paper A with another meaning. 
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To deduce formula (1) let us start by considering the following relation 
N n 
(A.6) PU [A,(t, Yo) > a]} < exp [—x 4/8, à + x + log SN]. 


Let us, then, call Q the hypersphere of the vectors of unit norm of the energy 
shell and (0, 7) the time interval during which the evolution of the system is con- 
sidered (assumed to be fixed, and, at least in the first instance, finite). 

Proceeding now in a way analogous to that followed in lemma (A.7) let C(Yo> 4) 
and D(t,a) be the set of values of the time interval (0, 7) and the part of © for 
which at least one of the following relations is verified 


(2) APE CE (= SCENDE 
for y=%, and t=t respectively. 


Then by definition we have 


Di di 4 = do 
(3) CIUTE a) -/3 E) PIU [Ay(t, Yo) > a)} -[£ , 


Co, 2) Dt, a) 
where do denotes the element of area of the spherical hypersurface © and £ its 


total area. 
By known theorems on integration we obtain 


Ti 
do dt dt [do 
(4) = —. 
Di ID ft Zi 
Lo) C(Yo> 4) 


0 D(t, a) 


And thus, from (A.6) 


7 
‘ do diga 
(5) da, Yo) = T P{U [A,(t, Yo) > a]} < exp [— x 4/8 à + x + log SN]. 
Li v=1 
a ò 


From which (1) follows immediately by noticing that the left hand side is clearly 
greater than 9P{Ha, po) > 8}. 
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Remarks on Neutral Pion Photoproduction in the High Energy Region. 


C. PELLEGRINI and G. STOPPINI 


Laboratori Nazionali del Comitato Nazionale per le Ricerche Nucleari - Frascati 


(ricevuto il 30 Giugno1960) 


J. J. SAKURAI (!) proposed the measurement of the recoil proton polarization 
in the process 


(1) alpino 


as a possible means of assigning the parity of the second resonance recently ob- 
served in the pion nucleon scattering and photo-production (?). In the energy 
region between the first and the second resonance, namely around Æ,—600 MeV, 
the polarization predicted was as high as 80% for negative parity and zero for 
positive parity. 

The experimental results (*) show indeed a high polarization, though less than 
80%, in the energy region (550--850) MeV (Table I) and this can be taken as a 
support to the negative parity assignment. 

We already pointed out (4) that the presence of a small s-wave amplitude can 
appreciably change Sakurai’s predictions, so to make the interpretation of the 
experimental results not completely unambiguous, if the large experimental errors 
are also taken into account. 


(3) J. J. SARURAI: Phys. Rev. Lett., 1, 258 (1958). 

(2) J. W. DE WIRE, H. E. JACKSON and R. LITTAUER: Phys. Rev., 110, 1208 (1958); P. C. STEIN 
and K. C. Rogers: Phys. Rev., 110, 1209 (1958); M. HEINBERG, W. M. MCCLELLAND, F. TURKOT, 
W. M. Woopwarp and D. M. Zipoy: Phys. Rev., 110, 1211 (1958); F. P. Dixon and R. L. WALKER: 
Phys. Rev. Lett., 1, 142 (1958). 

(3) R. QueRzOoLI, G. SALVINI and A. SILVERMAN: private communication, to be published. The 
‘errors given by these authors includes all the corrections. P. C. STEIN: Phys. Rev. Lett., 2, 473 (1959); 
L. BERTANZA, P. FRANZINI, I. MANNELLI, V. Z. PETERSON and G. V. SILVESTRINI: private com- 
munication; J. O. Manoy, J. I. FRIEDMAN, H. KENDALL, A. MANFREDINI, V. Z. PETERSON and 
G. SALANDIN: private communication. The two last are preliminary results, submitted to Rochester 
Conference on High Enegry Physics (1960). We thank the above mentioned authors for the com- 
munication of their data. 

(4) GC. PELLEGRINI and G. STOPPINI: reported by G. BERNARDINI at the Kiev Conference on 


High Energy Physics (1959). 
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TABLE I. — Polarization of the recoil protons emitted at 90° in the em. 
R. QUERZOLI, ta L. BERTANZA I. 0. MALOY 
E,(MeV) | G. SaLvin and PCS ERIN Bart | PT 
A. SILVERMAN 
550 — 0.30 +0.12 — = 
560 0.4 40.14 — = -- 
585 — — — 0.35 0.16 
610 0.63 10.27 = = | 
650 0.59 40.24 — — | — 
658 — | _ — 0.50 +0.18 
700 | 0.56 +0.11 0.59 +0.06 — | — 
720 == — 0.82-+0.29 | = | 
In the region (700+850) MeV preliminary results by R. QUERZOLI et al. seem to indicate 
a polarization yet equal to about 50%. 


As a matter of fact let us consider two models a) and b) for reaction (1) proposed 
respectively by R. R. WiLson (5) and R. F. PEIERLS ($) according to the transition 
matrices 


1 
(2) M, = 1H,(4)o-e — Fg Dl 4) + M, (3, $)]{2q(kxe)—i(o-kq-e—o-eq-k)}, 


vq 
1 
(3) M, = iB, (4)o-e Sa M,(%, 3){2q:(kxe) —i(o-ke-q —c'eq:k)} — 
. 1 9 
— 1B, (3, 3) SDAI dore), 


where £,(J, T) and M,(J, T) are the amplitudes for the absorption of an electric 
or magnetic multipole of order 1 leading to a final state with total angular mo- 
mentum J and total isotopic spin 7. 

o, e, k, q are respectively the spin of the proton, the polarization vector of 
the photon and the momenta of the photon and of the pion in the center of mass 
system. 

In model a) it is assumed that the second resonance is excited through a 
magnetic dipole absorption leading to a state T=}, J=3%, positive parity, while 
in model b) the multipole involved is supposed to be the electric dipole leading 
to the same state but with opposite parity. 

In both models £,(3) is the s-wave amplitude and M,(%, 3) is the multipole 
giving rise to the first resonance. 

Starting from (2), (3) and choosing as positive the polarization in the direction 
kxq, we obtain for the polarization of protons emitted at 90° in the center of 


(°) R. R. WILSON: Phys. Rev., 110, 1212 (1958). 
(5) R. F. PEIERLS: Phys. Rev. Lett., 1, 174 (1958). 
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mass system 


1 , 4 
(4) P, = = Im {ET (LM (3, 3) + ME, Dh, 
1 ‘i +k 1k x 
(5) Py = = Im {BG (ES, 3) + ABS, 4) M8, 8}. 


A is the isotropic term of the angular distribution for process (1) 


(6) DOr Per AE Jo Roc 
RE > COS U COS“ OU, POS — 
dQ o Di 


The term #,(4) introduced in (2) and (3) is necessary to obtain a non zero 
interference term at all energies in agreement with the experimental results (Fig. 1). 

Models a) and b) are not the only possible ones. Furthemore they are not able 
to explain facts such as the third resonance at about 950 MeV. However, we shall 
limit ourselves and consider only them. 


ile 1 1 4 IL 
200. 300 400 500 600 700 800 900 1000 
Ey(MeV) 


Fig. 1.- O W. S. McoDonanp, V. Z. PETERSON and D. R. Corson: Phys. Rev. 107, 577 (1957); 

O J. I. Verre: Phys. Rev., 111, 622 (1958); A K. BERKELMAN and J. A. WAGGONER: Phys. Rev., 

117, 1364 (1960); X R. M. WoRLOCK: Phys. Rev., 117, 537 (1960); G. E. MODESITT: * Thesis, 
unpublished; @ V. I. GoLDaNsSKy et al.: Nucl. Phys., 12, 327 (1959). 


We shall assume that they describe correctly the experimental situation around 
600 MeV. Thus, roughly speaking, the polarization below the second resonance 
is given by (4) and (5), while above, it becomes necessary to explain the large 
polarization, by introducing in the models an amplitude corresponding to the third 
resonance. 
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An alternative model, which gives respectively to the resonances the parity ++—, 
and explains all the polarization as due to the interference between the third and 
the first and second resonance has been proposed by L. MarsHALL and L. F. LAN- 
povirz (7). We can have a justification for not considering this model in the fact 
that up to 700 MeV there is no appreciable evidence for a cost 0 term. This term 
follows from the model owing to the fact that the third resonant state is a J=$ one. 

As was pointed out by J. 0. MaLoy (8) P, can reach a maximum value of — 31%. 
P,, when Æ,(4) is supposed to be negligible, can reach a maximum of 80%. 

If it is assumed that the addition of £,(3) does alter appreciably the value 
of P, we see that model b) is in better agreement with the measured values 
around 600 MeV. 

Anyway, due to the large experimental errors, it seems preferable to have a 
more stringent criterium of distinction between the two models. 

In fact it is difficult to obtain quantitative predictions from (4) and (5), be- 
cause, in spite of the simplification assumed it is impossible to derive from the 
experimental A, B and © all the complex multipole amplitudes. Actually in the 
energy region under consideration not even the multipoles phases are known. For 
this reason it seems that the experimental determination of the absolute polarization 
would serve more as additional experimental information useful for the extraction 
of the multipole amplitudes than for giving sharp answers of the kind required 
by Sakurai. 

The sign of the P in (4) and (5) is a priori undetermined because of the arbi- 
trariety in the sign of the multipole amplitudes. 

We want to show that the sign of P, and P, can be determined if B is assumed 
to be known and that, within the limits of reasonable hypotheses, only model b) 
allows to predict the sign experimentally observed. By writing 


His) = e, exp (te, |, Hi, ($, 4) = #, exp [26], 
Ma (és, +) = mM, exp lis, M(&, 3) = ms €xp io: |" 
we get 
ta Î IP = pi {ma sin (&33 —- 0) + M, sin (a3 — 21)} , 
| By = — 2e, {Ms cos (ass — 21) + my COS (dis — )}, 
| I = ene sin (0% j3 — %) + 4H, sin (og; — d)} 
D jo ae ce ee 
| B,= — 2my{e, cos (%3 — a1) + E cos (org — d)}. 


At Hy, ~ 450 MeV, m, and E, are negligible, respectively and it results from the 
experiment B>0. On the other hand cos (#3, —a,)<0(%). Then is has to be 
e,mg > 0. By choosing m, > 0 (all the results do not change if m,< 0 is assumed) 
we get ¢,>0. 

At By — 650 MeV, i.e. below the second resonance, where m, and e, are negligible 
respectively, B< 0 and 6<90°, x,3<90° so that cos (033 — 6) < 0, cos (a3 — 4) > 0. 


L. F. LANDOVITZ and L. MARSHALL: Phys. Rev. Lett., 3, 190 (1959). 
J. O. MALOY: private communication. 
B. PoNTECORVO: Proc. of the Kiev Conference on High Energy Physics (1959). 
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It follows E,<0, m, > 0 and 


iF Sal IP. (0) 


The polarization measurements give a negative value so excluding model a): 
on the ground of the above consideration. 

To obtain this result we made the hypothesis that e,, does not change sign 
going from low to high energy. Furthemore the s-wave term should be written 
E,(3)=e, exp |ix,]+e, exp [is]. As experimentally we know that «, >, we neg- 
lected the 7 —3% state. 

In the argument only for the s-wave part has it been assumed that the mul- 
tipole phase coincide with the corresponding scattering phase shift. The validity 
of this assumption has been justified by C. PeLLEGRINI and L. Tau (19) on the- 
ground of a quite resonable model. 

The choice E,<0 seems to be consistent with the observed 7° angular distri- 
bution. In fact in model b) 


A+ $C =|#,(4)|?+ 4e,H, cos (6 —«,), 


and at H.~ 700 MeV. C ~—A(") so A+(8)C <0. Because cos (0 —a,) > 0: 
this means that e, and Æ, have opposite Stan In the model a) A+8C =|#,(4)|?>> 0 
in contradiction with the experimental result © ~— A. It is worth being noted 
that the addition of the e, exp [ix3] term to the s-wave amplitude might change 
the sign of P, but not the sign of A+40, in the model a). Instead the same addition 
could have the effect of changing the sign of À +50 leaving the sign of P, unaltered, 
in model b). 

Furthermore it has to be noted that if E, < 0, e, > 0 the two terms in P tend 
to cancel each other so qualitatively explaining why the measured polarization 
is less than the value predicted by J. J. SAKURAI. 


We wish to thank Prof. M. Crt for very helpful discussions. 


(1°) C. PELLEGRINI and L. Tau: to be published in Nuovo Cimento. 
(1) K. BERKELMANN and J. A. WAGGONER: Phys. Rev., 117, 1364 (1960); R. 
Phys. Rev., 117, 537 (1960); J. I. VETTE: Phys. Rev., 111, 622 (1958). 
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Forward Angle Photoproduction of Single Positive Pions on Hydrogen. 


M. BENEVENTANO, G. FinoccHIARO, R. FINZI 
L. MezzertI, L. PAoLUZZI and C. SCHAERF 


Istituto di Fisica dell Universita - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 25 Luglio 1960) 


Measurements of the differential cross 
section of the reaction 


PAUSE En a 


for photon energies between (600— 
+1000) MeV have been performed by 
several experimenters (13). Systematic 
angular distributions have been given 
by Dixon and WALKER (7). As these 
distributions are of considerable theor- 
etical interest (45), we endeavoured to 
repeat the measurements and to extend 
them toward smaller angles. This is a 


() M. HEINBERG, W. M. MCCLELLAND, 
F. TURKOT, W. M. WooDWARD, R. R. WILSON 
and D. M. Ziroy: Phys. Rev., 110, 1211 (1958). 

() F. P. Dixon and R. L. WALKER: Phys. 
hea Lett, Age 142 (dose LS NO SSD Be ws 
Dixon: Ph. D. Thesis, California Institute of 
Technology (1959). 

C) E. TurKoT: Bull. Am. Phys. Soc., 4, 
22 (1959). 


(*) M. J. MoRAVOSIK: Phys. Rev, 104, 
1451 (1956). 

(©) PR: R. Witson: Phys. Rew, 110, 1212 
(1958). 

(OTREMEPRIEBREISEANPNVS IR COMETA 
174 (1958). 

() J. J. SAKURAI: Phys. Rev. Lett., 41, 
258 (1958). 

(5) A. M. WHETHERELL: Phys. Rev., 115, 


1722 (1959). 


preliminary report on measurements per- 
formed at angles between 10° and 90° 
in the center of mass system (c.m.). 

The experimental arrangement is 
sketched in Fig. 1. The bremsstrahlung 
beam of the Frascati electron synchro- 
tron, suitably collimated and cleaned 
from electron contamination, passed 
through a hydrogen target of the Illinois 
type (°) and was monitored by a Wilson 
type ionization chamber (1). 

Charged particles produced in the 
target were momentum analyzed by a 
magnetic spectrometer. This consisted 
of two equal magnets M, and M, (1:12), 
focusing in the horizontal and defocusing 
in the vertical plane, and two scintilla- 
tion counters S, and S,. 

The target, a vertical cylinder ~50 mm 
in diameter, was situated at the first 
focus of M,. Particles of the right 
momentum and angle emerging from it 


(*) Designed and constructed by the staff 
of the Cryogenic Laboratory at Frascati. 

(9) G. P. MURTAS: Laboratori Nazionali di 
Frascati, internal note n. 18 (1959). 

(11) P. G. Sona: Laboratori Nazionali di 
Frascati, internal note n. 3 (1958). 

(1) R. ToscHI; Laboratori Nazionali d 
Frascati, internal note n. 29 (1960). 
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were focused on counter S,, situated at 
the second focus of M,. 

The counter S, (which was 1 mm 
thick in order to reduce scattering) 


defined the angular resolution A0 of the 


— Experimental layout. The position 
is that corresponding so the 


Fig. 1. 


system, while the counter S, defined the 
vertical acceptance as well as the reso- 
lution in momentum Ap/p,. This system 
was designed so as to reduce the possi- 
bility that a particle scattered by the 
pole pieces be accepted by the counter 
telescope. The whole spectrometer could 
be rotated (12) about a pivot centered on 
the axis of the liquid hydrogen target. 

The floating wire technique has been 
used to determine the mean momentum 
of the spectrometer p, as a function of 
the magnet current setting and the 
acceptance diagram (Ap, A6) for dif- 
ferent target points. 

The results have been checked with 
measurements on pions and electrons 
and can be summarized in the following 
parameters 


Ap/po = + (0.0265 +0.005) 
(for a point source) , 
A@ = (3:08+0.1)°; 


2483 


275 


= (0.0106 em. 


The solid angle has not yet been 
measured directly; a preliminary esti- 


of the spectrometer relative to the y-ray beam 


smallest angle of measurement. 


mate, based on the above parameters 
and on the calculated vertical acceptance 
yields the value AQ ~ 0.47-10-3 sr. 

The counter OC, is a water Cerenkov 
counter, which is sensitive to pions but 
not to protons in the momentum range 
accepted by the spectrometer. Discrimi- 
nation against electrons, whose contri- 
bution is appreciable at small angles, 
was achieved by means of two CO, gas 
Cerenkov counters G, and G, having a 
threshold at a total energy #/mc?~ 15 (18). 

An additional scintillation counter Sy 
was placed behind G,, subtending about 3 
of the solid angle of the main tele- 
scope S,S,C,. 

A particle was identified as a pion 
if it caused a threefold coincidence 
S,S,0,G, within an effective resolving 
time of 27=2-10-8 s. Fourfold coinci- 
dences of the type S,S,0,S,G,G, were 
also registered. Accidentals for both 


(15) M. BENEVENTANO, E. DE AGOSTINO, 
C. A. GALTIERI, B. RISPOLI and A. SERRA: 


Nuovo Cimento, 12, 156 (1959). 
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telescopes were continuously monitored 
by means of suitable delayed coinei- 
dences. 

The lower intensity fourfold telescope 
was less sensitive to spurious events, 
involving low energy electrons and pho- 
tons, which were the main source of 
trouble at the smallest angles. 

Reversed field runs showed that the 
contribution of such spurious events to 
the true pion counting rates from hy- 
drogen was always negligible except for 
the threefold coincidences at 0< 15° 
(c.m.), were it amounted to not more 
than 20%. The empty target back- 
ground was on the average about 7%. 

With this arrangement we have so 
far performed measurements on the an- 
gular distributions at photon energies 
of 600, 700, 800, 900 MeV and on the 
energy dependence of the cross section 
at the two center of mass angles 6=20° 
and 60=90°. The points have been 
chosen with the criterion of comple- 
menting existing data (?). 

The synchrotron was always operated 
at an energy between 75 and 150 MeV 
above the y-ray energy concerned. 

The counting rates, after subtraction 
of empty target background and acci- 
dentals, were converted into centre of 
mass differential cross sections. The 
three fold telescope data were used 
throughout in this calculation, except 
for the two smallest angles (10° and 15° 
in the centre of mass system) where we 
preferred to use the fourfold counting 
rates, because their background of spu- 
rious events was negligible. To do this, 
the fourfold rates were normalized to 
the threefold rates taking into account 
the measured effect of scattering and 
nuclear absorption in the counters. 

The bremsstrahlung spectrum form 
used factor was that measured by Dram- 
BRINI et al. (14) with the same internal 


(15) G. DIAMBRINI, A. S. FIGUERA, B. RI- 
SPOLI and A. SERRA: Laboratori Nazionali di 
Ærascati, internal note n. 20 (1959). 
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target but a somewhat different colli- 
mation. 

No corrections have been applied for 
the pion decay in flight, for the muon 
contamination from x-u decay, and for 
possible pion scatterings from the magnet 
pole pieces. These corrections, however, 
are slowly varying functions of energy 
and angle, and should not change appre- 
ciably the shape of the distributions. 
A very rough estimate gives for the 
highest and lowest net decay correction 
10 and 6% respectively. 


pb 
sr 
(I 3 | 
20 3 Dixon Walker(Dixon thesis 1959) 
È Present work 
154 
104 
54 
(O 20 0 A o 0 o 0 
6,,(Pion CM. angle) 
Fig. 2. — Angular distributions at fixed photon 
energy. The solid curves represent the least 


square fits calculated by Drxon and WALKER 
(FP. Dixon: Thesis, 1959) following the suggestion 
of MORAVCSIK. 


The results are given in Figs. 2, 3, 
4, 5, 6; the errors indicated are purely 
statistical. Points calculated by aver- 
aging over experimental results obtained 
by Drxon and WALKER (2) with various 
arrangements are also shown in the 
angular distribution plots. To facilitate 
the comparison our points have been 
normalized to the best fit curves cal- 
culated by Dixon and Walker, averaging 
over all energies and over the angles 
between 20° and 90° (c.m.). The same 
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photon energy AF, is approximately 
+18 MeV (half width at half height) 
at the lowest energies and angles, and 
increases gradually with angle and energy 
to about +35 MeV at E,,=900 MeV, 


LU 


= 
"AR ] 
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109" 2075590) #40) 9 50) 60170) ‘807 O0 
6,(Pion CMangle) 
Fig. 3. — Angular distributions at fixed photon 
energy. The solid curves represent the least 


squares fits calculated by Dixon and WALKER 
(F. Dixon: Thesis, 1958) following the suggestion 
of MORAVCSIK. 


normalization factor was used for the 


points in Fig. 6. The resolution in 
FR | 
a, 
dy 
20+ + À $ Dixon-Watker(Dixon thesis 1989)! 
Present work 
| | | 
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8, Pion CM angle) 

Fig. 4. — Angular distributions at fixed photon 

energy. The solid curves represent the least 
squares fits calculated by Dixon and WALKER 
(F. Drxon: Thesis, 1959) following the suggestion 

of MORAVOSIK. 
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Fig. 5. — Angular distributions at fixed photon 

energy. The solid curves represent the least 

squares fits calculated by Drxon and WALKER 

{F, Dixon: Thesis, 1959) following the suggestion 
of MoRAVCSIK. 


According to the preliminary esti- 
mates of our solid angle and decay 
correction mentioned above, the values 
represented in Figs. 2, 3, 4 should be 
multiplied by a factor — 0.89 in order 
to get absolute cross sections; but this 
factor may change by as much as +15% 
as a consequence of future better eval- 
uations. 

The analysis of the results is in 
progress. As regards the energy depend- 
ence at fixed center of mass angle, it is 
interesting to note that the maximum 
around 700 MeV appears to be shifted 
towards somewhat smaller energies than 
previously reported. This is not in con- 
tradiction with previous measurements, 
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a ia 


cm. >| 90° 


+ + 
600 700 800 900 
E, (Photon Lab. Energy in MeV) 
Fig. 6. — Differential cross sections as a function 
od photon energy at fixed centre of mass angle 
of the pion. The solid curves are fitted by eye 
to the experimental results. 


but originates essentially from a nar- 
rower spacing of the experimental points. 
Our data seem to indicate also a some- 


what steeper dependence on energy of 
the differential cross sections in the 
region around 700 MeV, an effeet which 
has already been suggested by CORTEL- 
LESSA and REALE (!). Further measu- 
rements with better energy resolution 
are however needed, in our opinion, to: 
establish it with certainty. 


kK * 


Our deepest gratitude goes to the 
memory of our technician Mr. R. Bran- 
CHINI, killed in a road accident while at 
work for this experiment; his intelligent 
and enthusiastic cooperation was inval- 
uable for the preparation and good per- 
formance of the equipment. 

We wish to thank all the Synchrotron 
staff and in particular Dr. G. GHiGo 
for the operation of the machine and 
for their help in maintaining our equip- 
ment and solving some of our technical 
problems. 


(5) G. CORTELLESSA and A. REALE: Nwovo- 
Cimento, to be published. 
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LIBRI RICEVUTI E RECENSIONI 


L. MERTON — Methods of Experi- 
mental Physics. Vol. I: Classical 
Methods, pp. 596, Academic Press, 
New York, 1959, $ 12. 


È questo il primo volume di una 
serie di sei volumi che si propone di 
fornire al fisico sperimentale una rappre- 
sentazione concisa dei più importanti 
metodi sperimentali, in modo che possa 
essere guidato negli argomenti in cui egli 
è meno esperto e trarne ispirazioni e 
consigli utili nel proprio lavoro. Questo 
fine, certo encomiabile, ci sembra che 
sia stato solo in parte raggiunto. Già a 
prima vista la suddivisione della serie 
in sei volumi (il secondo dovrà trattare 
i metodi elettronici, e gli altri quattro 
la fisica molecolare, atomica, nucleare 
e dello stato solido rispettivamente) 
lascia un po’ perplessi. L'introduzione 
ci informa che con il termine « classici » 
si intendono le parti classiche della fisica, 
quali la meccanica, il calore, il suono 
la luce, l’elettricità ed il magnetismo. 
Allora c'è da chiedersi se tutto questo 
sia davvero solo un sesto della fisica 
sperimentale. Per rispondere a questa 
perplessità dobbiamo aspettare di vedere 
i volumi futuri. 

Questo primo volume ci sembra assai 
diseguale, nel senso che in certe parti 
è quasi un ricettario (ad es. vengono 
descritti con cura i metodi per attaccare 
un conduttore ad un dielettrico per misu- 
rare la costante dielettrica) mentre in 
altre parti è di una impressionante bre- 
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vita (ad es. la termometria con metodo 
magnetico è liquidata in due pagine) 
Naturalmente ci sono molte omissioni, 
come ci si può aspettare in un volume 
in fondo così piccolo che voglia trattare 
tanta fisica (ad es. non si parla di tecnica 
delle basse temperature), ma viceversa 
non si capisce la presenza di un lungo 
capitolo assai formale sui principi della 
termodinamica. 

È nostra impressione però che il libro 
possa essere utile come uno studio inter- 
medio di conoscenza dei metodi speri- 
mentali, quando proprio non si sappia 
nulla sull’argomento, o anche si voglia 
fare una rapida ricerca di dati biblio - 
grafici di cui l’opera è abbastanza ricca. 
Ma ci sembra che troppo spesso manchi 
il consiglio dell'autore su questi metodi 
stessi, per cui un profano può spingers 
in qualche direzione senza rendersi conto 
delle difficoltà che incontrerà. 


G. CARERI 


N. G. Parke III — Guide to the 
Literature of Mathematics and 
Physics. Dover Publications, Inc., 
New York, 1958, pp. XVIII-436. 
Prezzo $ 2.64. 


Questo libro costituisce un interes- 
sante tentativo di aiutare le ‘persone, 
già dotate di preparazione scientifica, 
ad iniziare il loro aggiornamento in 
quei campi della matematica e della 
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fisica con i quali non sono in profonda 
consuetudine. 

Per facilitare questo avvicinamento 
l’autore premette alla guida propria- 
mente detta una prima parte in cui 
cerca di dare al lettore un’orientamento 
e un metodo per un proficuo studio per- 
sonale. Non vengono perciò tralasciate 
informazioni minute per una efficiente 
ricerca bibliografica, accanto a chiari- 
ficazioni dei più comuni e importanti 
termini del linguaggio scientifico uniti 
a consigli sulla maniera più convenienti 
di apprendere direttamente dai libri 
quando non ci si possa valere del con- 
tatto diretto con insegnanti. 

La guida propriamente detta è ordi- 
nata alfabeticamente in un centinaio 
di argomenti che corpono in modo 
alquanto completo la matematica e la 
fisica. Ciascuno di questi argomenti, che 
sono ad esempio: equazioni differen- 
ziali, teoria degli errori, fisica nucleare 
ecc., inizia con un breve panorama del- 
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l’autore sulle opere principali dispo- 
nibili a cui segue un lungo elenco di 
opere raggruppate in voci, che sono ad 
esempio: teorie generali, applicazion 
tecniche, ecc., che facilitano al lettore 
la ricerca dei volumi che lo interessano 
specificatamente. 

La scelta è in genere oculata anche: 
se in qualche caso si è forse esagerato 
nel numero includendo opere di discu- 
tibile valore. Completa il volume un 
esauriente indice dei nomi unito a. 
quello analitico assai dettagliato. Si può 
ritenere questa opera assai ben riu- 
scita almeno nei limiti in cui un libro 
può aiutare a ‘supplire alla mancanza. 
di un buon insegnamento personale. 
Esso può risultare assai utile anche a. 
coloro che, nella preparazione di un 
corso universitario, desiderano confon- 
tare il proprio procedere didattico con 
quello seguito da altri in libri e trattati. 
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